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22nd  International  Symposium  on  Ballistics 
Vancouver  BC,  Canada 
14-18  November  2005 


Agenda 


Tuesday.  15  November  2005 

Invited  Presentation:  by  From  Columbia  To  Discovery:  Understanding  the  Impact  Threat  to  the  Space  Shuttle,  by  James  D.  Walker,  Southwest  Research  Institute 

General  Oral  Session  #1 

•  Plasma  Ignition  of  a  30mm  Cannon,  Richard  A.  Beyer,  Andrew  L.  Brant,  Joseph  J.  Colburn,  US  Army  Research  Laboratories 

•  Numerical  Computations  of  Subsonic  and  Supersonic  Flow  Choking  Phenomeana  in  Grid  Finned  Projectiles,  Nicolas  Parise,  SNC  Technologies,  Inc;  Alain 
Dupuis,  Precision  Weapons  Section,  Defense  Research  and  Development  Canada 

•  The  Use  of  Electric  Power  in  Active  Armour  Applications,  Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

•  Prevention  of  Sympathetic  Detonation  Between  Reactive  Armor  Sandwiches,  Andreas  Holzwarth,  Fraunhofer-Institut  Fur  Kurzzeitdynamik 

Terminal  Ballistics  Oral  Session  #7 

•  Bullet  Impact  on  Steel  and  Kevlar®/Steel  Armor  -  Experimental  Data  and  Hydrocode  Modeling  with  Eulerian  and  Lagrangian  Methods,  Dale  S.  Preece, 
Vanessa  S.  Berg,  and  Loyd  R.  Payne,  Sandia  National  Laboratories 

•  Progress  on  the  NDE  Characterization  of  Impact  Damage  in  Armor  Materials,  Joseph  M.  Wells,  JMW  Associates 

•  The  influence  of  sabot  threads  on  the  performance  of  KE  penetrators,  Nick  Lynch  and  John  Stubberfield,  QinetiQ 

•  The  Use  of  Electric  Power  in  Active  Armour  Applications,  Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

•  Prevention  of  Sympathetic  Detonation  Between  Reactive  Armor  Sandwiches,  Andreas  Holzwarth,  Fraunhofer-Institut  Fur  Kurzzeitdynamik 

Exterior  Ballistics  Oral  Session  #1 

•  Micro-Adaptive  Flow  Control  Applied  to  a  Spinning  Projectile,  Dr.  Jubaraj  Sahu,  U.S.  Army  Research  Laboratory 

•  Micro-Adaptive  Flow  Control  and  Nonlinear  Aerodynamics  (Combustion  Gas  Generators),  Dr.  Jubaraj  Sahu  and  Ms.  Karen  Heavey,  U.S.  Army  Research 
Laboratory 

•  Aerodynamic  Characteristics  of  a  Grid  Finned  Projectile  from  Free-Flight  Tests  at  Supersonic  Velocities,  Alain  Dupuis,  DRDC  -  Valcartier;  Claude  Berner, 
French-German  Research  Institute 

•  Recent  Computations  and  Validations  of  Projectile  Unsteady  Aerodynamics,  Roxan  Cayaz,  Eric  Carette,  Giat  Industires;  Remy  Thepot,  Patrick  Champigny, 
Office  National  D'Etudes  et  de  Recherches  Aerospatiales 

•  Research  of  Flight  Characteristics  of  Rod-Type  Projectile  with  Triangular  Cross-section,  Dr.  Wenjun  Yi,  Prof.  Xiaobing  Zhang,  and  Prof.  Jianping  Qian, 
Ballistic  Research  Laboratory  of  China 


Wednesday.  16  November  2005 


Exterior  Ballistics  Oral  Session  #2 

•  Impact  of  Nose-Mounted  Micro-Structures  on  the  Aerodynamics  of  a  Generic  Missile,  Dr.  Daniel  Corriveau,  Defence  R&D  Canada  (DRDC  -  Valcatier) 

•  Analyses  of  Gliding  Control  for  an  Extended-Range  Projectile,  Prof.  Zhongyuan  Wang,  Prof.  Houqian  Xu,  Dr.  Jinguang  Shi,  Dr.  Wenjun  Yi,  Prof.  Shaosong 
Chen,  Ballistic  Research  Laboratory  of  China 

•  Theoretical  Design  for  a  Guided  Supersonic  Projectile,  Pierre  Wey,  Claude  Berner,  Eckhart  Sommer,  Vollcer  Fleck,  Henry  Moulard,  French-German  Research 
Institute  of  Saint-Louis  (ISL) 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #1 

•  Ceramic  Gun  Barrel  Technologies,  Larry  Burton,  Jeff  Swab,  Rob  Carter,  Ryan  Emerson,  U.S.  Army  Research  Laboratory,  Weapons  &  Materials  Research 
Directorate 

•  M865  TID  Improvement  Study,  Kerry  Henry,  Army  Research  Development  &  Engineering  Center;  Jason  W.  Gaines,  General  Dynamics-OTS 

•  Two-Dimensional  Modeling  of  Mortar  Interal  Ballistics,  Clive  R.  Woodley,  David  Finbow,  QinetiQ;  Vladimir  Titarev,  Eleuterio  Toro,  Umeritek  Limited 

•  The  Mechanism  Analysis  of  Interior  Ballistics  of  Serial  Chamber  Gun,  Dr.  Sanjiu  Ying,  Charge  Design  Laboratory  of  China;  Prof.  Xiaobing  Zhang,  Prof. 
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Yaxiong  Yuan,  and  Dr.Yan  Wang,  Ballistic  Research  Laboratory  of  China 
Terminal  Ballistics  Oral  Session  #2 

•  Behind  Armor  Debris  Computations  with  Finite  Elements  and  Meshless  Particles,  Gordon  R.  Johnson  and  Robert  A.  Stryk,  Network  Computing  Services,  Inc. 

•  Experimental  and  Numerical  Study  of  the  Penetration  of  Tungsten  Carbide  Into  Steel  Targets  During  High  Rates  of  Strain,  Eva  K.  Friis,  Nammo  Raufoss  AS; 
Oyvind  Froyland  and  John  F.  Moxnes,  FFI  (Norwegian  Defence  Research  Establishment) 

•  Mine  Neutralisation  with  Small  Calibre  Projetile  Impact,  Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence  Safety  and  Security 

Vulnerability ,  Lethality  and  Wound  Ballistics  Oral  Session 

•  The  Application  of  Critical  Perforation  Analysis  (CPA)  to  Military  Personal  Armour  Research  and  Evaluation,  Catherine  H.  Crawford  and  Philip  Gotts, 
Defence  Clothing  Research  and  Project  Support 

•  Fragment  Patterns  Behind  Concrete  Structures  Caused  by  KE  Projectiles,  S.  Lampert,  Fene  Jeanquartier,  D.  Hoffmann,  and  B.  Lehmann,  Armasuisse 

•  Mine  Neutralisation  with  Small  Calibre  Projetile  Impact,  Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence  Safety  and  Security 


Thursday.  17  November  2005 


Warhead  Mechanisms  Oral  Session  #/: 

•  Soft-Recovery  of  Explosively  Formed  Penetrators,  David  Lambert  and  Matthew  Pope,  Air  Force  Research  Laboratory,  Munitions  Directorate;  Stanley  E. 

Jones  and  Jonathan  Muse,  Aerospace  Engineering  and  Mechanics,  University  of  Alabama 

•  The  Gurney  Velocity:  A  “Constant”  Affected  by  Previously  Unrecognized  Factors,  Joseph  E.  Backofen,  BRIGS  Co. 

•  The  Influence  of  Post  Detonation  Burning  Process  on  Blast  Wave  Parameters  in  Air,  Meir  Mayseless  and  I.  Belsky,  IDF,  Armor  Branch;  E.  Muzychuk,  IMI, 
Central  Laboratory 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #2 

•  3-D  Finite-Element  Gun  Launch  Simulation  of  a  Surrogate  Excalibur  155-mm  Guided  Artillery  Projectile  -  Modeling  Capabilities  and  its  Implications,  M.R. 
Chowdhury  and  A.  Frydman,  US  Army  Research  Laboratory;  J.  Cordes,  L.  Reinhardt  and  D.  Carlucci,  US  Army  ARDEC,  Analysis  and  Evaluation  Division 

•  Caseless  Ammunition  and  Advances  in  the  Characterization  of  High  Ignition  Temperature  Propellant  (HITP),  Paul  Shipley,  AAI  Corporation;  Erin  K. 
Hardmeyer,  US  Army  ARDEC;  and  Ben  Ashcroft,  Alliant  Technical  Systems 

•  Ballistic  Launch  to  Space,  Ed  Schmidt  and  Mark  Bundy,  Army  Research  Laboratory 

•  A  Novel  Launcher  for  Cavitating  Weapons,  Chris  Weiland  and  Pavlos  Vlachos,  Mechanical  Engineering  Department,  Virginia  Polytechnic  Institute  and  State 
University;  and  Jon  Yagla,  Engagement  Systems  Department,  Naval  Surface  Warfare  Center 

Warhead  Mechanisms  Oral  Session  #2 

•  The  Role  of  Rayleigh  Taylor  Instability  in  Shaped  Charge  Jets  Formation  and  Stability,  Dr.  Simcha  Miller,  Mr.  Gershon  Kliminz,  Rafael  Ballistic  Center, 

•  Simulation  of  Cylinder  Expansion  Tests  Using  an  Eulerian  Multiple-Material  Approach,  Laura  Donahue,  R.C.  Ripley,  Martec,  Ltd 

•  Application  of  Powder  Tantalum  Material  for  Explosively  Formed  Penetrator  (EFP)  Warhead,  Richard  Fong,  Mike  Hespos,  William  Ng,  and  Steven  Tang,  US 
Army  ARDEC 

•  Oilwell  Perforators:  Theoretical  Considerations,  Brenden  Grove,  Schlumberger  Reservoir  Completions  Center 

•  The  Study  on  Lethality  Simulation  Method  for  Fragmentation  Warhead,  Yang  Yunbin,  Qu  Ming,  and  Qian  Lixin,  Institute  of  Structual  Mechanics,  China 
Academy  of  Engineering  Physics 


Friday.  18  November  2005 


Terminal  Ballistics  Oral  Session  # 3 : 

•  Performance  Evaluation  of  Multi-Threat  Body  Armour  Systems,  B.  Anctil  and  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau,  M.  Bolduc,  and  D. 
Bourget,  Defence  R&D  Candada  -  Valcartier 

•  The  Residual  Damage  in  CFRP  Composite  After  Ballistic  Impacts  (Experiments  &  Simulations),  Koen  Herlaar,  TNO  Defence,  Security  and  Safety 

•  The  Effect  of  Boundary  Conditions  on  the  Ballistic  Performance  of  Textile  Fabrics,  Colin  R.  Cork,  University  of  Manchester,  School  of  Materials 

General  Oral  Session  #2 

•  Wind  Tunnel  Verification  of  the  Performance  of  a  Smart  Material  Canard  Actuator,  Paul  Weinacht,  William  F.  Drysdale,  Travis  Bogetti,  and  Rod  Don,  US 
Army  Research  Laboratory;  James  T.  Arters,  Jack  R.  Vinson,  Aaron  R.  Hickman,  University  of  Delaware;  Lamar  Auman,  US  Army  Aviation  and  Missile 
RD&E  Center;  Oded  Rabinovitch,  Technion  Israel  Institute  of  Technology 

•  The  Fragmentation  of  Metal  Cylinders  Using  Thermobaric  Explosives,  William  Andrews,  Royal  Military  College  of  Canada;  Michael  Dunning,  Defence  R&D 
Canada  -  Suffield;  and  Kevin  Jaansalu,  Montana  Tech  (University  of  Montana) 

•  A  Novel  Test  Methodology  to  Assess  the  Performance  of  Ballistic  Helmets,  B.  Anctil  and  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau  and  D. 
Bourget,  Defence  R&D  Candada  -  Valcartier 

•  Ballistic  Analysis  of  Bulgarian  Electroslag  Remelted  Dual  Hard  Steel  Armor  Plate,  William  Gooch,  Matthew  Burkins,  and  David  Mackenzie,  US  Army 
Research  Laboratory  Weapons  and  Materials  Research  Directorate;  Stefan  Vodenicharov,  Institute  of  Metal  Science,  Bulgarian  Academy  of  Sciences 
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International  Symposium  on  Ballistics  2005 

International  Symposium  on  Ballistics  2005  is  jointly  organized  and  supported  by  the  National 
Defense  Industrial  Association,  USA  in  conjunction  with  the  International  Ballistics  Committee 

Symposium  Co-Chairman:  William  Flis 
Symposium  Co-Chairman:  Brian  Scott 


PREVIOUS  INTERNATIONAL  SYMPOSIA  ON  BALLISTICS 


1st  Orlando,  Florida,  USA  1974 

2nd  Daytona,  Florida,  USA  1976 

3rd  Karlsruhe,  Germany  1977 

4th  Monterey,  California,  USA  1978 

5th  Toulouse,  France  1980 

6th  Orlando,  Florida,  USA  1981 

7th  The  Hague,  The  Netherlands  1983 

8th  Orlando,  Florida,  USA  1984 

9th  Shrivenham,  UK  1986 

10th  San  Diego,  California,  USA  1987 

11th  Brussels,  Belgium  1989 

12th  San  Antonio,  Texas,  USA  1990 

13th  Stockholm,  Sweden  1992 

14th  Quebec  City,  Canada  1993 

15th  Jerusalem,  Israel  1995 

16th  San  Francisco,  California,  USA  1996 

17th  Midrand,  South  Africa  1998 

18th  San  Antonio,  Texas,  USA  1999 

19th  Interlaken,  Switzerland  2001 

20th  Orlando,  Florida,  USA  2002 

21st  Adelaide,  South  Australia  2004 

22nd  Vancouver,  BC  Canada  2005 


SYMPOSIUM  SCOPE  AND  OBJECTIVES 

The  objective  of  the  22nd  International  Symposium  on  Ballistics  is  to  focus  on  potential  technical 
advances  and  break-throughs  in  the  21st  century  in  the  general  areas  of: 

Interior  Ballistics 

Launch  Dynamics 

Exterior  Ballistics 

Projectile  and  Warhead  Design 

Terminal  Ballistics 

Vulnerability 

Modeling  and  Simulation 
Wound  Ballistics 


Over  200  papers  will  be  presented  by  authors  from  26  countries. 
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SYMPOSIUM  PROGRAM 

Monday.  November  14.  2005 


2:00  pm  -  5:00  pm 

Registration 

5:00  pm  -  6:30  pm 

Reception  in  Exhibit  Area 

Tuesday.  November  15.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

7:00  am  -  6:00  pm 

Exhibits  Open 

8:00  am 

Opening  Remarks 

C.  Samuel  Campagna,  National  Defense  Industrial  Association 

8:10  am 

Welcome  and  Opening  Remarks 

William  Flis,  DE  Technologies,  Inc.  &  Brian  Scott,  US  Army  Research  Laboratory 

8:20  am 

Keynote  Address 

Dr.  Robert  Walker,  Director-General,  Research  and  Development  Programs 
(DGRDP),  Defense  Research  and  Development  Canada 

9:05  am 

Invited  Presentation 

From  Columbia  to  Discovery:  Understanding  the  Impact  Threat  to  the  Space  Shuttle 
James  D.  Walker,  Southwest  Research  Institute 

9:50  am 

Morning  Break 

General  Oral  Session  #1 

Chairpersons:  B.  Janzon  and  J.  Carleone 

10:20  am 

Plasma  Ignition  of  a  30mm  Cannon 

Richard  A.  Beyer,  Andrew  L.  Brant,  Joseph  J.  Colburn,  US  Army  Research 
Laboratories 

10:40  am 

Numerical  Computations  of  Subsonic  and  Supersonic  Flow  Choking  Phenomena  in 
Grid  Finned  Projectiles 

Nicolas  Parise,  SNC  Technologies,  Inc.;  Alain  Dupuis,  Precision  Weapons  Section, 
Defense  Research  and  Development  Canada 

11:00  am 

Multiple  Explosively  Formed  Penetrator  (MEFP)  Warhead  Technologies  for  Mine  and 
Improvised  Explosive  Device  (IED)  Neutralization 

Richard  Fong,  William  Ng,  Steve  Tang,  LaMar  Thompson,  U.S.  Army  Armament 
Research,  Development  and  Engineering  Center 

11:20  am 

The  Use  of  Electric  Power  in  Active  Armour  Applications 

Martin  van  de  Voorde,  R.  Boeschoten,  TNO  Defence,  Security  and  Safety 

11 :40  am 

Prevention  of  Sympathetic  Detonation  between  Reactive  Armor  Sandwiches 

Andreas  Holzwarth,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik 

12:00  pm 

Lunch 
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1 :30  pm 

1 :30  pm 

1 :50  pm 
2:10  pm 

2:30  pm 

2:50  pm 

3:10  pm 
3:40  pm 

3:40  pm 
4:00  pm 

4:20  pm 

4:40  pm 

5:00  pm 

5:20  pm 


3:1 0  pm  Exterior  Ballistics  Poster  Session 

Chairpersons:  Z.  Wang  and  P.A.  Karsten 

Terminal  Ballistics  Oral  Session  #1 
Chairpersons:  E.  Linden  and  C.  Anderson 

Bullet  Impact  on  Steel  and  Kevlar®/Steel  Armor  -  Experimental  Data  and  Hydrocode 
Modeling  with  Eulerian  and  Lagrangian  Methods* 

Dale  S.  Preece,  Vanessa  S.  Berg,  Mathew  A.  Risenmay,  Sandia  National 
Laboratories 

Progress  on  the  NDE  Characterization  of  Impact  Damage  in  Armor  Materials 
Joseph  M.  Wells,  JMW  Associates 

Design,  Analysis,  and  Testing  of  an  Unconfined  Ceramic  Target  to  Induce  Dwell 
Timothy  J.  Holmquist,  Network  Computing  Services,  Inc.;  C.  Anderson,  Jr., 
Southwest  Research  Institute;  Thilo  Behner,  Ernst-Mach-lnstitut 

The  Influence  of  Sabot  Threads  on  the  Performance  of  KE  Penetrators  against 
multiple  plate  targets 

Nick  J.  Lynch,  J.  Stubberfield,  QinetiQ 

Visualization  of  Wave  Propagation  and  Impact  Damage  in  a  Polycrystalline 
Transparent  Ceramic  -  AION 

Elmar  Strassburger,  Fraunhofer  Institut  fur  Kurzzeitdynamik;  Parimal  Patel,  James 
W.  McCauley,  US  Army  Research  Laboratory;  Douglas  W.  Templeton,  US  Army 
TARDEC 

Afternoon  Break 

5:20  pm  Terminal  Ballistics  Poster  Session  #1 
Chairpersons:  A.  Diederen 

Exterior  Ballistics  Oral  Session 
Chairpersons:  W.  Reinecke  and  A.  Dupuis 

Advanced  Time-Accurate  CFD/RBD  Simulations  of  Projectiles  in  Free  Flight 
Jubaraj  Sahu,  US  Army  Research  Laboratory 

Aerodynamic  Characterstics  of  a  Grid  Finned  Projectile  from  Free-Flight  Tests  at 
Supersonic  Velocities 

Alain  Dupuis,  DRDC  -  Valcartier;  Claude  Berner,  French-German  Research 
Institute 

Recent  Computations  and  Validations  of  Projectile  Unsteady  Aerodynamics 
Roxan  Cayaz,  Eric  Carette,  Giat  Industries;  Remy  Thepot,  Patrick  Champigny, 
Office  National  d’ Etudes  et  de  Recherches  Aerospatiales 

The  Derivation  of  Spin  Stabilised  Projectile  Yaw  Rates  and  Ballistic  Model 
Coefficients  Using  Conventional  CW  Doppler  Radar  Systems 

John  Tate,  FLEET 

Research  of  Flight  Characteristics  of  Rod-Type  Projectile  with  Triangular 
Cross-Section 

Wenjun  Yi,  Xiaobing  Zhang,  Jianping  Qian,  Ballistic  Research  Laboratory  of 
China,  Nanjing  University  of  Science  &  Technology 

Adjourn  for  the  Day 
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Wednesday.  November  16.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

7:00  am  -  5:00  pm 

Exhibits  Open 

8:00  am 

Administrative  Remarks 

8:10  am  -  9:50  am 

Terminal  Ballistics  Poster  Session  #2 

Chairpersons:  J.  Riegel  and  E.  Hirsch 

Exterior  Ballistics  Oral  Session  #2 

Chairpersons:  P.  Nel  and  E.  Schmidt 

8:10  am 

Impact  of  Nose-Mounted  Micro-Structures  on  the  Aerodynamics  of  a  Generic  Missile 
Daniel  Corriveau,  Defence  R&D  Canada  (DRDC  -  Valcatier) 

8:30  am 

Ballistic  Simulations  and  Wind  Tunnel  Testing  of  120  mm  Mortar  Bomb  Tail  Fin 
Geometries  -  In  Search  for  Extra  Range 

Jukka  Tiainen,  Ari  Makkonen,  Patria  Weapon  Systems  Oy;  Mikko  Korhonen, 
Timo  Sailaranta,  TKK/Laboratory  of  Aerodynamics 

8:50  am 

Bringing  Solid  Fuel  Ramjet  Projectiles  Closer  to  Application  -  An  Overview  of  the 
TNO/RWMS  Technology  Demonstration  Programme 

Ronald  G.  Veraar,  TNO  Defence,  Security  and  Safety  Research  Group  Rocket 
Technology;  Guido  Giusti,  Rheinmetall  Waffe  Munition  Schweiz  AG 

9:10  am 

Analysis  of  Gliding  Control  for  an  Extended-Range  Projectile 

Zhongyuan  Wang,  Houqian  Xu,  Jinguang  Shi,  Wenjun  Yi,  Shaosong  Chen, 
Ballistic  Research  Laboratory  of  China,  Nanjing  University  of  Science  &  Technology 

9:30  am 

Theoretical  Design  for  a  Guided  Supersonic  Projectile 

Pierre  Wey,  Claude  Berner,  Eckhart  Sommer,  Volker  Fleck,  Henry  Moulard, 

French-German  Research  Institute  of  Saint-Louis  (ISL) 

9:50  am 

Morning  Break 

10:20  am  - 12:00  pm 

Warhead  Mechanisms  Poster  Session 

Chairpersons:  R.  Fong  and  F.  Mostert 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #1 

Chairpersons:  C.  Candland  and  C.  Woodley 

10:20  am 

Ceramic  Gun  Barrel  Technology 

Lawrence  W.  Burton,  Jeffrey  J.  Swab,  Ryan  Emerson,  Robert  Carter,  US  Army 
Research  Laboratory,  Weapons  &  Materials  Research  Directorate 

10:40  am 

M865E3  Cold  Target  Impact  Dispersion  Study 

Kerry  Henry,  Army  Research  Development  &  Engineering  Center;  Jason  W. 
Gaines,  General  Dynamics-OTS 

11:00  am 

An  Alternative  Technique  to  Evaluate  and  Characterize  Pressure  Waves  in  Large 
Calibre  Guns 

Victor  Schabort,  Denel  Land  Systems  Western  Cape 

11:20  am 

Two-Dimensional  Modelling  of  Mortar  Internal  Ballistics 

Clive  R.  Woodley,  David  Finbow,  QinetiQ;  Vladimir  Titarev,  Eleuterio  Toro, 
Numeritek  Limited 
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1 1 :40  am 

12:00  pm 

1 :30  pm 
1 :50  pm 

2:10  pm 
2:30  pm 

2:50  pm 
3:10  pm 

3:40  pm 

4:00  pm 

4:20  pm 

4:40  pm 


The  Mechanism  Analysis  of  Interior  Ballistics  of  Serial  Chamber  Gun 
Sanjiu  Ying,  Charge  Design  Laboratory  of  China,  Nanjing  University  of  Science  & 
Technology;  Xiaobing  Zhang,  Qaxiong  Yuan,  Yan  Wang,  Ballistic  Research 
Laboratory  of  China,  Nanjing  University  of  Science  &  Technology 

Lunch 


Terminal  Ballistics  Oral  Session  #2 
Chairpersons:  M.  Mayseless  and  T.  Holmquist 

Behind  Armor  Debris  Computations  with  Finite  Elements  and  Meshless  Particles 
Gordon  R.  Johnson,  Robert  A.  Stryk,  Network  Computing  Services,  Inc. 

Expirmental  and  Numerical  Study  of  the  Penetration  of  Tungsten  Carbide  into  Steel 
Targets  During  High  Rates  of  Strain 

Eva  K.  Friis,  Nammo  Raufoss  AS,  Oyvind  Froyland,  John  F.  Moxnes,  FFI 
(Norwegian  Defence  Researdh  Establishment 

Fragmentation  Behavior  of  Tungsten  Alloy  Cubes  on  Normal  Aluminum  Plate  Targets 
Karl  Weber,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach  Institut 

The  Failure  Kinetics  of  High  Density  DEDF  Glass  Against  Rod  Impact  at  Velocities 
From  0.4  to  2.5  km/s 

Thilo  Behner,  V.  Hohler,  M.  Moll,  Fraunhofer  Institut  fur  Kurzzeitdynamik  (Ernst- 
Mach  Institut);  Ch.  E.  Anderson  Jr.,  Southwest  Research  Institute;  D.  L.  Orphal, 
International  Research  Associates,  Inc.;  D.  W.  Templeton,  US  Army 
RDECOM-TACOM 

Mine  Neutralisation  with  Small  Calibre  Projectile  Impact 

Mark  Dijkstra,  J.H.  Meulman,  TNO  Defence,  Safety  and  Security 

Afternoon  Break 

Vulnerability,  Lethality  and  Wound  Ballistics  Oral  Session 
Chairpersons:  R.  Vaziri,  A.  Persson 

The  Application  of  Critical  Perforation  Analysis  (CPA)  to  Military  Personal  Armour 
Research  and  Evaluation 

Catherine  H.  Crawford,  Philip  Gotts,  Defence  Clothing  Research  and  Project 
Support 

An  Efficient  Mechanistic  Approach  to  Modelling  the  Ballistic  Response  of  Multi-Layer 
Fabrics 

Ali  Shahkarami,  Reza  Vaziri,  Anounsh  Poursartip,  Composites  Group, 
Departments  of  Civil  Engineering  and  Materials  Engineering  The  University  of  British 
Columbia;  Navin  Tajani,  DuPont  Advanced  Fibers  Systems 

Pencilling  -  A  Novel  Behind  Armour  Blunt  Trauma  Injury 

ElunedA.  Lewis,  Defence  Clothing  Research  and  Project  Support;  Ian  Horsfall, 
Celia  Watson,  Engineering  Systems  Department,  Royal  Military  College  of  Science, 
Cranfield  University 

Scaling  the  Dynamic  Response  of  Armored  Vehicle’s  Floor  Subjected  to  a  Large 
Buried  Charge 

Avidov  Neuberger,  MOD,  Tank  Program  Management;  S.  Peles,  IMI,  Central 
Laboratory  Division;  D.  Rittel,  Technion,  Israel  Institute  of  Technology,  Faculty  of 
Mechanical  Engineering 


22nd  International  Symposium  on  Ballistics 


5:00  pm 

Fragment  Patterns  Behind  Concrete  Structures  Caused  by  KE  Projectiles 

Rene  Jeanquartier,  D.  Hoffmann,  S.  Lampert,  B.  Lehmann,  Armasuisse 

5:20  pm 

Adjourn  for  the  Day 

Thursday.  November  17.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

8:00  am  - 11 :00  am 

Exhibits  Open 

8:00  am 

Adminstrative  Remarks 

8:10  am  -  9:50  am 

Interior  Ballistics/Launch  Dynamics  Poster  Session 

Chairpersons:  C.  Woodley  and  C.  Candlant 

Warhead  Mechanisms  Oral  Session  #1 

Chairpersons:  M.  Murphy  and  P.Y.  Chanteret 

8:10  am 

Soft-Recovery  of  Explosively  Formed  Penetrators 

David  E.  Lambert,  Matthew  Pope,  Air  Force  Research  Laboratory,  Munitions 
Directorate;  Stanley  Jones,  Jonathan  Muse,  University  of  Alabama,  Aerospace 
Engineering  and  Mechanics 

8:30  am 

The  Gurney  Velocity:  A  “Constant”  Affected  by  Previously  Unrecognized  Factors 

Joseph  E.  Backofen,  BRIGS  Co. 

8:50  am 

Influence  of  Post  Detonation  Burning  Process  on  Blast  Wave  Parameters  in  Air 

Meir  Mayseless,  E.  Muzychuk,  IDF,  Mil.;  M.  Mayseless,  1.  Belsky,  IMI,  Central 
Laboratory 

9:10  am 

Steerable  Fragment  Masses 

Manfred  Held,  TDW/EADS 

9:30  am 

Penetration  Performances  of  Tungsten-Copper  Shaped  Charge  Liner 

Seong  Lee,  Eun  Pyo  Kim,  Youngmoo  Kim,  Sung  Ho  Lee,  Moon-Hee  Hong, 
Joon-Woong  Noh,  Agency  for  Defense  Development 

9:50  am 

Morning  Break 

10:20  am  -  12:00  pm 

Vulnerability/Lethality/Wound  Ballistics  Poster  Session 

Chairpersons:  W.  Gooch 

Interior  Ballistics/Launch  Dynamics  Oral  Session  #2 

Chairpersons:  B.  Burns  and  R.  Cayzac 

10:20  am 

3-D  Finite-Element  Gun  Launch  Simulation  of  a  Surrogate  Excalibur  155-mm  Guided 
Artillery  Projectile  -  Modeling  Capabilities  and  its  Implications 

M.R.  Chowdhury,  A.  Frydman,  US  Army  Research  Laboratory;  J.  Cordes,  L. 
Reinhardt,  D.  Carlucci,  US  Army  ARDEC,  Analysis  and  Evaluation  Division 

10:40  am 

Method  of  Calculating  Initial  Firing  Data  of  Artillery  Laser  Terminal-Guidance  Weapon 
System 

Feipeng  Zeng,  Liren  Liu,  Faculty  of  Artillery  Command,  Nanjing  Artillery  Academy 

22nd  International  Symposium  on  Ballistics 


11:00  am 

Caseless  Ammunition  &  Advances  in  the  Characterization  of  High  Ignition 
Temperature  Propellant 

Patricia  M.  O’Reilly,  Erin  Hardmeyer,  Chad  Sensenig,  US  Army  ARDEC;  Ben 
Ashcroft,  Alliant  Techsystems;  Dave  Cleveland,  The  Johns  Hopkins  University, 
OApplied  Physics  Laboratory;  Bo  Engel,  Paul  Shipley,  AAI  Corporation 

11:20  am 

Ballistic  Launch  to  Space 

Edward  Schmidt,  M.  Bundy,  US  Army  Research  Laboratory 

1 1 :40  am 

A  Novel  Launcher  for  Cavitating  Weapons 

Chris  J.  Weiland,  Pavlos  P  Vlachos,  Dept  of  Mechanical  Engineering 

Virginia  Tech;  Jon  J  Yagla,  Mechanical  Engineer  Engagement  Systems  Department 
Naval  Surface  Warfare  Center 

12:00  pm 

Lunch 

Warhead  Mechanisms  Oral  Session  #2 

Chairpersons:  R.  Brown  and  M.  Held 

1 :30  pm 

The  Role  of  Rayleigh  Taylor  Instability  in  Shaped  Charge  Jets  Formation  and  Stability 
Simcha  Miller,  Gershon  Kliminz,  Rafael  Ballistic  Center 

1 :50  pm 

Simulation  of  Cylinder  Expansion  Tests  Using  an  Eulerian  Multiple-Material  Approach 

Laura  K.  Donahue,  R.C.  Ripley,  Martec,  Ltd. 

2:10  pm 

Application  of  Powder  Tantalum  Material  for  Explosively  Formed  Penetrator  Warhead 
Richard  Fong,  William  Ng,  Steven  Tang,  Michael  Hespos,  US  Army  Armament 
Research,  Development  and  Engineering  Center 

2:30  pm 

Oilwell  Perforators:  Theoretical  Considerations 

Brenden  M.  Grove,  Schlumberger  Reservoir  Completions  Center 

2:50  pm 

Planar  Cutting  Jets  from  Shaped  Charges 

Geoffery  EB  Tan,  T.K.  Lam,  Y.K.  Tham,  DSO  National  Laboratories 

3:10  pm 

The  Study  on  Lethality  Simulation  Method  for  Fragment  Warhead 

Yang  Yunbin,  Qu  Ming,  Qian  Lixin,  Institute  of  Structural  Mechanics,  China 
Academy  of  Engineering  Physics 

3:30  pm 

Adjourn  for  the  Day 

4:00  pm  -  5:30  pm 

Reception 

Friday.  November  18.  2005 


7:00  am 

Continental  Breakfast  and  Registration 

8:00  am 

Adminstrative  Remarks 

Terminal  Ballistics  Oral  Session  #3 

Chairpersons:  1.  Cullis  and  D.  Nandlall 

8:10  am 

Performance  Evaluation  of  Multi  threat  Body  Armour  Systems 

B.  Anctil,  M.  Keown,  Biokinetics  and  Associates  Ltd.;  G.  Pageau,  M.  Bolduc  D. 
Bourget,  Defence  R&D  Canada  -  Valcartier 

22nd  International  Symposium  on  Ballistics 


8:30  am 

8:50  am 

9:10  am 

9:30  am 

9:50  am 

10:20  am 

10:40  am 

11 :00  am 

11:20  am 

1 1 :40  am 

12:00  pm 

12:15  pm 
12:25  pm 


Finite  Element  Simulations  and  Experiments  to  Determine  the  Residual  Damage  of  a 
CFRP  Composite  Material  After  Ballistic  Impacts 

Koen  Herlaar,  M.  Van  der  Jagt-Deutekom,  TNO  Defence,  Security  and  Safety 

The  Effect  of  Boundary  Conditions  on  the  Ballistic  Performance  of  Textile  Fabrics 
Colin  R.  Cork,  University  of  Manchester,  School  of  Materials 

Terminal  Ballistic  Effects  of  Low  Density  Materials  Used  as  Confinement  Plates  for 
Explosive  Reactive  Armours 

Hanspeter  Kaufmann,  RUAG  Land  Systems;  Andre  Koch,  Armasuisse 

Quantification  of  the  Effect  of  Using  the  Johnson-Cook  Damage  Model  in  Numerical 
Simulations  of  Penetration  and  Perforation 

Charles  E.  Anderson  Jr.,  T.  R.  Sharron,  Southwest  Research  Institute;  Timothy  J. 
Holmquist,  Network  Computing  Services,  Inc. 

Morning  Break 


General  Oral  Session  #2 
Chairpersons:  V.  Sanchez-Galvez  and  P.  Cuniff 

Comparisons  of  Internal  Ballistics  Simulations  of  the  AGARD  Gun 
Clive  R.  Woodley,  QinetiQ;  Alain  Carriere,  Patrice  Franco,  Dieter  Hensel,  Julien 
Nussbaum,  Institut  Franco-Allemand  de  Recherches  de  Saint-Louis  (ISL);  Tatjana 
Groge,  Ernst-Mach-lnstitut  (EMI);  Stefan  Kelzenberg,  Fraunhofer-lnstitut  fur 
Chemische  Technologie  (ICT),  Baptiste  Longuet,  DGA/DCE/ETBSr3 

Wind  Tunnel  Verification  of  the  Performance  of  a  Smart  Material  Canard  Actuator 
Paul  Weinacht,  William  F.  Drysdale,  Travis  Bogetti,  Rod  Don,  US  Army  Research 
Laboratory;  James  T.  Arters,  Jack  R.  Vinson,  Aaron  R.  Hickman,  University  of 
Delaware;  Lamar  Auman,  US  Army  Aviation  and  Missile  RD&E  Center;  Oded 
Rabinovitch,  Technion  Israel  Institute  of  Technology 

Fragmentation  of  Metal  Cylinders  Using  Thermobaric  Explosives 
M.R.  Dunning,  Defence  Research  and  Development  -  Suffield,  IMS.  Andrews, 
Department  of  Chemistry  and  Chemical  Engineering,  Royal  Military  College  of 
Canada;  K.M.  Jaansalu,  Department  of  Metallurgical  and  Materials  Engineering,  The 
University  of  Montana 

A  Novel  Test  Methodology  to  Assess  the  Performance  Ballistic  Helmets 
B.  Anctil,  M.  Keown,  Biokinetics  and  Associates  Ltd.;  D.  Bourget,  G.  Pageau, 
Defence  R&D 

Ballistic  Analysis  of  Bulgarian  Electroslag  Remelted  Dual  Hard  Steel  Armor  Plate 
William  Gooch,  Matthew  Burkins  and  David  Mackenzie,  US  Army  Research 
Laboratory,  Weapons  and  Materials  Research  Directorate;  Stefan  Vodenicharov, 
Institute  of  Metal  Science,  Bulgarian  Academy  of  Sciences 

Presentation  of  Awards 

The  Rosalind  and  Pei  Chi  Chou  Award  for  Young  Authors 
The  Neil  Griffiths  Memorial  Award 
The  Louis  and  Edith  Zernow  Award 

Invitation  to  the  23rd  International  Symposium  on  Ballistics,  Tarragona,  Spain,  2007 
Closing 

William  Flis,  DE  Technologies,  Inc.  &  Brian  Scott,  US  Army  Research  Laboratory 


POSTER  SESSIONS  START  HERE 

Exterior  Ballistics  Poster  Session 
1 :30  pm  -  3:10  pm  Tuesday,  November  15 

1913  Fractional  Calculus  for  Design  of  Aerodynamic  Missile’s  Autopilot  and  Digital  Realization 

Bangchu  Zhang,  Chenming  Li,  Zipeng  Han,  Zou  Yun,  Fuming  Xu,  Ballistic  Research  Laboratory  of  China,  Nanjing 

University  of  Science  &  Technology 

1915  The  Simulation  of  Rocket  Trajectory  in  Simulink 
Xin  Changfan,  Nanjing  University  of  Science  &  Technology 

1924  Establishing  a  Pitch  Damping  Test  Capability  at  CSIR  Defencetek 
Fabrizio  Dionisio,  CSIR,  Defencetek 

1944  The  Investigation  About  Using  Different  Guidance  Laws  on  Improving  Impact  Point  Deviation  of  a  Rocket 
Handong  Zhao,  Fang  Wang,  Qingshang  Liu,  Key  Laboratory  of  Instrumentation  and  Dynamic  Measurement,  North 
University  of  China 

1951  Numerical  Integration  Method  Based  on  4th  Largrange  Polynomial  of  Strap-Down  INS  System 
Guoguang  Chen,  Xiaoli  Tian,  Changfan  Xin,  Yaqi  Bao,  North  University  of  China 

1952  Research  on  Real  Time  Trajectory  Measure  Device  of  Range 
Changfan  Xin,  Guoguang  Chen,  Xiaoli  Tian,  North  University  of  China 

1953  Research  on  Attitude  Control  Strategy  of  Glide  Range  Extend  Rocket 
Xiaoli  Tian,  Guoguang  Chen,  Changfan  Xin,  North  University  of  China 

1954  Optimal  Algorithm  of  Glide  Range  Extend  Rocket’s  Trajectory 
Guoguang  Chen,  Xiaoli  Tian,  Changfan  Xin,  North  University  of  China 

1978  Practical  Propulsion  by  Directed  Energetic  Processes 

Joseph  P.  Backofen,  BRIGS  Co. 

2014  Investigating  the  Method  of  Obtaining  Ammunition  Roll  Attitude  by  Detecting  the  Geomagnetic  Vector 
Hongsong  Cao,  Guoguang  Chen,  Department  of  Mechatronics  Engineering,  North  University  of  China 

2065  External  Ballistic  Trajectory  Computations  for  Direct/Indirect  Fire  Weapon  Systems 
David  J.  Norton,  General  Dynamics  Canada 

2087  The  Influence  of  Laser  Rangefinder  Parameters  on  the  Hit  Probability  in  Direct  Tank  Fire 
Vladimir  Cech,  OPROX,  Inc.,  Jiri  Jevicky,  Department  of  Mathematics,  University  of  Defense 

2123  Flight  Dynamics  Modeling  and  Experiment  for  Composite  Concepts.  Application  to  Ribbon  Aerodynamic 
Stabilization 

Christopher  Grignon,  S.Heddadj,  Giat  Industries 

21 28  Onboard  Measurements  with  Magnetic  Sensors:  Determination  of  the  Attitude  and  the  trajectory  Position 

V.  Fleck,  E.  Sommer,  S.  Changey,  French-German  Research  Institute  (ISL);  D.  Beauvois,  Ecole  Superirure  d’Electrecite 

(Supelec) 

2133  Aerodynamic  Characteristics  of  a  Long  Range  Spinning  Artillery  Shell  Obtained  from  3D  Magnetic  Sensors 
V.  Fleck,  E.  Sommer,  C.  Berner,  French-German  Research  Institute  (ISL);  A.  Dupuis,  DRDC 

2143  Experimental  Testing  and  Numerical  Simulation  of  Separation  Disturbances  for  Two-Stage  Kinetic  Energy  Missiles 
Nicolas  Parise,  SNC  Technologies,  Inc.;  Richard  Lestage,  Francoise  Lesage,  Precision  Weapons  Section,  Defense 
Research  and  Development  Canada  Valcartier 

2147  Numerical  Study  on  the  Base  Drag  Characteristics  of  a  Base  Bleed  Projectile  with  a  Central  Propulsive  Jet 
Chang-Kee  Kim,  Agency  for  Defense  Development;  J.Y.  Choi,  Pusan  National  University,  Department  Aerospace 
Engineering 


2169  Solid  Fuel  Ramjet  (SFRJ)  Propulsion  for  Artillery  Projectile  Applications  -  Dynamic  Testing  Progress 

Anton  Stockenstrom,  Dynax 

3010  Pitch  and  Bending  During  In-Flight  Extension 

W.  G.  Reinecke.  Institute  for  Advanced  Technology;  M.  G.  Miller,  Physical  Sciences,  Inc. 

4004  Improvements  in  Aerodynamic  Design  for  KE  Less-Lethal  Projectiles 
Jamie  H.  Cuadros,  Arts  &  Engineering 


Terminal  Ballistics  Poster  Session  #1 
3:40  pm  -  5:20  pm  Tuesday,  November  15 

1 001  Numerical  Simulations  of  Silicon  Carbide  Tiles  Impacted  by  Tungsten  Carbide  Spheres 

Constantine  G.  Fountzoulas,  Jerry  C.  LaSalvia,  Bryan  A.  Cheeseman,  Weapons  and  Materials  Research  Directorate; 
Michael  J.  Normandia,  Ceradyne,  Inc. 

1007  Shock  Mitigation  for  Blast  Protection  Using  Hertzian  Tapered  Chains 

Robert  Doney,  US  Army  Research  Laboratory;  Su raj it  Sen,  Department  of  Physics,  State  University  of  New  York  at 
Buffalo 

1011  A  Predictive  Model  for  the  Dwell/Penetration  Transition  Phenomenon 
Jerry  C.  LaSalvia,  US  Army  Research  Laboratory 

1 01 2  Effect  of  Ceramic  Thickness  on  the  Dwell/Penetration  Transition  Phenomenon 
Jerry  C.  LaSalvia,  US  Army  Research  Laboratory 

1 01 4  The  Development  of  Hybridized  Thermoplastic-Based  Structural  Materials  with  Applications  to  Ballistic  Helmets 
Shawn  Walsh,  Brian  R.  Scott,  David  M.  Spagnuoio,  AMSRD-ARL-WM-MB 

1 01 5  Time  Resolved  Observation  of  the  Deformation  and  Surface  Strain  of  a  Textile  Fabric  Subject  to  Ballistic  Impact 
Brian  Scott,  Peter  Dehmer,  US  Army  Research  Laboratory;  Timothy  Schmidt,  Trilion  Quality  Systems 

1 01 6  Analytic  Design  Trends  of  Fabric  Armor 

Brian  Scott,  Chian-Fong  Yen,  US  Army  Research  Laboratory 

1018  High-Speed  Photographic  Study  of  Wave  and  Fracture  Propagation  in  Fused  Silica 

Elmar  Strassburger,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach-lnstitut  (EMI);  Parimal  Patel,  James  W. 

McCauley,  US  Army  Research  Laboratory;  Douglas  W.  Templeton,  US  Army  TARDEC 

1022  Low  Velocity  Ballistic  Properties  of  Shear  Thickening  Fluid  (STF)-Fabric  Composites 

M.  J.  Decker,  R.  G.  Egres,  N.  J.  Wagner,  University  of  Delaware,  Dept,  of  Chemical  Engineering  and  Center  for 
Composite  Materials;  E.  D.  Wetzel,  U.S.  Army  Research  Laboratory 

1 023  An  Approximate  Solution  of  the  Long-Rod  Penetration  Equations 
William  Walters,  Cyril  Williams,  ARL,  Terminal  Effects  Division 

1901  Tubular  Projectile  Interaction  with  Stationary  and  Moving  Oblique  Plates 

Olof  Andersson,  Swedish  Defence  Research  Agency  (FOI),  Weapons  and  Protection  Division 

1911  A  Study  on  the  Moving  Features  of  Double-Layer  Explosive  Reactive  Armor  with  Definite  Angle  by  Numerical 
Simulation  and  Experiments 

Zhengxiang  Huang,  Xiant eng  Zhang,  Gang  Li,  School  of  Mechanical  Engineering,  Nanjing  University  of  Science  & 
Technology 

1927  Mechanics  of  Structural  Design  of  EPW  Warhead 

X.W.  Chen,  Institute  of  Structural  Mechanics,  China  Academy  of  Engineering  Physics 

1928  Armour  Qualification  Utilizing  Maximum  Likelihood  Ballistic  Limit  Calculation 
Moshe  Ravid,  Shlomo  Galperin ;  Rimat  Advanced  Technologies,  Ltd. 


1934  Perforation  of  Concrete  Targets  by  an  Eroding  Tungsten-Alloy  Rod 

Stephan  Lampert,  Rene  Jeanquartier,  Armasuisse 

1955  Ballistic  Properties  of  Single-Melt  Titanium-6Aluminum-4Vanadium  Alloy  Plate 

Brij  J.  Roopchand,  US  Army  Tank-Automotive  and  Armament  Command,  Armament  Research,  Development,  & 
Engineering  Center 

1987  Preliminary  Investigations  of  Potential  Light  Weight  Metallic  Armour  Applications 
Martin  van  de  Voorde,  A.M.  Dierderen,  K.  Herlaar,  TNO  Defence,  Safety  and  Security 

2001  Oblique  Warhead  Penetration  and  Perforation  of  Multi-Layered  Metallic  Targets 

Yongxiang  Dong,  Feng  Shunshan,  Wang  Fang,  State  Key  Laboratory  of  Explosion  Science  &  Technology,  Beijing 
Institute  of  Technology 

2019  Influence  of  Projectile  Material  on  Yawed  Long  Rod  Projectiles  Penetrating  Oblique  Plates 
Ewa  Liden,  Swedish  Defence  Research  Agency  (FOI),  Weapons  and  Protection  Division 

2035  Advanced  Aliphatic  Polyurthane  Resins  for  High  Durability  and  Superior  Ballistic  Performance  Laminated  Glass 
Francisco  Folgar,  INTER  Materials,  LLC 

2037  Impact  and  Penetration  of  B4C  Ceramic,  Aluminum,  and  Berlyllium  by  Depleted  Uranium  Rods  at  2.0  KM/S 
Scott  A.  Mullin,  James  D.  Walker,  Carl  E.  Weiss,  Southwest  Research  Institute;  Paul  O.  Leslie,  Los  Alamos  National 
Laboratories 

21 22  A  Comparison  of  Some  Analytical  and  Empirical  Models  for  Kinetic  Energy  Penetration  of  Semi-Infinite  and  Finite 
Thickness  Steel  Targets 
Nick  J.  Lynch,  J  T Mills,  QinetiQ 

2181  Computed  Tomography  of  High-Speed  Events 

Karsten  Michael,  Philip  Heiberg,  Fraunfoher  Institute  for  High  Speed  Dynamics,  Ernst-Mach-lnstitut 

2186  Characterization  of  Behind-Armor  Debris  Particles  from  Tungsten  Penetrators 
Brad  A.  Pedersen,  S.  Bless,  Institute  for  Advanced  Technology 

Terminal  Ballistics  Poster  Session  #2 
8:10  am  -  9:50  am  Wednesday,  November  16 

2050  On  the  Critical  Thickness  of  Ceramic  to  Shatter  WC-Co  Bullet  Cores 

Paul  J.  Hazell,  Engineering  Systems  Department,  Cranfield  University,  Royal  Military  College  of  Science;  C.  J.  Roberson, 
Advanced  Defence  Materials  Limited 

2060  The  Effect  of  Spaced  Armour  on  the  Penetration  of  Shaped  Charge  Warheads 
James  D.  Shattock,  Cranfield  University 

2072  Modeling  Impact  and  Penetration  Using  a  Deterministic  and  Probabilistic  Design  Tool 

David  S.  Riha,  Jason  B.  Plemlng,  Ben  H.  Tucker,  Scott  A.  Mullin,  James  D.  Walker,  Carl  E.  Weiss,  Southwest 
Research  Institute;  Edward  A.  Rodriguez,  Paul  O.  Leslie,  Los  Alamos  National  Laboratories 

21 06  On  the  Ballistic  Efficiency  of  the  Three  Layered  Metallic  Targets 

Stanislav  Role,  Military  Technical  Institute  of  Protection;  Jaroslav  Buchar,  Mendel  University;  Giovanni  Cozzani,  OTO 
MELARA  S.p.A  Vojtech  Hruby,  University  of  Defence 

2107  Effect  of  the  Temperature  on  the  Ballistic  Efficiency  of  Plates  Made  From  Cast  Iron 
Stanislav  Role,  Military  Technical  Institute  of  Protection;  Jaroslav  Buchar,  Mendel  University 

21 1 3  Displacement  Device  to  Measure  the  Acceleration  of  the  Bulge  of  RHA  Plates  Under  Anti-Tank  Mine  Blast 
Manfred  Held,  TDW/EADS;  Peter  Heeger,  WTD;  Josef  Kiermeir,  CONDAT 


21 1 6  Defeating  Mechanisms  of  Explosive  Reactive  Armour  Sandwiches 
Manfred  Held,  TDW/LFK/EADS 

2121  Comparisons  of  Unitary  and  Jacketed  Rod  Penetration  into  Semi-Infinite  and  Oblique  Plate  Targets  at  System 
Equivalent  Velocities 

John  Stubberfield,  NJ  Lynch,  QinetiQ;  /  Wallis,  QinetiQ  Farnborough 

2126  Finite  Element  Simulations  and  Experiments  of  Ballistic  Impacts  on  High  Performance  PE  Composite  Material 
Koen  Herlaar,  M.  Van  der  Jagt-Deutekom,  TNO  Defense,  Security  and  Safety 

2129  The  Use  of  Foam  Structures  in  Armoured  Vehicle  Protection  Against  Landmines 

David  A.  Cendon,  Vincente  Sanchez-Galvez,  Francisci  Galvez,  Alejandro  Enfedaque,  Departamento  de  Ciencia  de 
Materiales,  E.T.S.I.  Caminos,  Canales  y  Puertos,  Universidad  Politecnica  de  Madrid,  Spain 

2136  The  Numerical  Simulation  of  the  Impact  of  an  Aluminum  Cylinder  into  a  Steel  Cone 
Izak  M.  Snyman,  Defencetek  Landwards  Programme,  CSIR 

2139  Characterization  of  Al  6061 -T6  using  Split  Hopkinson  Bar  Tests  and  Numerical  Simulations 

Amal  Bouamoul,  Manon  Bolduc,  DRDC  -  Valcartier 

2145  Finite  Element  Modeling  of  Light  Armoured  Vehicles  (LAV)  Welds  Heat  Affected  Zones  Sublected  to  an 
Anti-Vehicular  (AV)  Blast  Landmine  Loading:  A  Summary  of  the  Numerical  Model  and  Experiment 
Patrice  Gaudreault,  Defence  Research  &  Development  Canada;  Amal  Bouamoul,  Robert  Durocher,  Benoit  St-Jean, 
DRDC  Valcartier 

2151  Designer  Projectiles  by  Density  Variation:  Towards  the  Nano-Projectile 

John  P.  Curtis,  QinetiQ 

2160  Numerical  Simulation  for  the  Front  Section  Effect  of  Missile  Warhead  on  the  Target  Perforation 

Ho  Soo  Kim,  Ki-Sun  Yeom,  Seong  Shik  Kim,  Agency  for  Defense  Development  (ADD);  Larry  Sotsky,  US  Army  ARDEC 

2180  The  Penetration  Process  of  Projectiles  into  Long  Bars  in  the  Axial  Direction 

Dan  Yaziv,  G.  Gans,  Y.  Reifen,  RAFAEL 

2199  The  Electromagnetic  Launch  Trends  Utilization  for  Shaped  Charge  Jets  Penetration  Depth  Decrease 
S.V.  Demidkov,  Effective  Soft  Ltd. 

3003  Simulation  of  the  Perforation  of  Low  Mass  Long  L/D  Rods  Against  Finite  RHA  Plates 

P.  Church,  I.  Cullis,  A  Bowden,  D  Gibson,  QinetiQ,  Ltd. 

3011  Deflection  and  Fracture  of  Tungston  Rods  by  Yawed  Impact 

S.  Bless,  R.  Russell,  Instutite  for  Advanced  Technology;  K.  Tarcza,  US  Army  ARDEC;  E.  Taleff,  Department  of 
Mechanical  Engineering,  The  University  of  Texas  at  Austin;  M.  Huerta,  The  University  of  Texas  at  El  Paso 

301 2  Anomalies  in  the  Strength  of  Alumina  under  Dynamic  Compression 

T.  Beno,  S.  Bless,  S.Nichols,  Instituite  for  Advanced  Technology 

301 4  On  the  3D  Visualization  of  Ballistic  Damage  in  TI-6AL-4V  Applique  Armour  with  X-Ray  Computed  Tomography 
J.M.  Wells,  JMW  Associates;  W.H.  Green,  N.L.  Rupert,  US  Army  Research  Laboratory,  Weapons  &  Materials  Research 
Division;  John  M.  Winter,  Jr.,  ORISE  Contractor  at  WMRD;  S.J.  Cimpoeru,  DSTO  Melbourne 

401 2  Analytical  Models  for  Foam,  Ice  and  Ablator  Impacts  into  Space  Shuttle  Thermal  Tiles 
James  D.  Walker,  Sidney  Chocron,  Walt  Gray,  Southwest  Research  Institute 

4013  CTH  Simulations  of  Foam  and  Ice  Impacts  into  the  Space  Shuttle  Thermal  Protection  System  Tiles 
Sidney  Chocron,  Walt  Gray,  James  D.  Walker,  Southwest  Research  Institute 

4016  Damage  Created  in  Composite  Sheet  by  Explosives  -  Effects  of  Fibre  Type,  Explosive  Mass  and  Attenuating 
Material 

M.  R.  Edwards,  R.  Unwin,  Centre  for  Materials  Science  and  Engineering,  Cranfield  University 


Warhead  Mechanisms  Poster  Session 
10:20  am  - 12:00  pm,  November  16 

1914  Experimental  Investigation  of  Equivalent  Blast  Characteristics  for  Aluminiferous  Explosive  in  Shallow  Underwater 
Wenbin  Gu,  Jianqing  Liu,  Qingli  Su,  Weiping  Zhou,  Ballistic  Research  Laboratory  of  China,  Nanjing  University  of 
Science  &  Technology 

1 935  Break-up  of  Copper  Shaped  -  Charge  Jets:  A  Combined  Experimental/Numerical/Analytical  Approach 
Jacques  Petit,  Centre  d’Etudes  de  Gramat;  V.  JeanClaude,  C.  Fressemgeas,  Laboratoire  de  Physique  et  Mecanique  des 
Materiaux,  Universite  de  Metz/CNRS 

1 949  A  Theoretical  Analysis  for  Initial  Fragment  Velocity  and  Peak  Overpressure  of  a  Blast  Fragmentation  Device 
Jin  Jianming,  Institute  of  Structural  Mechanics,  China  Academy  of  Engineering  Physics 

1 963  Scaling  the  Dynamic  Response  of  Armored  Vehicle’s  Floor  Subjected  to  a  Large  Buried  Charges 

Avidov  Neuberger,  IMOD,  MANTAK,  Tank  Program  Management;  S.  Peles,  IMI  Central  Laboratory  Division;  D.  Rittel, 

Technion,  Israel  Institute  of  Technology,  Faculty  of  Mechanical  Engineering 

1984  High-Speed  Flash  X-Ray  Computed  Tomo-Cinematography 

Philip  Heiberg,  Karsten  Michael,  Fraunhofer  Institute  for  High-Speed  Dynamics,  Ernst-Mach-lnstitut 

2002  The  Influence  of  Parameters  Other  Than  Liner  Velocity  on  Shaped  Charge  Jet  Coherence 

Frederik  Mostert,  CSIR  Defencetek;  C.  J.  Terblanche,  Denel  Land  Systems  -  Western  Cape;  M.  F.  Maritz,  Department  of 

Applied  Mathematics,  University  of  Stellenbosch 

2053  Comparison  of  Vulnerability  and  Performances  of  Insensitive  Munitions  (IM)  and  Non  IM  Directed  Energy 
Warheads 

Frederic  Peugeot,  MSIAC,  NATO  HQ 

2059  Trumpet  Shaped  Liners’  Influence  on  Slug  Properties 
Eitan  Hirsch,  Consultor;  Meir  Mayseless,  IMI  Central  Laboratory 

2077  A  Study  on  the  Structure  of  Small  Caliber  EFP 

Chen  Zhigang,  Zhao  Taiyong,  Hou  Xiucheng,  Dong  Surong,  The  Research  Institute  of  Explosive  Demolition  & 

Defence  Technology,  North  University  of  China;  You  Zheng,  Department  of  Precision  Instruments  and  Mechanics, 
Tsinghua  University 

2078  A  Framework  for  the  Analyses  and  Visualization  of  X-Ray  Computed  Tomography  Image  Data  using  a  Compute 
Cluster 

Jeffrey  R.  Wheeler,  US  Army  Research  Laboratory;  William  H.  Green,  US  Army  Research  Laboratory,  Weapons  & 
Materials  Research  Division;  Michael  Schuresko,  Baskin  School  of  Engineering,  University  of  California;  Michael  Patrick 
Lowery,  Computational  Mathematics  Department,  University  of  California 

2137  An  Artillery  Shell  for  Anti-Bunker  Applications  (155  ABS) 

Remi  Boulanger,  Giat  Industries;  Anders  Vangen  Jordet,  Dagfinn  Hoff,  NAMMO  Raufoss 

2141  Development  of  a  TMRP-6  Surrogate  Mine 

Yves  Baillargeon,  A.  Sirois  and  G.  McIntosh ;  DRDC  Valcartier 

2148  Use  of  Foamed  TNT  Mixtures  as  a  Dispersion  Charge  of  Submunitions 

Jun  Sik  Hwang,  S.-K.  Kwon,  C.-K.  Kim,  S.-W.  Kwon,  S.-S.  Kim,  S.-H.  Moon,  Explosive  Trains  and  Gun  Propellant 
Team,  Agency  for  Defense  Development 

2152  Wall  Breaching  Tandem  Warhead 

Andreas  Helte,  Torgny  Carlsson,  Hakan  Hansson,  Svante  Karlsson,  Jonas  Lundgren,  Lars  Westerling,  Hakan 
Ornhed,  Swedish  Defence  Research  Agency,  FOI,  Weapons  and  Protection  Division 

2156  An  Analytical  Penetration  Model  for  Jets  with  Varying  Mass  Density  Profiles 

Milton  F.  Maritz,  Stellenbosch  University;  Klaus  D.  Werneyer,  TTP  Products;  Frederik  J.  Mostert,  Defencetek 


2168  Peripheral  Initiation  Technology  Development 

Arthurs.  Daniels,  Ernest  L.  Baker,  William  J.  Poulos,  Vladimir  M.  Gold,  B.  Fuchs,  US  Army  ARDEC 

2175  Shaped-Charge  Jet  Stability  Calculations:  The  Role  of  Initial  and  Boundary  Conditions 
James  S.  Stolken,  S.  Christian  Simonson,  Mukul  Kumar,  Lawrence  Livermore  National  Laboratory 

3007  Enhanced  Focused  Fragmentation  Warhead  Study 

Richard  Fong,  William  Ng,  Peter  Rottinger,  Steve  Tang,  US  Army  Armaments  Research,  Development  and  Engineering 
Center 

401 5  Microstructure  and  Properties  of  the  Explosively  Formed  Petals  in  Aluminium  Alloys 
M.  R.  Edwards,  J.  M.  Cassar,  Centre  for  Materials  Science  and  Engineering,  Cranfield  University 

Interior  Ballistics/Launch  Dynamics  Poster  Session 
8:10  am  -  9:50  am  November  17 

1 01 7  In-Bore  Mechanics  Analysis  of  the  M855  Projectile 
Joseph  T.  South,  James  F.  Newill,  US  Army  Research  Laboratory 

1020  Dynamic  Strain  Measured  in  a  105-mm  Composite  Gun  Barrel  -  A  Fiction  or  Reality 
Jerome  T.  Tzeng,  US  Army  Research  Laboratory 

1916  A  Vector  Way  for  Calculating  Propellant’s  Combustion  Performance 

Wei  Zhifang,  Department  of  Mechanical  and  Electronic  Engineering,  North  University  of  China 

1931  Interior  Ballistics  Code  Applied  to  ETC  Concept:  Computations  and  Validations 
Gilies  Legeret,  Dominique  Boisson,  Giat  Industries 

2004  The  FHIBS  Internal  Ballistics  Code 

Clive  R.  Woodley,  Steve  Billett,  QinetiQ;  Caroline  Lowe,  Department  of  Applied  Mathematics  and  Theoretical  Physics, 
Centre  for  Mathematical  Studies,  University  of  Cambridge;  William  Speares,  The  Cylinders;  Eleuterio  Toro  Laboratory  of 
Applied  Mathematics,  Faculty  of  Engineering,  University  of  Trento 

2005  Modelling  the  Ignition  of  40mm  Gun  Charges 
Clive  R.  Woodley,  QinetiQ 

2020  Thermo-Mechanical  Erosion  Study  of  the  120mm  Chromium  Coated  Gun  Barrel:  Computation  and  Validation  of  the 
Heat  Exchange  Boundary  Condition 

Dominique  Boisson,  Gilies  Legeret,  Roxan  Cayzac,  Giat  Industries 

2031  MOBIDIC-NG:  A  1 D/2D  CFD  Code  Suitable  for  Interior  Ballistics  and  Vulnerability  Modelling 
Baptiste  Longuet,  Pascal  Millet,  Eric  Taiana,  ETBS;  Patrick  Della,  Pieta  Christiane  Reynaud,  SNPE  Materiaux 
Energetiques  CRB;  Patrice  Franco,  Alain  Carrere,  Institut  Franco-Allemand  de  Recherches  de  Saint-Louis  (ISL) ;  Gilies 
Legeret,  Dominique  Boisson,  Giat  Industries;  Alexandre  Papy,  ERM  ABAL  30 

2070  Barrell  Life  Results  of  the  5.56  mm  XC77A1  Cartridge 
Etienne  Munger,  SNC  Technologies,  Inc. 

2117  Further  Investigation  of  the  Effect  Known  as  Electrothermal  Pyrolysis 

Steve  R.  Fuller,  M.J.  Taylor,  QinetiQ 

2150  Determination  of  Force  and  Temperature  Impact  on  Missile’s  Fuel  Charge  in  Process  of  Ignition 
Dmitriy  Orlov,  GDT  Software  Group 

2161  Unsteady  Intermediate  Ballistics:  2D  and  3D  CFD  Modelling,  Application  to  Sabot  Separation 
Roxan  Cayzac,  Eric  Carette,  Giat  Industries,  Division  Munitions;  Thierry  Aiziary  de  Roquefort,  Universite  de  Poitiers, 
Laboratoire  d’Etudes  Aerodynamiques;  Philippe  Bidorini,  Emmanuel  Bret,  Pascal  Delusier,  Serge  Secco,  DGA/ETBS, 
Direction  de  I’Expertise  Technique 


2166  Large  Caliber  Firing  with  Electro  Thermal-Chemical  Ignition  (ETI) 

Jonathan  D.  Shin,  John  J.  O’Reilly,  David  T.  Keyser,  US  Army  Research,  Development  and  Engineering  Center  - 
TACOM;  Jahn  Dyvik,  United  Defense  L.R 

3006  Rail  Gun  Test  Projectile  for  Improved  Developmental  Testing  of  Precision  Munition  Electronics 

T.  Myers,  D.  Carlucci,  J.A.  Cordes,  US  Army  ARDEC,  Analysis  and  Evaluation  Division,  Fuze  and  Precision  Munitions 

Technology  Directorate 

4010  Improved  Mortar  Barrel  Thermal  Model 

M.  Pocock,  C.  Guyott,  Frazer-Nash  Consultancy  Ltd;  P.  Locking,  BAE  Systems,  Land  Systems 

Vulnerability/Lethality/Wound  Ballistics  Poster  Session 
10:20  am  - 12:00  pm  November  17 

1855  On  Incorporating  XCT  into  Predictive  Ballistic  Impact  Damage  Modeling 
Joseph  M.  Wells,  JMW  Associates 

1878  New  Soft-Target  Failure  Criteria  for  System-Analytical  Considerations 

Markus  J.  Estermann,  RUAG  Defence,  Warhead  Division;  Beat  P.  Kneubuhl,  Aramasuisse 

1941  Protecting  Vehicles  from  Landmine  Blasts 

Sheri  L.  Hlady,  Denis  Bergeron,  Defence  R&D  Canada  -  Suffield;  Rene  Gonzalez,  US  Army,  PM  Light  Tactical  Vehicles 

1957  Office  of  Naval  Research  Limb  Protection  Program 

Graham  K.  Hubler,  NRL 

1980  Survivability  and  Lethality  Assessment  Software  Based  on  Virtual  Mode  Technology 
Lu  Yonggang,  Qian  Lixin,  Yang  Yubin,  Liu  Tong,  Institute  of  Structural  Mechanics,  CAEP 

1981  “TBM-Xpert”  -  A  New  Endgame  Code:  Features  and  Validation 

Werner  Arnold,  EADS-TDW  Gesellschaft  fur  verteidigungstechnische;  E.  Rottenkoiber,  NUMERICS  GmbH 

1989  The  Use  of  Ballistic  Knowledge  in  Ammunition  Safety  Cases 
Martin  van  de  Voorde,  TNO  Defence,  Safety  and  Security 

2011  A  Note  on  the  Roecker-Ricchiazzi  Model  of  Penetrator  Trajectory  Instability 
William  J.  Flis,  DE  Technologies,  Inc. 

2022  Numerical  Calculation  and  Simulation  of  Missile  Jet-Airplane  Interaction 
Feipeng  Zeng,  Faculty  of  Artillery  Command,  Nanjing  Artillery  Academy 

2111  Need  for  Enhanced  Protection  Against  Blast  Threats  for  Soldiers  Exposed  to  Roadside  Improvised  Explosive 
Devices  (lEDs) 

Frangois-Xavier  Jette,  Jean-Philippe  Dionne,  Aris  Makris,  Med-Eng  Systems,  Inc.;  Karl  Masters,  PEO  Soldier; 
Christine  Perritt,  PM  Soldier  Equipment 

2119  WitnessMan:  The  Software  Tool  to  Design,  Analyse  and  Assess  a  Witness  Pack  with  Respect  to  Military  and 
Medical  Effects  on  an  (Un)protected  (Dis)mounted  Soldier. 

Theo  L.A.  Verhagen,  R.  Kemper,  H.  Huisjes,  S.G.  Knijnenburg,  A.  Pronk,  M.H.  van  Klink,  TNO  Defence,  Security  and 
Safety 

2154  Injury  Risks  Resulting  from  Deminer  Position 

Frangois-Xavier  Jette,  Jean-Philippe  Dionne,  Ismail  El  Maach,  Aris  Makris,  Matt  Ceh,  Med-Eng  Systems,  Inc.;  Denis 
Bergeron,  Defence  R&D  Canada  Suffield 

21 63  RPG  Mitigation  for  Military  Vehicles 

Karl  Pfister,  Dipl.  Ing  (FH)  Armatec  Survivability  Corporation 


21 64  Protection  Against  Closely-Spaced  Impacts  by  Small  Arms  Bullets 

Michael  J.  Iremonger,  Cranfield  University,  Royal  Military  College  of  Science;  Abdullah  Alsalmi 

301 3  Vulnerability  Evaluations  of  30mm  Airburst  Ammunition 

Quoc  Bao  Diep,  Eimund  Smedstad,  Nammo  Raufoss  AS;  Nick  Rogers,  System  Design  Evaluation  (SDE) 

3018  Challenges  and  a  Solution  in  Determining  Land  Mine  or  I  ED  Neutralization  Effectiveness 
Robert  Colbert,  Mark  Majerus,  William  Clark,  DE  Technologies,  Inc. 

4011  Numerical  and  Experimental  Analysis  of  the  Detonation  of  Sand-Buried  Mines 
N.  Heider,  A.  Klomfass,  Fraunhofer-lnstitut  fur  Kurzzeitdynamik,  Ernst-Mach-lnstitut 


The  symposium  registration  fees  are: 


Regular  Late/Onsite 
after  10/28/05 

$950  (US)  $1045  (US) 


IBC  Committee  Appreciation  Dinner 
Guest  at  Both  Receptions 
Guest  at  One  Reception 


$75  (US) 
$75  (US) 
$50  (US) 


The  symposium  registration  fee  includes  attendance  at  all  sessions,  bound  symposium  proceedings  with  CD, 
continental  breakfasts,  coffee  breaks,  lunches,  receptions,  and  administrative  costs.  The  registration  fee  will 
also  include  a  compact  disc  (CD)  which  contains  a  cumulative  database  of  titles,  authors  and  abstracts  of  all  of 
the  22  Ballistics  Symposia. 

To  register  online  for  this  conference  visit:  http://reaister.ndia.ora/interview/reaister.ndia7~Brochure~6210.  You 
can  also  visit  the  NDIA  website  at  www.ndia.org  and  select  “Schedule  of  Events”.  Then  select  2005  November 
and  scroll  down  to  the  22nd  International  Symposium  on  Ballistics.  Once  there,  select  the  blue  “Register”  link  in 
the  lower  left  hand  corner  of  your  screen.  Review  your  information  and  then  select  “submit”  one  time  only  and 
then  select  “confirm”.  On-line  registration  will  close  at  5:00  pm  EST  on  October  28,  2005.  You  must  register 
on-site  after  this  date. 

-or- 

You  may  fax  the  completed  registration  form  contained  in  this  brochure  to  (703)  522-1885. 

-or- 

You  may  mail  the  completed  registration  form  contained  in  this  brochure  to:  Event  #  6210,  National  Defense 
Industrial  Association,  2111  Wilson  Boulevard,  Suite  400,  Arlington,  VA  22201 -3061. 

Payment  must  be  made  at  the  time  of  registration.  Registrations  will  not  be  taken  over  the  phone. 

Cancellations  and  Refund  Policy 

Registrants  who  cannot  attend  the  22nd  International  Symposium  on  Ballistics  must  provide  written  notification 
via  email  to  bbommelje@ndia.org  or  fax  to  (703)  522-1885  on  or  before  September  16,  2005  to  avoid  a 
cancellation  fee. 

Cancellations  received  between  September  16,  2005  and  October  28,  2005  will  receive  a  refund  minus  a  $75 
cancellation  fee.  No  refunds  will  be  given  to  cancellations  recieved  after  October  28,  2005  however, 

SUBSTITUTIONS  ARE  WELCOME  IN  LIEU  OF  CANCELLATIONS. 

You  must  have  a  government  picture  identification  (drivers  license,  passport,  military  ID,  etc.)  to  receive  a 
symposium  badge.  Badges  must  be  worn  at  all  times  during  the  symposium. 

Special  Needs 

NDIA  supports  the  Americans  with  Disabilities  Act  of  1 990.  Attendees  with  special  needs  should  call 
(703)  522-1820  prior  to  October  3,  2005. 

Hotel  Accommodations 

A  limited  block  of  rooms  have  been  reserved  at  the  Fairmont  Waterfront  Hotel.  The  industry  room  rate  is  $219 
Canadian  (approximately  $180  US).  The  government  symposium  room  rate  is  approximately  $114  Canadian 
($94  US).  Please  call  (604)  691-1991  to  make  reservations. 

In  order  to  ensure  the  discounted  NDIA  rate,  please  make  reservations  early  and  ask  for  the  NDIA  room  block. 
Rooms  will  not  be  held  after  Tuesday,  October  11 , 2005  and  may  sell  out  before  than.  Rates  are  also  subject 
to  increase  after  this  date. 


*The  government  room  rate  applies  only  to  active  duty  military  and  civilian  government  employees.  It  is  not 
available  to  government  contractors,  retired  military  or  retired  civilian  government  employees.  ID  cards  and/or 
travel  orders  will  be  required  at  check-in  to  verify  rate  eligibility. 


GENERAL  INFORMATION 


Symposium  Attire 

Appropriate  dress  for  this  symposium  is  business  for  civilians  (coat  and  tie)  and  class  A  uniform  for  military. 

Inquiries 

For  more  information  regarding  the  symposium  contact  Britt  Bommelje,  Meeting  Planner  at  (703)  247-2587  or 
bbommelje@ndia.org. 

Promotional  Partnerships 

Increase  your  company  or  organization  exposure  at  this  premier  event  by  becoming  a  Promotional  Partner.  A 
Promotional  Partnership  ($5,000)  will  add  your  company  name  to  the  back  cover  of  the  on-site  brochure  as  well 
as  main  platform  recognition  throughout  the  conference,  signage  at  all  events  including  the  opening  reception 
and  a  350-word  organization  description  in  the  conference  agenda.  For  more  information,  please  contact  Sam 
Campagna  at  703-247-2544  or  scampagna@ndia.org. 

www.defensejobs.com 

The  Defense  Industry’s  leading  employment  website;  find  a  job,  post  a  job  listing,  post  a  resume,  and  search 
resumes.  For  more  information  please  contact  info@defensejobs.com  or  (703)  247-9461 .  Please  visit 
www.defensejobs.com 

IBC  Committee  Appreciation  Dinner 

The  IBC  Committee  Appreciation  Dinner  will  be  held  on  Friday,  November  18,  2005  in  Vancouver.  This  dinner 
is  open  to  IBC  Committee  members  and  their  guests  only.  If  you  and  your  guest  would  like  to  attend,  please 
make  a  note  of  it  on  the  registration  form.  There  is  a  $75  charge  per  person  to  attend. 


“The  Department  of  Defense  finds  this  event  meets  the  minimum  regulatory  standards  for  attendance 
by  DoD  employees.  This  finding  does  not  constitute  a  blanket  approval  or  endorsement  for  attendance. 
Individual  DoD  component  commands  or  organizations  are  responsible  for  approving  attendance  of  its 
DoD  employees  based  on  mission  requirements  and  DoD  regulations.” 
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22nd  International  Symposium  on  Ballistics 

Vancouver  Convention  Center  ~  Vancouver,  BC,  CANADA 

November  14-18,  2005 

Event  #  621 0 

Organization  Information 


Organization  Name  (as  it  should  appear  on  booth  sign  --  limited  to  40  characters  and  spaces) 

Point  of  Contact  (for  fees  and  exhibitor  service  kit) 

Title 

Street  Address 

City 

State 

Zip  Code 

Telephone 

Fax 

E-mail 

Reserve  your  booth  and  register  your  exhibit  staff  at  http://exhibits.ndia.org 


Exhibit  Space  Information: 

Corporate  Members  and  bona-fide  government  organizations: 
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Plasma  Ignition  of  a  30-mm 

Cannon 

Richard  A.  Beyer 
Andrew  L.  Brant 
Joseph  J.  Colburn 


US  Army  Research  Laboratory 
Aberdeen  Proving  Ground,  Maryland 

22nd  International  Symposium  on  Ballistics 
14-18  November '  2005 
Vancouver BC 


Goal  of  This  Study 


•  Apply  the  knowledge  learned  in  our  earlier 
Plasma-Propellant  Interaction  research 
to  the  Ignition  of  a  30-mm  cannon 

•  Start  to  Develop  "Design  Rules"  for 
plasma  ignition 


Ultimately:  Make  plasma  ignition  in  large- 
caliber  guns  more  efficient 


Hardware 


•  Power  Supply 

-  Standard  PFN  600  mF,  3  kv,  3  kj 

-  Short  (300  ms)  or  Long  (900  ms)  Pulse 

•  Propellants 

-  JA2,  M30,  Graphite-Free  JA2 

•  Standard  Ignition  for  30-mm  cannon 

-  M52  primer  w/1.5  grams  Benite  in  bayonet 


Small-Scale  Experiments 


Plasma  into  bayonet  igniter  tube 


Schematic  of  Test  “Igniter” 


pressure  gauge 


vent  holes 


anocle  ca™ode  (4  sides) 


Designed  to  look  like  an  igniter  we  can  put  into  30-mm  cannon. 

32  total  1/16-inch  vent  holes  (8  per  side) 
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Response  with  1.5  grams 
Standard  JA2  in  tube 

Parameter  =  discharge  energy 
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Response  with  1.5  grams 
Graphite-Free  JA2  in  tube 

Parameter  =  discharge  energy 
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(1.5  g  JA2  »  8.5  kj  chemical  energy) 


30 -mm  Cannon 


30 -mm  Stub -Case  Cartridge 
for  Plasma  Injection 


cathode 


capillary 


stub  case 


polycarbonate  cartridge  case 


propellant  in  primer  tube  propellant  bed 


Plasma  Discharge  Energy:  2.1  to  2.3  kj 

Propellant  in  primer  tube:  1.6  to  1.9  grams  (8.1  to  9.6  kj  chem  energy) 

Chamber  volume:  135  cc 

Projectile  mass:  550  g 

Shot -Start  Pressure:  very  low 

Conventional  ignition  baseline:  M 52  primer  and  1.5  grams  Benite  in  bayonet 


M30  Propellant  Charge 

( GF-JA2  in  plasma  igniter  tube) 
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time  (jits) 


JA2  Propellant  Charge 

GF-JA2  vs.  Standard  JA2  in  igniter  tube 


pressure  (MPa) 


JA2  Propellant  Charge 
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JA2  Propellant  Charge 

Soda  straw  w/  SF-JA2  "box" 
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JAZ  Propellant  Charge 

Long  plasma  pulse 
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Early  Time  Behavior 

with  short-pulse  plasma  ignition 
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KEY:  red  line:  breech  blue  line:  projectile  gray  line:  discharge 
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Early  Time  Behavior 

with  long-pulse  plasma  ignition 


Summary 


•Graphite  Free  JA2  provides  higher/faster  pressurization  inside 
igniter  tube 

•But ...  The  differences  not  as  dramatic  in  gun  firings 

•Plasma  ignition  more  prompt  than  conventional  ignition 

•At  low  (2  kj)  plasma  energy 
•With  M30  “boosted”  with  JA2  igniter 

•Relative  quickness  of  ignition  follows  our  intuition  regarding 
strength  of  plasma-propellant  interactions 

•GF  faster  than  STD  JA2 

•Filled  igniter  tube  faster  than  empty 

•Short  pulse  faster  than  long  pulse 

•Is  faster  better? 

•Short  pulse:  Fast  and  more  PPI 

•Long  pulse:  Only  slightly  slower;  smoother  in  our  cannon 


•Key  to  Plasma  Ignition:  Rapid  Pressurization  of  Charge 
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Recent  Computations  and  Validations 

of 

Projectile  Unsteady  Aerodynamics 

R.Cayzac,  E.Carette,  R.Thepot  and  P.Champigny 


G  I A  T  industries 


ONERA 


CFD  Computations  of  Projectile  Unsteady  Aerodynamics 


♦  OBJECTIVES 

♦  THEORETICAL  AND  NUMERICAL  APPROACHES 

♦  WIND  TUNNEL  TESTS 

♦  VALIDATION  RESULTS 

>  YAWING  AND  SPINNING  PROJECTILES 

>  KINETIC  PROJECTILES 

♦  CONCLUSIONS 
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Main  Aerodynamic  Coefficients  Concerned 


•  FORCE  COEFFICIENTS 

>CA  =  CA(a=0)  +  ACA  (M,  a,  P) 

>CN  =  CNa(M)a  +  CNq(M).qD/V  +  ACN  (M,  a,  4>) 

>CY  =  Cypa(M).p.a.D/V  +  ACY  (M,  a,  (|>) 

•  MOMENT  COEFFICIENTS 

>Cm  =  Cma(M).a  +  ACm  (M,  a,  §)  +  Cmq(M).qD/V 
>Cn  =  Cnpa(M).p.a.D/V  +  ACn  (M,  a,  <|>)  +  Cnr(M).rD/V 
>CI  =  Clo(M)  +  ACI  (M,  a,  4>)  +  Clp(M).pD/V 


•  DYNAMIC  COEFFICIENTS^  DAMPINGS,  MAGNUS  AND 
PSEUDO-MAGNUS  EFFECT 
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ONERA 


•  Theoretical  and  Numerical  Approaches  (FLU3M  and  elsA) 

>  3D  Navier-Stokes  equations,  RANS  and  URANS 

>  Baldwin-Lomax,  Spalart-Allmaras,  k-oo,..,  turbulence  models 

>  Finite  volume  method 

>  Cell-centered  discretization  in  an  absolute  frame 

>  Fully  implicit  Gear  scheme  (1st  and  2nd  order) 

>  (Pulliam  under-iteration  technique  (URANS),  grid  movement) 

•  Grid  :  multiblocks  structured  with  hexahedral  cells 

>  Wall  cell  »  1  )nm 

>  30  to  50  cells  in  the  boundary  layer 

>  Stretching  factor  <1.2 

>  Y+ »  1  (a  posteriori  criterion) 

>  Up  to  5,000,000  cells 

•  Computational  Performances 

>  »  1  jLis/cel I/iteration  NEC  SX-6 

>  »  45  (is/cel I/iteration  SGI  Octane  RISC  12000 

>  »  20  jLis/cel  I/iteration  Cluster  of  Xeon  (2.2  GHz) 
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ONERA 


Main  Aerodynamic  Coefficients  Concerned  and  CFD  Generalities 


I 

•  Cell-centered  discretization  of  RANS  equations  expressed  in  an 
absolute  framework  R,  grid  partition  in  blocks  of  rigid  hexahedral  cell 


Validation:  Wind  Tunnel  Tests 


ONERA:  S3MA 

•  Static  coefficients 

>CA,  CN,  Cy,  Cl,  Cm,  Cn 

•  Dynamic  coefficients 

>  Cy  (Cypa),  Cn  (Cnpa) 
>CNq+CNa,  Cmq+Cma 
>Clp 


•Mach  number 
•Reynolds  number 
•Roll  rate  p  (P*=p.D/vj 

•  Pitch  rate  q 

•  Damping  frequency 
•Angle  of  attack 
•Roll  position  § 
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Example  of  Yawing  and  Spinning  Projectile  Grid 


ON  E  R  A 
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Example  of  Kinetic  Projectile  Grid 


90  Blocks 
£,700,000  cells 


Y+  « 1 


40  cells  in  the  boundary  layer 
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ON  E  R  A 


Example  of  Y+ A  Posteriori  Verification  Criterion 


mm  m 

: 

*  y+  distribution 

(Mach  =  0.90,  a  =  3°,  Pi  =  1.20E+05  Pa,  Ti  =  300.0  K) 
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Yawing  and  Spinning  Projectiles:  Magnus  Validation  Results 


•  ^'Navier-Stokes  Computations  and  Validations  of  a  Yawed  Spinning  Projectile", 

18th  International  Symposium  on  Ballistics,  15-19  November,  San  Antonio,  Texas,  USA,  1999. 

•  2"Recent  Developments  on  Aeroballistics  of  Yawing  and  Spinning  Projectiles:  Part  I,  Wind 
Tunnel  Tests  ",  20th  International  Symposium  on  Ballistics,  7-11  October,  Orlando,  USA,  2002. 

•  3"Recent  Developments  on  Aeroballistics  of  Yawing  and  Spinning  Projectiles:  Part  II,  Free 
Flight  Tests  ",  20th  International  Symposium  on  Ballistics,  7-11  October,  Orlando,  USA,  2002. 

•  4"Recent  Developments  on  Aeroballistics  of  Yawing  and  Spinning  Projectiles:  Part  III, 

Validation  Results",  20th  International  Symposium  on  Ballistics,  7-11  October,  Orlando,  USA, 
2002. 

•  5"Analysis  of  Static  and  Dynamic  Stability  of  Spinning  Projectiles",  21st  International 
Symposium  on  Ballistics,  Adelaide,  Autralia,  19-23  April ,  2004. 


•  Agreement  between  computations  and  experiments  is  satisfactory  at  moderate 
angles  of  attack  on  the  Magnus  dynamic  coefficients 

•  Strong  difficulties  at  high  incidence  (up  to  10°)  and  in  the  transonic  regime 
•RANS  ->  MILES,  DES,  etc.,  simulations  (ARL,  De  Spirito,  etc.) 
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Example  of  Validation  Results:  Cmq,  Pseudo-URANS 


Supersonic  Conditions:  Wind  Tunnel  C4 /  FLU3M  Baldwin-Lomax 
(M^  2.89,  a  =  0°  (1  °for  static),  Pi  =  3.224  bar,  Ti  =  299  K,  0  ±  1  °,  f  =  9  to  13  Hz) 


Transonic  Conditions:  Wind  Tunnel  S3Ma/  elsA  k-co 
(M^  0.9,  a  =  3°,  Pi  =  1.2  bar,  Ti  =  300  K,  0  ±  1  °,  f  =  6  a  10  Hz) 
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Example  of  Validation  Results:  Magnus  of  Kinetic  Projectiles 


e  I V 1 


Q 


•  Magnus  effect  origin  on  fins 

•  Interaction  with 
asymmetric  fuselage  wake 

•  Modifications  of  the  local 
incidences  induced  by 
spin  (apex  shock  &  tip 
vortex) 

•  Flow  description 

•  I  Apex  shocks 

•  II  Tip  vortex 

•  aiocai=  f(a,5,r,0)  alocal  =  a  ±  8  ±  r.  Q/Vro 

•  aN  =  tan1  (tan  alocal/cos  A) 

•  MachN  =  Mach^  V(  1-cos2  a,oca| .  sin2  A) 


•  URANS,  Baldwin-Lomax,  Mach  4.3,  p*=  0.041,  Pi  =  7.7  Bar, 
Ti  =  295  K,  a  =  4.22°,  L/D  =  12.5. 
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ONERA 


Force  Coefficient  C 


Example  of  Validation  Results:  Magnus  of  Kinetic  Projectile 


0  90  180  270  360 

Fin  Position  0  ° 


CN 


Normal  Flow 


Fin  azimuthal  position  influence  (variations  of  the  local  incidence)  on  force  coefficients 
URANS,  Baldwin-Lomax,  Mach  4.3,  p*=  0.041,  Pi  =  7.7  Bar,  Ti  =  295  K,  a  =  4.22°,  L/D  =  12.5. 
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Example  of  Validation  Results:  Magnus  of  Kinetic  Projectile 


POLARS  IN 
ANGLES  OF  ATTACK 


experiment 


POLARS  IN 
SPIN  RATES 


•  URANS  FLU3M  Baldwin-Lomax 

>  Mach  4.3,  p*  =  0  ^  0.041.  Pi  =  7.7  Bar,  Ti  =  295  K,  a  =  0°->  4,22°,  L/D  =  12.5. 
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Example  of  Validation  Results:  Normal  force  and  Magnus  Coefficients  of  Kinetic  Projectile 


Global  normal  force  is  independent  of  the  rotation 

Mach  4.3,  p*=  0.041  (100  rounds/s),  Pi  =  7.7  Bar,  Ti  =  295  K,  a  =  4.22°,  L/D  =  12.5. 


Different  contributions : 


—  —  —  Fuselage  (ogive/cylinder) 

—  - - Fuselage  (tail) 

pin  surface 

—  —  —  Leading  edge  &  lip  surface 

- - -  All  surfaces 


Numerical  convergence 

•  7200  iterations/round 

•  6  under-iterations/iteration 
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Different  contributions  : 


—  —  —  Fuselage  (ogtve/cylinder) 

- - —  Fuselage  (tail) 

.  Fin  surface 

- Leading  edge  &  tip  surface 

-  All  surfaces 
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Example  of  Validation  Results:  Cmq,  Clp  and  Magnus  of  Kinetic  Projectiles 
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Example  of  Validation  Results:  Cmq,  Clp  and  Magnus  of  Kinetic  Projectiles 


CFD  (URANS  FLU3M+Balwin-Lomax  and  Elsa+k-co)  /  Wind  Tunnel  S3Ma  (L/D  =  30) 

Cmq:  4.5,  a  =  0°^  5°,  Pi  =  6  bar,  Ti  =  299  K,  9  ±  1.5°,  f  =  2.2  4.2  Hz  (FLU3M) 

Clp:  M „=  4.5,  a  =  0°^  5°,  Pi  =  6  bar,  Ti  =  299  K,  p  =  1.5  ^  55  Rd/s  (elsA) 

Magnus:  4.5,  a  =  -1  5.5°,  3°,  Pi  =  6  bar,  Ti  =  356  K,  p  =  10->  90  Rd/s  (65)  (FLU3M) 


Mach  number  field 


TESTS 

CFD 

Errors 

(%) 

Clp 

-15 

-15.35 

2.3 

Cmq 

-2800 

-2619 

6.5 

Cy 

0.045 

0.0474 

-5.3 

Cn 

-1.071 

-1.112 

-4 

•  Periodic  behaviour  of  the  Magnus  and  roll  moment 
coefficients  during  unsteady  computation 

•  Mach  4.4,  p*=  0.02,  Pi  =  6  Bar,  Ti  =  360  K,  a  =  2°,  L/D  =  30. 
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ONE  R  A 


1000 


£000 

iteration 


3000 


-0.45 

-0.475 

-0.5 

-0.525 

-0.55 

-0.575 


4000 


o  10 


•  Extensive  wind  tunnel  database  for  validation,  need  free 
flight  results 

•  Progress  in  CFD  allows  us  to  predict  projectile  unsteady 
aerodynamics.  With  respect  to  spinning  and  kinetic 
projectiles,  a  demonstration  of  the  capability  of  the 
numerical  approach  was  carried  out 

•  Satisfactory  agreement  between  computations  and 
experiments  on  the  pitch  and  roll  dampings  and  on  the 
Magnus  effect 
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Example  of  Validation  Results:  Clp  of  Spinning  Projectile 


#  Clp  -  ( C I (p))/ ( P - L ref/ V 

CFD  (RANS,  Baldwin-Lomax),  p*  «  0.2,  a  =  0°) 
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ON  E  R  A 


CFD  Prediction  of  the  Pitch  Damping  Coefficient:  Pseudo-URANS 


•  Cmq  =  (Cm(q)  -  Cm(q=0))/(q.Lr6f/VJ 

(Mach  =  0.90,  a  =  3°,  Pi  =  1.20  Bar,  Ti  =  300  K) 


>  q  =  10  rad/s 
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DELTA  WING:  LEEWARD  SIDE  FLOW  TOPOLOGY 


■tr A3L.J 


CLASSIFICATION  OF  LEESIDE  FLOWFIELDS  FOR 
THICK  DELTA  WINGS  WITH  SHARP  LEADING  EDGES 


Q  Center  line 
shock  region 

0  Leading  edge 
separation 

0  Leading  edge 
separation 

0  Separation 
with  shock 

0  Shock  induced 
separation 

0  Leading  edge 
separation 
with  shock 


A=  70°,  Mach  4.3,  p*=  0.041,  Pi  =  7.7  Bar,  Ti  =  295  K,  a  =  4.22°^  MN  =  1.52  et  aN  =  0.27 
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Boundary  Conditions 


jtesrjaiie.-esfc 


•  Fuselage  state 

>no  slip  condition 


>adiabatic  wall 


> normal  pressure  gradient 


•  Fin  state 

>  slip  or  no  slip  conditions 


0  ^ 


•  Outer  boundary  condition  is  obtained  from  the  Theory  of 
Characteristics 
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^Aerodynamic  Characteristics  of  a  Grid  Finned 
Projectile  from  Free-Flight  Tests 
at  Supersonic  Velocities 


A.  Dupuis 
DRDC- V  alcartier 


Canada 


C.  Berner 
ISL 
France 


OUTLINE 


•  Background 

•  Model  Configuration 

•  Experimental  Tests 

•  Results 

•  Grid  finned 

•  Planar  (Air  Force  Finner) 


Conclusions 


BACKGROUND 


Grid  fins  are  being  studied 
by  several  organizations 

Favorable  lift  characteristics  at  high  a 
Low  hinge  moments 
Good  storability 

High  drag  (can  be  minimized  through  web  tailoring) 

Prior  research  include  wind  tunnel  tests  and  recently  CFD 
•  Results  presented  at  last  3  ISBs 


OBJECTIVES 


•  Obtain  Aero.  Coeffs.  &  Stability  Derivatives 

•  From  Free-flight  Tests 

•  DRDC  Aeroballistic  Range 

•  Establish  Data  base  to  compare  WT  andCFD  results 
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CONFIGURATION 


R  6.53 


\ 


R  0.007 


1 

*1  set  of  fins  canted  at  2.0°  to  produce  roll 

1 

0.7 

.1  set  of  fins  canted  at  0.5°  to  produce 

trim  forces  &  moments 
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PHYSICAL  PROPERTIES 


(A/B  model) 


d 

(mm) 

m 

(g) 

lx 

(g-cm2) 

Jy 

(g-cm2) 

1 

(mm) 

CG  from  nose 

(XCG  /l) 

30.0 

881.2 

954.55 

64517.07 

300.0 

0.405 

•  Length:  250  m 

•  Model  size:  5.56  mm  to  155  m 


Speeds:  up  to  Mach  7.0 

Temperature  et  humidity 
controlled 


DRDC  AEROBALLISTIC  RANGE 


►Instrumented  length:  220  m 


Section:  6  m  x  6  m 


54  Stations:  Indirect  orthogonal 

shadowgraphs 


4  Schlieren  stations 


at  least  4  firings/day 


RjL)  drdc  aeroballistic  range^3^ 

\' 

ARFDAS  -  Aeroballistic  Range  Facility 
Data  Analysis 


Physical  Properties 

L,M,D,I  ,1  ,/ 

*  y  z 

^xy’^xz’^yz 


Single  & 
Multiple  Fits 


6DOF  Dynamic  Data 

t,x,y,z,<f>,0,  Y 


ARFDAS 

Startup 


V 
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Linear  Theory 
Analysis 


/  6DOF  Symmetric  \ 
or 

\  6DOF  Asymmetric  / 


Cnja,Cma3,Cmaf5,C  mq  ,  C  m q  2  >  C  np  a  )  C  np  a  3  5  C  tip  a  5  ■ 

Cn/3t>Cnp'i)C  nr  iCnr2iCncciCmya  3  ?  C  n  ycc  3 . 

C\P  ’  C\P  a  2  ’  CXS  S’  ClJa  3 . 


Atmospheric 

Conditions 

P,P,T,/li 


Fit  Theoretical 
to  Experimental 


Aerodynamic  Forces  &  Moments  vs. 

Mach  No.  &  Angle  of  Attack  &  Roll  Angle 

c  Xo  ,Cxa2,Cxa4,CNa,CNaT,,CNa5,CYp  a  ,  C  Yp  a  3 


A/B  TESTS 


11  projectiles  fired 
-  Mach  1.4  to  3.5 


110  mm  Smooth  Bore  Gun 


First  Max  Yaws  ranged  from  2.4°  to  10.5 


Down  range  spin  rate  20.0  °/m 


A/B  TESTS 

VMUz  =  1215  m/s 


r \$y  A/B  TESTS 


M  =2.12 


RESULTS 


•  Grid  Finned  Results 

•  Compare  with  Air  Force  Finner 

•  Tested  in  DRDC  and  Eglin  Air  Force  Base 
Free-Flight  Ranges 
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CONCLUSIONS 


:V 


•11  Projectiles  Fired  Mach  1.4  to  3.5 

•  All  Main  Aerodynamics  well  Determined 

•  First  Time  for  Reliable  CMq  &  Clp 

•  Some  nonlinear  ones  also 


Comparison  with  Planer  Fin  Model 
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Multiple  Explosively  Formed  Penetrator 
(MEFP)  Warhead  Technologies  for  Mine  and 
Improvised  Explosive  Device  (I ED) 

Neutralization 
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Steve  Tang 
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OUTLINE 


AMC 

•  Overview 

•  Current  Neutralization  Techniques 

•  Program  Objective 

•  Technical  Approach 

•  MEFP  Warhead  Technology 

-Mine  Neutralization 
-IED  Neutralization 

•  Summary 


PAN  -  Percussion 
Actuated  Neutralizer 


M221  Clipped  Demo  Charge 


RE-70/MK40 


M82A1  Sniper  Rifle 
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Overview 

Current  lED/Bomb/Mine  Neutralization 
Techniques  &  Devices 


Ml  12  Demo  Blocks 


EOD  Robots 


No  Universal  Mine/IED 
Neutralizer  Exists  to 
Effectively  Defeat 
Threats  in  All 
Environments 


■  •  ■ 
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PROGRAM  OBJECTIVE 


Concealed  In  Dirt/Sand 


Concealed  In  Cloth/Sheathing 


Miniaturized 


To  understand  the  critical  MEFP 
parameters  required  to  cause  various 
reactions  (detonation,  deflagration, 
damage)  against  mines  and  Improvised 
Explosive  Devices  (lEDs). 


Concealed  In  Concrete 


Concealed  In  Plaster 


External  Components 
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Typical  Threats 


82mm  Mortar 


120mm  Mortar 


152mm  Projectile 
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TECHNICAL  APPROACH 


_ WHY  MEFP? 

•  Standoff  Capability 

•  Increase  Hit  Probability 

•  Design  Flexibility 

•  Lightweight  and  Compact 

•  Versatile  Deployment 


MEFP  Warhead  Pattern 
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DEMONSTRATED  CONTROLLED  FRAGMENTATION  PATTERNS 


Happy  Face  Simulation 
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Mine  Neutralization 


Design  Approach 

-Preliminary  Warhead  Evaluation 

-MEFP  Optimization  Study 

•  Mass 

•  Velocity 

•  Shape 

-Prototype  Design  &  Test 


Preliminary  Warhead 

7  MEFP  Configuration 
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MEFP  WARHEAD 

-  Effective  Range:  1-150  Meters 

-  Spheroids  to  Long  Rods  (L/D  1-3) 

-  On-Target  Pattern  Control 

-  Multiple  EFP’s  from  a  Single  Warhead 

-  Proven  Against  Many  Different  Light  Targets 


Case  / 
Explosive 


Liner 


Booster  /  Detonator 
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MEFP  Optimization 


Warhead  with  different  L/D 
MEFP  liner  design 


Stripper  Plate 


Single 

Penetrator 

X-rav 


Stripper  Plate 


Sand  Box  (Overburden) 
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MEFP  Optimization  Test  Results 

“5  Warhead  Designs  Tested” 


Flash 

Radiograph  of 
EFP 

Recovered 

MEFPs 
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No  Residual  Pen. 

□ 
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No  Residual  Pen. 
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3.85 

Perforated  Sandbox 
With  1/4”  Residual 
Penetration  in 
Armor 

AMC 


Prototype  Design 

“MEFP  Warhead” 


•  61  Multiple  EFPs  From  A  Single  Warhead 

•  On-Target  Pattern  Control  ~  Uniformly  Spaced 
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■p*  X 


\  I  i  /. 

\  %  i  /  # 

*  n+  >**  *- 


*  {  O'*' 
iS;s 


Simulated  Target  Pattern 


Prototype  Warhead  Test  Results 

“Residual  Penetration  Through  Sand  Overburden 
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Prototype  Warhead 

“Test  Setup” 


BLAST  SHIELD 


Standoff 


ANTI-TANK  MINE 


GRAVEL  OR  SAND 


SAND  OR  GRAVEL 


ANTI-TANK  MINE 


RESIDUAL  TARGET 
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System 

Demonstrator 

Requirement 


Prototype  Warhead  Test  Results 


Depth 

Sand  Overburden 

Gravel  Overburden 

A 

* 

* 

B 

* 

Detonation 

C 

Detonation 

* 

D 

* 

* 

E 

* 

* 

F 

* 

* 

G 

Detonation 

Detonation 

H 

Deflagrate 

* 

I 

* 

Deflagrate 

J 

No  Reaction 

* 

K 

Deflagrate 

Deflagrate 

L 

No  Reaction 

Deflagrate 

M 

No  Reaction 

No  Reaction 

N 

* 

No  Reaction 

*Not  tested 
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IED  Neutralization 


•  MEFP  Parametric  Study 

-  L/D  Study 

-Four  Prototype  Warhead  Designs 
Evaluated 
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Parametric  Study  Test  Matrix 


Overburden 

(Sand) 

Threat 

81mm 

120mm 

155mm 

*-  “  r'  -  -  : 

Cl 

$  ■  ^ 

o 

*5  O 

-  1 

L/D 

KE 

(kj) 

A 

D 

0" 

A 

D 

O" 

A 

D 

O" 

A 

D 

1 

29.6 

X 

X 

34.0 

X 

X 

X 

X 

X 

X 

X 

34.9 

X 

X 

X 

X 

48.5 

X 

X 

X 

X 

X 

1.5 

34.0 

X 

X 

X 

X 

X 

3 

26.0 

X 

X 

X 

X 

34.0 

X 

X 

X 

X 

X 

3+ 

34.1 

X 

X 

X 

X 

X 

X 

**  Shaded  Boxes  correspond  to  the  L/D  Study 
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120  mm  Surrogate  Simulations 

“Penetration  Analysis” 

Parameters  Used  To  Evaluate  L/D 


Constant  Kinetic  Energy  of  34  kJ 


V  =  1.25  krrVs 


V  =  1.0  krrVs 


V  =  0.75  krrVs 
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Penetrator  L/D  Study 


Overburden  Test  Setup 
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Overburden  Test  Results 
For  Residual  K.E. 


Liner  Edge  Material 


Residual  Penetrator 


Penetrator 
L/D=  1 


Overburden  “A”  (Sand) 
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Overburden  Test  Results 
For  Residual  K.E. 


PenetratorX  Overburden 
L/D=1  \  “D”  (Sand) 

Edge  Material 


Residual  Penetrator 


Make  Screens  that  triggered  X-Rays 


Overburden  Test  Results 
For  Residual  K.E. 


Make  Screens  that  triggered  X-Rays 
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Warhead  Test  Summary 

"  Effects  of  L/D” 


Over- 

burden 

(sand) 

Threat 

81mm 

120mm 

155mm 

Warhead 

Design 

A 

D 

O" 

A 

D 

O" 

A 

D 

O" 

A 

D 

KE  34.  KJ 

L/D  =1 

X 

■ 

H 

K 

X 

X 

X 

L/D  =  1.5 

X 

■ 

□ 

L 

X 

X 

L/D  =3 

X 

x 

□ 

X 
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Prototype  Warheads  Evaluated 


Single  (1) 
Penetrafeor 
Warhead 


Seven  (7) 
Penetrafeor 
Warhead 


Nineteen  (19) 
Penetrafeor 
Warhead 


Sixty-One  (61) 
Penetrafeor 
Warhead 
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Typical  Test  Setup 


AMC 


CTH  Simulations 

Simulation  vs.  Test  Data 


DETONATION 
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Penetrator  K.E.  =  48.5  KJ 
Overburden  =  A 


155mm  Round 
Behind  Overburden 


Residual  155mm 


155mm  Round 
with  Base  Plate 


t  A 


hir 


Frontal  View 


Residual  of  Base  Plate 


DEFLAGRATION 


AMC 


K.E.  =  34.1  KJ 


No  Overburden 


120  mm  Round  Setup 


1 20  mm  Round  Residual  1 20mm 


K.E.  =  34.1  KJ 
Overburden  =  A 


8 1  mm  Round  Deflagrated  Round 

Frontal  View  Setup 


Residual  81mm 


155  mm  Round  Setup  Damaged  Round 


Penetrator  K.E.  =  48.5  KJ 
“D”  Overburden 


8 1  mm  Round  Setup 


Damaged  Round 
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Overburden  Test  Fixture 


155mm  Post-Test 
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Prototype  Warhead  Test 

Summary 


Overburden  (Sand) 

Threat 

81mm 

120mm 

155mm 

Penetrator  L/D 

KE  (K|  ) 

A 

D 

O" 

A 

D 

O" 

A 

D 

O" 

A 

D 

1 

29.6 

X 

X 

34.9 

X 

X 

48.5 

X 

X 

X 

X 

3 

26.0 

X 

3+ 

34.1 

X 

X 

X 

X 

X 

MEFP  Warhead  Test  Summary 

Kinetic  Energy  Vs  Target  Reaction 


No  OB 

OB  A 

OB  D 

• 

• 

o 

• 

• 

1  1  C tatfe 

I  1  Dferde 

II  Brag; 

I  |  NbRattn 

•  L/D=1 

A  L/D=1 .5 

-  L/D=3 

▲  • 

—  A  o 

•  * 

—  A  o 

• 

—  o 

O 

o 

u 

o 

o 

81 


120 


155 


81  w/A  1 20  w/A  1 55  w/A 


81  w/D  1 20  w/D  1 55  w/D 


AMC 


Summary 


•  Established  MEFP  warhead  design 
parameters  required  to  cause  various 
reactions  against  mines  and  lEDs. 

-  MEFP  mass  &  velocity 

-  MEFP  L/D 


Demonstrated 
to  neutralize  m 


prototype  warheads  that  can  be  used 
ne  and  lEDs. 


•  Results  from  this  study  supports  two  Army  S&T 
programs 

-  Mine  Neutralization  Program 

-  Extended  Area  Protection  System  (EAPS)  Program 


Prevention  of  Sympathetic  Detonation  between 

Reactive  Armor  Sandwiches 


A.  Holzwarth 


22nd  International  Symposium  on  Ballistics 
Vancouver  BC,  November  14-18,  2005 


Fraunhofer 


Institut 

Kurzzeitdynamik 

Ernst-Mach-lnstitut 


ontents 


Introduction  and  objectives 

Experiments 

>  Test  set-ups 

>  Materials 

>  Tests  with  foams 

>  Tests  with  layers 

>  Flash  x-ray  photographs 

>  Experimental  results 

Summary  and  conclusions 
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Linctioning  of  Reactive  Armor  Sandwiches 


Reactive  Armor  Sandwich 


Base 

Armor 
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ubjects  of  Investigations 


Appropriate  materials 

-  Absorption  of  kinetic  energy 

-  Damping  of  shock  waves 

Useful  and  necessary  densities 
Arrangement  of  materials 
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?st  Set  -  ups:  Tests  of  Foams  and  Layers 


am  between  2  ERA-SWs 

SW  2  SW  1 


4 


Foam 


D 


V 


\ 


FP  3  FP  2  FP  1 

X-Ray 


Layer  on  Acceptor-SW 

SW  2  SW  1 


E 

03 

O 

H — 

o 

H — ' 

CO 


50 


\ 


7*\ 


FP  A U FP  3  FP  2  FP  1 


X-Ray 


> 


Foam-Layer  on  Acceptor- 

SW  SW  2  SW  1 
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xperiments  1 :  Foams  between  2  ERA-Sandwiches 


SW  2 


SW  1 


Foam 


D 


V 


]< 


FP  4  FP  3  FP  2  FP  1 


X-Ray 
3  Channels 


Test  set-up  data 


Sandwiches:  <2  /  4  /  2>,  1 60  x  40  mm2 


Flyer  plates:  Mild  steel,  2  mm  dick 

Explosive:  Seismoplast  (PETN  +  binder) 

4  mm  thick;  42,5  g 


Distance  D:  D  =  50  mm  -  80  mm 

Standard  value:  D  =  50  mm 
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ested  Foams 


>S  28,  p=28  kg/m3 
>S  50,  p=50  kg/m3 

>P  50,  p=50  kg/m3 

*P  300,  p=255  kg/m3 


J  -  Foam,  p=330  kg/m3 


p  <  400  kg/m3 


uminum-Honeycomb,  p=152kg/m3 
hite  Foam,  p=137  kg/m3 
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bservations  with  flash  x-ray 


Hocity  of  flyer  plate  FP  4: 

o  Strength  of  reaction  in  SW  2 

-  Detonation  /  Partial  Detonation 
Deflagration  /  No  Reaction 
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SW  1 


V 


Foam 


D 


FP  4  FF 


FP  3 


l< 


FP  2  FP  1 


X-Ray 
3  Channels 
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reliminary  Tests:  No  Material  between  Reactive  Sandwiches 


SW  2  reactive  -  detonation 


i/V  2  SW  1 


\ 


X-Ray 
3  Channels 


120  jus 


80  jis 


1 00  jus 


SW  2  inert  -  no  reaction 
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lash  X-Ray  Photographs:  Foams 

EPS  28  -  D=50  mm 


1390  m/s 


80  jis 


140  \is 


200  jis 
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:lash  X-Ray  Photographs:  Foams 


EPP  50  -  D=50  mm 


80  JlS 


120  JiS 


1 60  jus 
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White  Foam  -  D=50  mm 


80  [iS 


200  JLlJ 
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:lash  X-Ray  Photographs:  Foams 

EPP  300  -  D=50  mm 


80  JlS 

120  \XS 
1 60  jis 
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PU-Foam  -  D=50  mm 


100  \XS 


1  60  [iS 


220  jus 
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(FP  4)  versus  Areal  Density  mF  of  Foam  Fillings 


c/) 

E 


SW  2 


SV 


0 


10 


15 


20 


25  30 

mF  [kg/m2] 
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xperiments  2:  Layers  on  Acceptor-Sandwich 


SW  2 


SW  1 


03 

o 

H — 

o 

CO 


50 


\ 


7~\ 


FP  4  FP  3  FP  2  FP  1 

^  X-Ray 
3  Channels 


l< 
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Distance  between  SWs:  D  =  50  mm 
Thickness  of  layer:  d  =  variable 
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ested  Layers 

ibber  (Perbunan),  p=1 450  kg/m3 


lyurethane  (PU),  p=1 260  kg/m3 


ganic-plastic  material  (OGK),  p=1 100  kg/m3 


IP,  p=  1 950  kg/m3 


vlar,  p=  1 200  kg/m3 
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:lash  X-Ray  Photographs:  Layers 


Rubber 
6=5,6  mm 


--  - _ _ 

^/FP4:  800  m/s 

rlash  times: 

30,  140,  190  jas 
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PU 

d=10  mm 


VFP4:  585  m/s 

Flash  times: 

80,  140,  200  ns 


OGK 

d=14  mm 


i 

I 


VFP4:  490  m/s 

Flash  times: 

120,  160,  200  jas 
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PU 

d=20  mm 


\ 


V 


OGK 

d=22  mm 


VFP4:  350  m/s 

Flash  times: 

80,  140,  200  jus 


VFP4:  265  m/s 

Flash  times: 
140,  180,  220  ja: 
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(FP  4)  versus  Areal  Density  mF  of  Layers 


SW  2  SW 


mF  [kg/m2] 
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omparison:  v  (FP  4)  versus  Areal  Density  mF  of  Layers  and 


aams„ 

C f) 

E 

> 


SW  2  SW  1 


mF  [kg/m2] 
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periments  3:  Foams  as  Layers  on  the  Acceptor  Sandwiche 


SW  2  SW  1 


-P  4  FP  3  FP  2  FP  1 


Tested  foams:  PU-foam,  EPP  03 


Distance  between  SWs:  D  =  50  mm 
Thickness  of  layer:  d  =  5  mm  -  30  mm 


X-Ray 
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overview:  Foams,  Layers,  Foam-Layers 


SW  2 


SW  1 


”  1400 
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“  1200 


1000 
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ummary 

rhe  experiments  didn’t  exhibit  a  clear  distinction  between  detonation  and  no  reaction. 

3orous  foams  (p  <  400  kg/m2)  are  more  effective  than  layers  (p  >  1 1 00  kg/m2)  without 
Dorosity,  regardless  of  the  way  they  are  placed  between  reactive  sandwiches. 

rhe  areal  densities  necessary  for  a  certain  prevention  of  sympathetic  detonation  were  : 

-  In  case  of  foams:  Areal  density  >8-10  kg/m2 

-  In  case  of  layers:  Areal  density  >20-25  kg/m2 

-or  the  prevention  of  sympathetic  detonation  not  only  the  density  of  a  material  is  important  I 
also  other  material  properties,  e.g.  porosity. 
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End  of  Presentation 


2005  Fraunhofer  EMI 


Fraunhofer  |nstHut 

Kurzzeitdynamik 

Ernst-Mach-lnstitut 


Seite  24 


A.  Holzwarth 


uture  Work 


^/lore  accurate  investigation  of  the  material  properties  important  for  the  avoidance  of 
sympathetic  detonation 


Decrease  of  protection  power  of  reactive  armor  caused  by  damping  materials 


nvestigations  with  larger  reactive  armor  systems 
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QinetiQ 


The  influence  of  sabot  threads 
on  the  performance  of  KE 
penetrators 


Nick  Lynch,  John  Stubberfield:  QinetiQ  Ltd,  Fort  Halstead, 
UK 

22nd  ISB,  Vancouver,  November  2005 


Introduction 


•  The  majority  of  fin  stabilised,  kinetic  energy  (KE)  projectiles 
use  threads  along  the  interface  with  the  sabot  to  launch  the 
penetrator  from  the  gun 

•  The  threads  are  generally  undesirable  at  impact  on  a  target 
since  the  thread  root  forms  a  stress  concentration 

•  If  the  number  of  threads  could  be  reduced,  would  this 
improve  penetration  performance  ? 

•  Are  threads  needed  in  hydrocode  simulations  of  impact 
events  and  a  possible  cause  of  discrepancies  between 
experiment  and  simulation  ? 


QinetiC 


Scope  of  the  work 


•  Forward  ballistic  tests  (40mm  calibre) 

-  Four  designs  of  L/D  15  penetrator 

-  Two  types  of  multi-plate  target 
-1 600  m/s 

•  Reverse  ballistic  tests  (40mm  calibre) 

-  Two  designs  of  L/D  30  penetrator 

-  Oblique  plate  target  fired  at  pitched  attitude  penetrators 

-  1 650  m/s 


Projectiles 


Plain  finish  and  full  thread 


Projectiles 

•  h 

Half  thread  and  double  thread 

c  55  „  , 

45  PELL  THREAD 

36  THREADS  PER  25.4mm 

OO 

"St 

i 

f 

1 

. :TrTrTZ> 

^  10  90  FUL’L  THREAD 

36  THREADS  PER  25.4mm  ; 
(ONE  THREAD  HOT  CU 

l  START 

T) 

GO 

{ 

f 

-  - -  -  - - 

ze . : . ::: . 

QinetiC 

Target  designs 


1  plate,  t=8mm  RHA 
up  to  13  plates,  t=4mm  RHA 


V=1600  m/s  nominal 


Results  against  Target  1 


•  3  rod  types  tested 

•  Assessment  of  results 
made  difficult  by  variation 
in  impact  pitch  angle 

•  The  results  can  be  ranked 
by  pitch 

•  Allowing  for  this,  no 
apparent  difference  in 
penetration 


14 

12 


Total  penetration  (mm) 


Penetration  into  Target  2 


Dble  thread  ■  Full  thread  ▲  Half  thread  x  Rain 


Impact  KE(KJ) 

•  Unthreaded  rods  with  highest  energy 
went  no  deeper  than  other  designs 

•  1 61 5  m/s  unthreaded  rod  has  5%  greater 
KE  than  full  thread  design  at  1633  m/s 


•  Impact  pitch  less  than  0.5° 

•  All  except  one  result  perforated  6  plates 

•  Need  to  compare  line  of  sight  penetration 


Unthreaded  1/2  threaded  Fully  Double  start 

threaded 


QinetiC 


Average  crater  width 


•  Crater  width  reduces  due 
to  projectile  deceleration 

•  Crater  width  for  un¬ 
threaded  rods  increased 
from  plate  1  to  plate  2  - 
widest  craters  in  most  of 
the  plates 

•  Full  thread  rod  tends  to 
have  a  narrower  crater 


Crater  width  (mm) 


0  1  2  3  4  5  6  7 

Hate  number 


QinetiC 


How  much  change  in  penetration  could  be 
expected  ? 


•  If  the  threaded  rods  lost  one  pitch  per  plate  due  to  shear  at 
break-out,  what  effect  would  this  have  on  penetration  ? 

•  This  was  assessed  using  an  analytical  penetration  model, 
deleting  part  of  the  rod  at  plate  exit 

-  Nil  deleted  (plain  rod) 

-  0.7mm  deleted  (standard  thread) 

-  1 ,4mm  deleted  (double  pitch) 

•  What  effect  could  be  expected  just  from  the  difference  in 
effective  rod  diameter  ? 


QinetiC 


Penetration  depth  (mm) 


Penetration  vs.  time  for  rod  loss  options 
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Predicted  effect  of  rod  loss 


Rod  type 

Impact  Velocity 

(m/s) 

Rod  diameter 
(mm) 

Rod  loss  per  plate 
(mm) 

Total  penetration 
(mm) 

Unthreaded 

1625 

8 

Nil 

114 

Vz  thread 

1625 

8 

112  (interpolated) 

Full  thread 

1625 

8 

0.7 

110 

Double  thread 

1625 

8 

1.4 

106 

Full  thread 

1625 

7.70 

Nil 

112.85 

Double  thread 

1625 

7.77 

Nil 

112.8 

QinetiC 


1 .2  mm  change  in  penetration  predicted  due  to  effective  diameter 
8  mm  change  in  penetration  predicted  due  to  pitch  loss 
8  mm  difference  would  be  observed.  No  evidence  that  this  is  occurring 


Comparison  of  X-rays  -  Target  2 


Unthreaded  rod 


Full  thread  rod 


QinetiC 


Reverse  ballistic  experiments 


93%  Tungsten  alloy 

RHA  target  fired  projectiles 

at  1 650  m/s 


Projectiles  pitched  at  4° 


Dimensions  in 
millimetres 

QinetiC 


Comparison  of  L/D  30 
unthreaded 


Threaded  region 


Unthreaded 

rod 


VS. 


Conclusions 


•  Four  variants  of  L/D  1 5  threaded  penetrator  showed  no 
significant  difference  in  penetration  depths  against  two 
multiple  plate  targets 

•  In  contrast  there  was  a  marked  difference  in  the  fracture 
behaviour  of  L/D  30  pitched  attitude  rods  with  and  without 
threads 

•  Conclude  that  representing  threaded  rods  with  plain 
surfaces  in  simulations  is  valid  for  multiple  plate  targets  but 
not  for  more  disruptive  targets 


QinetiC 


QinetiQ 
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Introduction 


Grid-fin  configurations  offer  interesting  alternative 
to  classical  fin  designs  : 

-  Low  Hinge  Moment, Easy  storage,  Good  performance  at 
high  AoA. 

-  Drawback  :  Higher  drag  penalty 

Experimental  studies  conducted  at  DRDC  (in 
collaboration  with  ISL)  demonstrated  an 
aerodynamic  choking  phenomena  on  two 
configurations  : 

-  Thick  fin  model :  flow  is  choked  over  a  large  range  of 
Mach  numbers 

-  Thin  fin  model :  flow  choking  occurs  at  specific  Mach 
numbers 

Choking  effect  has  not  been  reproduced  by  CFD  or 
wind  tunnel  tests. 
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Introduction 


1.000  1.20G  1400 
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Figure  2.  Zero  yaw  drag  coefficient. 
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Figure  3.  Pitch  moment  coefficient. 
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Figure  4.  Normal  force  coefficient.  Figure  5.  Center  of  pres  sure  location. 

Choking  phenomena  measured  in  aeroballistic  range  (Eglin  AFB) 
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Model  Configuration  (Thin  fins) 


a.  aa 


All  dimensions  in  calibers 


Air  Force  Finner  body,  with  four  grid  fins 
9  grid  cells  per  fin. 

No  cant  angle  for  fins 
1  caliber  =  20mm 
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11,14 
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Numerical  Modeling 

L^ 


•  CFD  work  done  with  ANS YS  CFX  5.7 

-  Navier-  Stokes  equations 

-  2nd  order  advection  scheme 

-  k-epsilon  turbulence  model  was  used 

•  Hexahedral  mesh  built  with  ANSYS  ICEM  CFD  Hexa 


-  1.6  million  elements  for  supersonic  flow  domain 

-  2.1  million  elements  for  subsonic  flow  domain  (with  base  flow) 
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Numerical  Modeling 


CFD  computations  performed  at  following  conditions  : 

-  Baseline  Mach  no.  0.6,  0.8,  1.5,  2.0,  2.5,  3.0 

-  Subsonic  runs  conducted  at  small  increments  between  Mach 
0.7  and  0.85  for  flow  choking  simulations 

-  Supersonic  runs  conducted  at  small  increments  between  Mach 
1.8  and  2.5  to  explain  sudden  change  in  aerodynamic 
coefficients. 


All  runs  made  a  2°  angle  of  attack 

Cx,  Cn,  Cm  and  Xcp  aerodynamic  coefficients  computed 


Defence  R&D  Canada  •  R  &  D  pour  la  defense  Canada 

^  SNC  TEC 


Results  -  Normal  Force  Coefficient 


Peak  in  CNa  at  M=0.78. 
Wind  Tunnel  results  do  not 
reproduce  this  peak. 

Theoretical  Choking  Mach 
for  fin  cells  is  M=0.744. 

Subsonic  choking  captured 
in  CFD  results. 

Discontinuous  variation  of 
CNa  between  Mach  1.8 
and  2.4. 

Good  agreement  between 
CFD  and  Wind  Tunnel 
results. 
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Results  -  Pitch  Moment  Coefficient 


Similar  peak  at  M=0.78. 
Wind  Tunnel  results  do  not 
reproduce  this  peak. 

Effect  of  subsonic  choking 
on  CMa  is  to  increase 
stability. 


Again,  discontinuous 
variation  of  CMa  between 
Mach  1.8  and  2.4. 

Good  agreement  between 
CFD  and  Wind  Tunnel 
results. 
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Results  -  Axial  Force  Coefficient 


CA  vs.  Mach  -  CFD  and  Experimental  Results 
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Good  agreement  between 
CFD  and  Wind  Tunnel 
results. 

CFD  results  in  supersonic 
regime  had  to  be  corrected 
for  base  drag. 
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Results  -  Center  of  Pressure 


ao  0.5  LO  L  5  ID  2.5 

Mach  na 


3.0  3.5  4.0 
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Analysis  -  Subsonic  choking 


CNa  contribution 
breakdown :  peak  effects 
come  from  fins. 

Peak  value  at  Mach  0.78  due 
to  choking  effects. 

Confirmed  by  apparition  of 
Mach  waves  inside  fin  cells. 
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J\l_y  Analysis  -  Subsonic  choking 
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a)  Freestream  Mach  =  0.72 
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b)  Freestream  Mach  =  0.75 
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Contour  plots  of  sonic 
regions  between  Mach 
0.95  and  1.05. 

As  Mach  number 
moves  from  0.72  to 
0.78,  normal 
shockwaves  appear  in 
cells. 

Completely  choked 
state  reached  at  Mach 
0.78. 

Good  correlation 
between  theoretical 
choking  Mach  number 
and  CFD  results. 
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f  hj  Theoretical  Grid  Fin  Flow  Field  Model 

L^ 


MACH  NUMBER 

^NF<*  Versus  Mach  Number, 


CEILWM1 

*.  0OtJNOAflY  U-  *°RW<lL 

layer  *  Smock 

*  \ _ 

rtjefrrrrtrtxi 

a)  GhofcwJ  Flow  {M  <1} 


c)  Shock  Mtactunsnl 
andRefteaton 


f^rnasm 


tf&rri 


b)  Choked  Fk>w  (M  >  i) 


Grid  Fin  Flow  Field. 


Ref.  Washington  et  al. 
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Analysis  -  Supersonic  choking 


CNa  vs.  Mach  ■  Supersonic  Regime  (2deg  AoA) 


6.0 
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CNa 
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A 

CNAb 
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.  1 

2.5 

3.0 

3.5 

CNa  contribution 
breakdown :  peak  effects 
come  from  fins. 

Discontinuous  behavior 
between  Mach  1.8  and  2.4. 

Change  in  fin  shock-wave 
configuration  explains 
unusual  variation  in 
coefficients. 


Mach 
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Analysis  -  Supersonic  choking 


M  =  2.30 


Mach  Number 

(Contour  Flrrj  Mach) 
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M  =  1.80 


Contour  plots  of 
sonic  regions 
between  Mach 
0.95  and  1.05 


•  Flow  completely  choked  at  Mach  1.5,  with  transition 
to  complete  “unreflected”  state  at  Mach  2.3. 

•  Irregular  behavior  between  Mach  1.8  and  2.40  due  to 
flow  topology  change  from  fully  choked  (shock 
outside  of  grid  cells)  to  “non-reflecting”  states. 
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Conclusion 


•  A  detailed  CFD  analysis  of  flow  inside  grid  fins 
gave  insight  on  non-conventional  behavior  of  the 
main  aerodynamic  coefficients. 

•  Flow  choking  effects  successfully  predicted  in 
subsonic  and  supersonic  flow  regimes. 

•  Subsonic  choking  occurred  at  a  specific  Mach 
number.  Important  offset  between  wind  tunnel  and 
CFD  results  in  subsonic  regime. 

•  Supersonic  flow  transitions  from  fully  choked  to 
non-reflecting  states  over  a  large  range  of  Mach 
numbers. 
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Conclusion 


•  Theoretical  flow  choking  model  by  Washington 
(1993)  demonstrated  by  CFD  calculations. 

•  Experimental  work  Eglin  AFB  revealed  loss  of 
stability  at  subsonic  choking  conditions.  With  the 
present  configuration,  effect  is  opposite  :  gain  in 
stability. 

•  CFD  proved  to  be  an  essential  tool  in  the  study  of 
this  complex  flowfield. 
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Normal  force  coefficient  derivative, 


Figure  1  Generic  tail  control  missile  (GTCM)  with  grid  fins 
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Introduction 

♦  Projectile:  Lead/Copper  Partitioned  (A-Frame) 
Hunting  Bullet 

♦  .338  Winchester  Magnum 

♦  Chronographed  Muzzle  Velocities 

842.16  m/s  (2763  ft/s)  Mild  Steel  Impacts 
854.35  m/s  (2803  ft/s)  Armor  Impacts 
16  tests:  mean  =  847.98  m/s  (2782  ft/s) 

std  dev  =  ±7.44  m/s  (24.4  ft/s)  =  ±0.877  % 

♦  Witness  Plate:  Mild  Steel 

♦  Armor:  Kevlar®  and  Kevlar® /Steel 


Lead 


Copper 


32.5  mm 
(1.28  in) 


yC 


Copper 


ji/ 


1  Equation  of  State 

Shock 

I  Reference  Density  (g/cm3) 

8.93  1 

1  Gruneisen  Coefficient 

1.99  1 

1  Parameter  C1  (m/s) 

3.94E03  | 

1  Parameter  S1 

1.489  1 

1  Strength  Model 

Von  Mises 

1  Shear  Modulus  (KPa) 

4.5E7  1 

[  Yield  Strength  (KPa) 

7.0E4 _ 

■*_  m.  _  J  A  ■'TL, 

■I  . ■ 

^^9 

-  \  f  r  -  p  r-~  7 

?  i  V 

Equation  of  State 

Shock 

Strength  Model 

Johnson-Cook 

Reference  Density  (g/cm3) 

7.896 

Gruneisen  Coefficient 

2.17 

Parameter  C,  (m/s) 

4.569E03 

Parameter  S1 

1.49 

Reference  Temperature  (K) 

300 

Shear  Modulus  (kPa) 

8.18E07 

Yield  Stress  (kPa) 

5.17106E05 

Hardening  Constant  (kPa) 

2.75E05 

Hardening  Exponent 

0.36 

Strain  Rate  Constant 

0.022 

Thermal  Softening  Exponent 

1.0 

Melting  Temperature  (K) 

1.81 1E03 

Lead/Copper  Bullet  Penetration  of  Mild  Steel 
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Color  Represents  Pressure 


AUT0DYN-2D  Version  4.3.01a 


Century  Dynamics  Incorporated 


MATERIAL 

LOCATION 


LEAD 
COPPER 
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VOID 
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X 


Sea  I  e 
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BULL-B;  225  GR  0,339  CAL  FACTORY  HUNTING  BULLET 


Penetration 
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Color  Represents  Absolute  Velocity 
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Century  Dynamics  Incorporated 


MATERIAL 
LOCAT I  ON 


LEAD 
|  COPPER 
1006  STEEL 
VOID 
Y 

A 

X 

- ? 

Eca  I  e 

2.3Q0E+01 
AX  (rrm.mg  .ms) 


BULL-B;  225  GR  0,333  CAL  FACTORY  HUNTING  BULLET 
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Where: 

1 

d 

1 71 

t 

e 


=  Projectile  length 
=  Projectile  diameter 
=  Projectile  mass 
=  Target  thickness 

=  Impact  Angle 


=  3.25  cm 
=  0.8585  cm 
=  14.578(g) 

=  1.427  cm 

=  0.0  deg 


V,  =  939.12- 

1  s 

Compared  to  848  m/s  =>  almost  penetrates 


Ref:  Introduction  to  Terminal  Ballistics  -  Course  Notes,  Donald  R.  Carlucci,  2004 


Range  of  mixed  results 


600  800  1000  1200 
Striking  velocity  (m/s) 


Ref:  Impact  Dynamics,  Zukas  et  al 
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Kevlar®^ 
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Equation  of  State 

Puff 

Reference  Density  (g/cm3) 

1.29 

Parameter  A1  (kPa) 

8.21  E06 

Parameter  A2  (kPa) 

7.036E07 

Parameter  A3  (kPa) 

0.0 

Gruneisen  Coefficient 

0.35 

Expansion  Coefficient 

0.25 

Sublimation  Energy  (J/Kg) 

8.23E06 

Parameter  T1  (kPa) 

0.0 

Parameter  T2  (kPa) 

0.0 

Reference  Temp  (K) 

0.0 

Specific  Heat  (C.V.)  (J/kgK) 

0.0 

Strength  Model 

Von  Mises 

Shear  Modulus 
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Lead/Copper  Bullet  Penetration  of  Kevlar®  and  Mild  Steel 

Color;  Represents  Absolute  Velocity 

AUT0DYN-2D  Version  4.3.01a  Century  Dynamics  Incorporated 

MATERIAL 
LOCAT I  ON 

LEAD 
COPPER 
1006  STEEL 
KEVLAR-VM 
VOID 


Eca  I  e 

2 , 6GQE+01 
AX  (rrm.mg  .ms) 


B-K-VM;  225  GR  0,333  CAL  FACTORY  HUNTING  BULLET 


a.  6  jjls 


a.  76  jjls 

* 


b.  88  jjls 


A05 , VEL 
(rn/s) 


IS,50E+02 
7,fS5E+02 
e.aoE+02 
5.95E+02 
5, 10E+02 
4.25E+02 
3.40E+02 
2,55E+02 
1 , 70E+02 
S,50E+01 
0 , OOE+OO 


■  '  Iv-J 

ii  v'i'Uii  _d. 


AUT0DYN-2D  v5.0  from  Century  Dynamics 


Material  Status 


Elastic 


kevlar 
Cycle  0 

Time  0.000E+000  ms 
Units  mm,  mg,  ms 
Axial  symmetry 


225  V  ]  'BB  CAL  FACTORY  HUNTING  BULLET 
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Kevlar  Model 


225  Grain  0.338  Cal  Hunting  Bullet 
Kevlar  Board 
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Color  Represents  Absolute  Velocity. 


Reasonable  Bullet  Erosion 


Lead/Copper  Biillet  Penetration  of  Kevlar®  and  Mild  Steel 
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Conclusions 

♦  Hydrocode  and  ballistics  lab  results  match  well  for 
lead/copper  A-Frame  bullets  impacting  on  mild  steel. 

♦  A  reasonable  match  between  hydrocode  predictions 
and  ballistics  lab  results  is  obtained  for  Kevlar®  armor 
using  a  Von  Mises  strength  model.  This  “engineering” 
model  was  used  to  complete  the  armor  design. 

♦  A  Von  Mises  strength  model  misses  some  of  the 
essential  physics  of  bullet  penetration  into  a 
transversely  isotropic  material  like  Kevlar®. 

<►  Use  of  a  Lagrangian  transversely  isotropic  material 
model  yields  more  consistent  results  for  the  Kevlar. 
However,  modeling  the  lead/copper  bullet  is  more 
appropriate  in  the  Eulerian  Frame  of  Reference. 

♦  Simulations  using  an  Eulerian  Bullet  and  Lagrangian 
Kevlar  have  been  attempted  with  reasonable  success. 

The  results  are  dependent  on  the  Kevlar  erosion  rate  selected 
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SCORPION  -wmm 

Self-CORrecting  Projectile  for  Infantry  OperatioN  ^ 


GOAL:  Demonstrate  a  Guided  Spinning  Projectile  Using  MAFC  Technology 


2. 


3. 


Pacing 

Technologies: 

'/MAFC  actuators 

'/Flow  control  concept 
for  spinning  projectiles 

'/Flight  control  algorithm 

'/Initialization  and  INS  for  spinning  projectile 

'Z Compact,  g-hardened  electronics  and  packaging 


Objectives-. 

1.  Demonstrate  microadaptive  flow 
control  (MAFC)  authority  and 
guidance  algorithm  for  a  medium- 
caliber  munition  at  subsonic  speeds. 

Provide  a  suite  of  validated 
advanced  design  tools. 

Establish  technology  transitioning 
pathways  for  tactical  systems. 


'/Design  tools:  integrated  computational  fluid  dynamics 

(CFD)  and  flight  dynamics 


Micro-jet  CFD  Flow  Visualization 

U,  =  37  m/s,  a  =  0°,  Ujet  =  31  m/s,  f=  1000  Hz,  no  spin 


Lift 

A 


0  10 


Investigate  “placement”of  jets  for 

•  lift  control  to  improve  accuracy 

•  side  force  for  divert  authority 


30  40 


100950  101250 


Asymmetric 

Flow 

Separation 


Flow 

Turning 

mv 


Flow 

Separation 


Particle  Traces  Visualization 

U  =  82  m/s,  a  =  0°,  Ujet  =  31  m/s,  f  =  1000  Hz,  Spin  =  67  Hz 


Clip  Plane  -  Turbulence 

TIME=0.520300 


0.0 


Turbulence 

136*3 


272.6 


68,2 


204,5 


For  Two  Spin  Cycles 
1st  cycle:  Red 
2nd  cycle:  Blue 


Clip  Plane  -  Turbulence 

TIME=0.520300 


SCORPION  Open-Loop  Test 

Flight  Test  Results  Compared  With  Simulated  Impacts 


Experimental  Muzzle  Velocity  Distribution 
Used  for  Simulation 


Experimental 

Impacts 

^  Right  Divert 

,  Command 
^No 

Command 
^  Left  Divert 
Command 


ONLY  a  few  jet 
cycles  are 
available 


Cross  Range  (m) 

Activate  for  1/4  revolution  (about  4  jet  cycles)  such  that  force 
generated  will  be  horizontal  (left  or  right,  as  selected). 


Flight  Data  Confirms  Control  Force 


il  Georgia 


Control  Force 
Experimentally  Derived 
From  Flight  Test  D; 


Sabot/Pusher 

Assembly 


Projectile 


ARL  Unsteady  CFD 
(Spinning) 


GTRI  Wind  Tunnel  Experiments 
(Non-spinning) 


UNSTEADY  AERODYNAMICS/FLIGHT  DYNAMICS 

DOD  HPCMP  Grand  Challenge  Project 


GOAL:  Virtual  Fly-Out  of 
Projectiles 

Coupled  CFD/Rigid  Body 
Dynamics  (RBD)  Simulations  of 
in-flight  spinning  and  finned 
projectiles 


Guided  Medium  Caliber  Munitions 


r.\ 


•  Demonstrated  Micro-Adaptive  Flow  Control  for  divert  of  subsonic  guided  40  mm  grenade 

•  Demonstrated  Multi-disciplinary  physics  modeling  -  flew  munition  through  the  computer 
using  High  Performance  Computing 

•  First  divert  ever  of  a  spin  stabilized  munition  system  at  60  hertz  spin  rate 

•  Developed  a  miniature,  G  hard,  on  board  flight  control  system 

•  Demonstrated  initialization  at  muzzle  exit  -  Velocity  -  Orientation 

•  Demonstrated  open  loop  divert 

•  Demonstrated  closed  loop  guidance  to  the  target  on  major  error  source  -  Velocity 

•  Cut  on  target  dispersion  due  to  muzzle  velocity  variation  to  one  third  of  the  system  value 


Micro  Adaptive 
Flow  Control  Actuation 


Georgia 

1 


Experimentally  Demonstrated  Novel 
Aerodynamic  Control  Methodology 
Capable  of  Diverting  Medium  Caliber 
Munitions 


Technologies  developed  in  this  research  were  critical  in  the  design  of  flight  control 
systems  required  for  subsequent  flight  tests  that  demonstrated  the  ability  of  MAFC  to 
divert  the  trajectory  of  a  spinning  projectile  in  flight 


SCORPION  CFD  VIDEO 
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PTIVE  FLOW  CONTROL  AND 
LINEAR  AERODYNAMICS  ^ 


~^^*{60MBUSTION  GAS  GENERATORS) 

Dr.  Jubaraj  Sahu 
Ms.  Karen  Heavey 

U.S.  Army  Research  Labor 


GTRI,  Atlanta,  GA,  23  Sept  2005 


Weapons  <£  Materials  Research  Directorate 
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GAS  COMBUSTOR  J  ETS 


■■I 


Georgia  ^ifletacomp  fH 

J  ^TbchjiXE3U®@ja?  Technologies  M  M  m 


Demonstrate  Adequate 
Aerodynamic  Steering 
Forces  using  Flow  Control 


Objective:  Provide  fundamental  understanding  of  the  flow 
phenomena  associated  with  synthetic  jet  micro-adaptive  flow 
control  (MAFC)  and  assess  its  effectiveness  to  provide 
adequate  aerodynamic  forces  to  divert  or  guide  a  small  caliber 
projectile  to  its  target 


Pacing  Technologies: 

•  Combustor  gas  generators 

•  High  Performance  Computing 

•  Advanced  Visualization 

•  Unsteady  Aerodynamics 


Warfighter  Payoffs: 

•  Improved  lift  control 

•  Precision-guided  munitions 

•  Increase  lethality 


High  Fidelity  Computational  Tool  for  Improved  Performance  of  Army  Munitions 


COMBUSTOR  JET  CFD 


**ERIC4* 


•  3-D  Navier-Stokes  computational  technique 

•  Steady  jet  CFD 

•  Unsteady  (time-accurate)  jet  CFD 

•  Fast  convergence  to  steady  state 

•  Fast  computation  of  unsteady  flows 

•  Dual  time-stepping 

•  Special  jet  boundary  conditions  for  unsteady  CFD 

•  Validation  of  both  non-spinning  and  spinning  cases 
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U-Velocity  Contours 
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STEADY  JET  CFD 

M  =  0.8,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 
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Velocity 

Contours 


STEADY  JET  CFD 

M  =  0.8,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Surface 

Pressure 

Contours 


Temperature 

Contours 


STEADY  JET  CFD 

M  =  0.8,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Body  jet-on,  Fy  =  -9.1  N 
Body  jet-off,  Fy  =  0.0 
Jet  Force,  Fjet  =  -3.8  N 


Amplification  Factor  =  (Fjet  +Fjj)  /  Fjet 
or,  AF  =  ( -3.8  +  (-9.1  -  0.0) )  /  (-3.8)  =  3.3 


COMBUSTION  JET  CFD 

Jet  Size  =  0.6  mm  axial, 

0.3  mm  circumferential 
Jet  location  =  31  mm  ahead  of  the  step 


STEADY  JET  CFD 

M  =  0.8,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Pressure 

Contours 


Temperature 

Contours 


STEADY  JET  CFD 

M  =  0.8,  a  =  0°,  P0  =  20  atm,  T0  =  4000  K 


Pressure 

Contours 


T913e-n2 


Temperature 

Contours 


COMBUSTION  JET  CFD 

80mm 


STEADY  JET  CFD 

U  =  30  m/s,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Pressure 

Contours 


Temperature 

Contours 


STEADY  JET  CFD 

U  =  30  m/s,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Experiment 


CFD 


STEADY  JET  CFD 

U  =  30  m/s,  a  =  0°,  P0  =  5  atm,  T0  =  2000  K 


Body  jet-on,  Fy  =  0.06  N 
Body  jet-off,  Fy  =  0.0 
Jet  Force,  Fjet  =  -0.08N 


Amplification  Factor  =  (Fjet  +Fjj)  /  Fjet 

or,  AF  =  ( -0.08  +  (0.06  -  0.0) )  /  (-0.08)  =  0.25 


The  use  of  electric  power  in  active 
armour  applications 

IBS  2005 
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P  reservation  Overview 


•Vehicle  survivability 


•  Countermeasure  with  pulsed  power  supply  and 
directional  control 


>  DAS 


•  S  ensor  development  I 


•Conclusion 
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P  reservation  Overview 


•Vehicle  survivability 


•Countermeasure  with  pulsed  power 
supply  and  directional  control 

J  >  DAS 


•  S  ensor  development  J 

•Conclusion 
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Advantages  of  EM  launched  plates 


On/Off  switching  capability 

No  energetic  material  on  the  outside  of  the  vehicle 

•  Directional  launching 

•  Multi-hit  capability 

Possibly  effective  against  CE  (shaped  charge)  and  KE 
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EM  launched  plate  (3-D) 
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EM  launched  plate  (2-D) 


Plate  velocity:  150-200  rrVs 
Plate  mass:  0.5  kg 
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EM  launched  plates  vs  KE: 

static  experiments 
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Effectiveness  against  30*173  mmAPFSDS 
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EM  launched  plates  vs  KE: 

numerical  simulations 
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EM  launched  plates 
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Launched  plates 


Plate  velocity:  50  rrVs 
Plate  mass:  3.0  kg 
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X-ray  and  camera  images 


Reduction  of  DOP: 
>  60% 
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P  reservation  Overview 


•Vehicle  survivability 


•  Countermeasure  with  pulsed  power  supply  and 
directional  control 

r  DAS 


•  Sensor  development 

« Conclusion 
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Conclusion 


Countermeasure  steering 
•  3-D  possible 

•  Countermeasure  effectiveness 

Initiates  premature  detonation  of  shaped  charge 
Break  up  or  erosion  of  KE :  Decreasing  DOP 

•  Promising  sensor  design 

•Accuracy  of  3D  positioning  and  timing  is  a  real 
challenge 
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P  reservation  Overview 


•Vehicle  survivability 


•  Countermeasure  with  pulsed  power  supply  and 
directional  control 

•Sensor  development 

•  Sensor  design 

•  RPG-7  signatures  fromTNO  experiments 


'V 


r  DAS 
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•Conclusion 
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Experimental  set-up 
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From  Columbia  to  Discovery : 
Understanding  the  Impact  Threat 
to  the  Space  Shuttle 
(Abridged) 

James  D.  Walker 
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The  work  described  in  this  talk  was  paid  for  by  the 
American  taxpayer  through  NASA  and  the 
Columbia  Accident  Investigation  Board. 
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Lots  of  People 


Photograph  during  Columbia  investigation 


USA 


Lots  of  people  worked  on  the 
Columbia  accident  investigation 
and  return  to  flight. 

Some  who  directly  helped  with  the 
impact-related  material  in  this  talk: 

Donald  J.  Grosch  (SwRI) 

Sidney  Chocron  (SwRI) 

Walt  Gray  (SwRI) 

Justin  Kerr  (NASA/JSC) 

Freeman  Bertrand  (Jacobs/Sverdrup) 
Paul  Parker  (Boeing) 

Mike  Dunham  (Boeing) 
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Southwest  Research  Institute 
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The  Columbia  Accident 
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STS-107 


•  1 13th  flight  of  the  Space  Shuttle  Program 

•  28th  flight  of  Columbia 

•  Columbia  was  the  first  Shuttle  to  fly  (April, 
1981) 

•  Launched  January  16,  2003 

•  Disintegrated  on  re-entry,  February  1,  2003 

•  Crewed  by 

-  Rick  Husband 

-  William  C.  McCool 

-  Michael  P.  Anderson 

-  David  M.  Brown 

-  Kalpana  Chawla 

-  Laurel  Clark 

-  Ilan  Ramon 
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CAIB  Report,  Volume  I 


COLUMBIA 

J [UK Hi  IMtfl$i|S([!iH  H II I H □ 


R  efort  Volume  1 

A  .in  :JflT  2QQ3 


Volume  1  of  the  Columbia  Accident 
Investigation  Report  was  published 
August  26,  2003 

It  can  be  found  at  the  CAIB  web  site: 
www.caib.us 

The  conclusions  were  that  there  were 
two  causes  of  the  loss  of  the 
Columbia : 

-  A  physical  cause 

-  An  organizational  cause 

SwRI  was  involved  in  identifying  the 
physical  cause,  discussed  in  Chapter  3 
of  the  Volume  I  with  more  details  in 
Volume  II  and  a  NASA  report. 
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Columbia  Accident  Investigation  Board  \ 
Report,  Volume  II,  Appendix  D.12  (Oct.  2003) 
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Bipod  Ramp  Foam  Insulation 
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Foam  Impact  on  STS- 107 


STS-107  ET208  Video  Fie 


'W 


SfcKK) 


-mm 


-m s© 


16:15:40:21.692  to  .858  UTC 
Unable  to  track  debris  in  two  fields  (.775  &  .792) 


•  81.7  seconds  into  flight,  the  bipod 
ramp  foam  insulation  broke  away. 

•  The  shuttle  was  at  66,000  ft. 

•  The  shuttle  velocity  was  mach  2.46 
(1580  mph). 

•  The  foam  traversed  the  58  feet 
between  the  bipod  connection  and 
the  left  wing  in  0.16  seconds. 

•  The  refined  estimate  of  the  impact 
velocity  was  775  ft/s  (528  mph). 

•  Estimated  foam  mass:  1.67  pounds. 

•  Estimated  size:  21”  x  11.5”  x  5.5” 
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Possible  Impact  Locations 


Wing  Acreage 


Main  Landing  Gear  Door 
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Shifting  Focus 


Modular  Auxiliary  Data  System  Recorder 


Early  telemetry  showed 
first  anomalies  in  left 
main  landing  gear  bay. 

Upon  recovery  of  the 
Modular  Auxiliary  Data 
System  recorder  on 
March  19,  2003, 
analysis  of  the  data  led 
attention  to  the  leading 
edge. 
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Main  Landing  Gear  Door  from  Enterprise 
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SwRI  Large  Compressed  Gas  Gun 


500  gallon  tank,  275 
psi  working  pressure, 
10-inch  diameter  by 
35-foot  long  barrel. 


Typically  used  to 
launch  large, 
irregularly-shaped 
projectiles. 


Helium  or  Nitrogen  is 
used  as  the  driver  gas. 
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Launching  Foam 
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Main  Landing  Gear  Door  Tests 


•  Performed  5  impact  tests  on 
the  left  Main  Landing  Gear 
Door. 

•  These  impacts  demonstrated 
that  at  the  speed  (775  ft/s)  and 
angles  of  impact  on  the  tiles 
(up  to  13°)  the  foam  barely 
damages  the  tiles. 
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Tile  and  Foam  Modeling 
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Foam  Insulation  Compression  Tests 


Stress-strain  curve  for  sample  1 


Stress-strain  curve  for  sample  2 


Stress  (nsi ) 


Transverse  compression  tests  on  roughly 
2”  cube  tile  samples  taken  from 
tile  MISC-7 94-400- 1 20 


Failure  begins 


Crashing  begins  at  345  kPa  (50  psi) 
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Stress  (psi) 


New  EOS  Model  for  Foam  Insulation 


Enlargement  of  low  strain  region, 
blue  and  black  are  data,  red  is  model 


Full  stress-strain  curve,  blue  and 
black  are  data,  red  is  model 


Foam  density  is  0.03844  gm/cm3  (2.4  lb/ft3),  initial  crushup  stress  is  220  kPa  (32  psi) 
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Stress  (psi) 


New  EOS  Model  for  Thermal  Tile 


Strain  (%)  Strain  (%) 


Enlargement  of  low  strain  region, 
blue  and  black  are  data,  red  is  model 


Full  stress-strain  curve,  blue  and 
black  are  data,  red  is  model 


Tile  density  is  0.18  gm/cm3  (11.2  lb/ft3),  initial  crushup  stress  is  345  kPa  (50  psi) 
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Stress  (psi) 


Computations 
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LUO)  A 


Computations 


YY  Stress 

Stress  YY 
(d/cm1) 

I — I  2.WQKTO*  • 
I  1.747*10* 

1  1.494*10* 
1.241*10* 
9.673*10® 
7-342*10* 
4.B1Q*1Q5  # 

-2.&32x10* 
-2.7WK103 
M  -5.31&*1Q* 
-7.B4Bx1Q5  • 
1 1  -1030x10* 
1-1291*10* 
H-tS44KlOe 
I  -1797*10B 


-4  0  4  8  12 

3DR  X  (cm)  Z=0. 

Foam  insulation  into  tile  impact,  BOQfps,  30dag,  tl21  5/3/03(CTH) 
ECKDDY  5/03/03  11:49:25  CTH  351  Time=2.00713x10'  + 


Most  computations  were  performed 
in  2D  plane  strain  -  however,  there 
was  excellent  agreement  between  2D 
plane  strain  and  3D. 

The  normal  stress  was  examined 
along  the  surface  of  the  tile  gages 
(“tracers”)  spaced  0.5  cm  apart. 

The  gage  readings  show  the  normal 
stress  in  the  y  direction  at  the 
respective  location. 

Computations  performed  in  CTH 
with  new  material  models. 
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Y  (cm 


No  Damage  for  700  ft/s  Impact  at  15° 

Impact  Angle 


YY  Stress 


YY  Stress  versus  time 


-4 

2DR  X  (cm) 

Foam  Insulation  into  tile  impact,  700fps,  15deg,  tl16  5/2/03(CTH) 
EBTCFF  5/02/03  19:26:30  CTH 


641  Time=5. 00001x10 


-4 


TIME  fas) 

roam  insulation  into  tile  impact,  70Qfps,  15deg,  tl16  5/2/03(CTH) 
IBTCFF  5/02/03  19:26:30  CTH 


Tile  surface  mostly  flat  after  impact;  normal  stresses  are 
below  345  kPa  (50  psi)  tile  crushup  level 


500 
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wo)  A 


Crater  for  700  ft/s  Impact  at  23°  Impact 

Angle 


YY  Stress 


2DR  X  (cm) 

Foam  insulation  into  tile  Impact,  700fps,  23deg,  tl17  5/2/03(CTH) 
EBUBPU  5/02/03  20:19:48  CTH  638  Time=5.00619x10"4 


YY  Stress  versus  time 

6  0  0  - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - r 


_ 10  0  _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 1 _ 

0  100  200  300  400  500 

TIME  G*s) 

sam  insulation  into  tile  impact,  700fps,  23deg,  tl17  5/2/03(CTH) 
iLIBPU  5/02/03  20:19:48  CTH 


Crater  (indention)  seen  in  target;  normal  stresses  are  above  345  kPa 
(50  psi)  tile  crushup  and  400  kPa  (58  psi)  tile  yield/failure 
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Damage/No  Damage  Transition 
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Many  CTH  computations 
were  performed  using 
new  EOS  with  1”  cubes 
of  foam  insulation 
impacting  a  tile. 

These  computations 
allowed  a  determination 
of  where  the  transition 
from  no  damage  to 
damage  occurred,  in 
terms  if  impact  velocity 
and  impact  angle. 

Figure  shows  results  of 
computations. 


Inpact  Velocity  (ft/s) 


Theory 


•  It  is  possible  to  determine  the  impact  velocity  at  90  degrees  (flyer 
plate)  at  which  the  thermal  tile  just  begins  to  crush.  The  velocity 
is  given  by 


V crush  -  Uet  +  11  ef  + 


^ crush-t  ® crush- f 

Pef  (C1  /  Uef  ) 


where  a  is  stress,  c  is  sound  speed,  subscript  t  refers  to  tile, /to 
foam,  e  to  elastic,  and 


Ue  -  ° crush  1  A/'O  *  Pe  =  Po  /(  l  ~  Ue  1 Ce  ) 


When  computed,  Vcrush  =  68.2  m/s  (224  ft/s).  This  value  was 
confirmed  by  ID  CTH  computations. 
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Y  (cm 


Theory 


i  o 


YY  Stress 


8  - 


6  - 


-2 


-4 
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8 


Stress  YY 


(d/cm1) 
m  2,000*10* 

•  When  foam  impacts  the  tile  at 

1  1.747x10* 

1  1.494x10* 

an  angle,  the  horizontal  and 

1.241x10® 

downward  components  of  the 

9.S73x10G 

7.342x1Q5 

velocity  interact  nearly 

4.B10X105 

separately  if  the  interface 

-2.552x10^ 

remains  relatively  flat. 

-2.7Wx10a 

•  Thus,  crushing  of  the  thermal 

-5.316x10* 

-7.B4Bx10E 

-t03Bx10* 

tile  would  be  expected  when 

-1201x10* 

-1544*10* 

-1797x10b 

y  ^ crush 

1  2 


X  (cm) 

Foam  insulation  into  tile  impact,  BQQfps,  30deg,  tlx  5/2/03(CTH) 
ECNEGS  G  5/03/03  13:47:58  CTHGEN  0  Time=G. 


sin(0) 
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Damage/No  Damage  Transition 


When  the  computed 
Vcrash  =  68.2  m/s  (224 
ft/s)  is  used,  good 
agreement  is  obtained 
with  the 

computational  results. 
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Damage/No  Damage  Transition 


1999  foam  impact  data  with  1”  square  cross  section 
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Damage/No  Damage  Transition 


Five  impact  tests  during  Columbia  investigation  (minimal  damage) 


Damage/No  Damage  Transition 


400  800  1200  1600 

Impact  Velocity  (ft/s) 


All  the  foam  data  is 
plotted  from  Rand’s 
1979  report,  SwRI’s 
1999  report  and  the 
Columbia 

investigation;  there  is 
good  agreement  with 
the  theoretical  curve. 
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Final  Estimates  of  Foam  Path 
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Leading  Edge  Structural  Subsystem 


Test  Article 


■Support 

structure 


Pone I  6 
slot 


-Seal 


Panel  7 


Panel  3 


Panel  9 


Panel  10 
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Interior  Cameras  and  Gages 


Up  to  16  high¬ 
speed  video 
cameras  were 
used  per  test  (up 
to  8  outside  of 
the  target,  up  to 
8  inside  the 
target). 

Up  to  250 
channels  of 
strain  gage, 
accelerometer 
and  load  cell 
data  were 
collected. 
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RCC  Panel  6  Test 


•  RCC  Panel  6  had  flown  30  flights  on  Discovery. 

•  The  test  resulted  in  a  cracked  rib. 

•  Damage  thought  insufficient  to  cause  the  loss  of  the  vehicle. 
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RCC  Panel  8 


•  RCC  Panel  8  had  flown  26  missions  on  Atlantis. 
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RCC  Panel  8 


16.0”+/-0.1 


9.6  +/-0.1 


8.3”+/-0.1 


16.75  +\-0.1 


JSC200e46792 
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Experimental  Results  RCC  Panel  6 
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Replicating  RCC  Panel  6  Test 


H 

DB:  rccb01002. pdb  ^ 

Cycle:  3  002  T  .  me :  0.000225046 


Pseudocolor  plot 
Vsr =  SmaA 
DB:  r ccb01 002 . pdb 
Cycle:  1002 
T.we:  0.0002250 

Pc  levels  M0e61 


36. 2629 
31  ■  1  944 

1  ’£  .  1259 
S  I  .0575 
15.9890 
10.9205 


It 


8621 

7836 

2849 


l>axi  3 . 626e*07 
I*  , n:  -4. 285e*06 


X 


z _ fr 


0.225  ms 


user ■cthjdw 

Mon  Ju I  28  22:05:42  2003 


DB:  rccb06411 . pdb 

Cycle:  6413  T . me :  0.00350016 


1.5  ms 


user icthjdw 

Mon  JUI  28  22:23:01  2003 


The  impact  point  0.83”  left  from  the  5-6  T-seal,  18.7”  up  from  the  carrier  panel. 
Flight  direction:  a=5.5°  (bottom  to  top),  P=2.5°  (away  from  center).  J 

The  impact  velocity  was  768  ft/s,  computation  carried  out  to  5  ms. 
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Replicating  RCC  Panel  6  Test 


DB:  rccbl 1501 . pdb 

Cycle:  1  1  Ei2  1  7  ,  me :  0_  002*70008 


S 

M 


Pseudocolor  plot 
Var  s  Smax 
DB:  r ccb 1 1 50 1 . pdb 
Cycle:  1150] 

T  I  we:  0 . 002700 

Pc  levels  f]0e9l 


2. 0000 
H-  1.7500 
_n=l _ 1  -  5000 

H-  1.2500 

H-  1.0000 

-  0.7500 


It 


0. 5000 
0.2500 
0  .  0000 


Max:  1 . 955e*09 

M.n:  -  1 . S42e*08 


Mesh  plot 
Mesh:  rcc_panel 
DB:  r ccbl 1 501 . pdb 
Cyc l e :  11 501 

T i we:  0.002700 


2.7  ms 


3.7  ms 


DB:  rccbl 5743. pdb 
Cyc I e:  15743  7 , me : 

Pseudocolor  plot 
Var :  Smax 
DB:  r  ccb 1 574  3 . pdb 
Cycle:  15743 
T i we:  0. 003700 

Pc  levels  f]0e9i 

2. 0000 
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Experimental  Results  RCC  Panel  8 
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Replicating  RCC  Panel  8  Test 


DB:  rccb00S46.pdb 

Cycle:  54-6  1  .  me :  0.00010007 


Pseudocolor  plot 

Var  :  Smax 

DB:  r ccb00546. pdb 

Cycle:  546 

T i me:  0.0001001 

Pc  levels  M0e6i 


-  40.0075 

■  34.4122_ 

■  28.816 

■  23.221 

■  17.626 

■  12. 

-  6. 


Mesh  p I ot 
Mesh:  rcc 
DB:  r 
Cyc l e : 

T i me:  0 . 


H  DB:  r ccb06046. pdb 
^  Cycle:  6046  T . me :  0.00120007 


Pseudocolor  plot 

Var :  Smax 

DB:  r ccb06046 . pdb 

Cycle:  5046 

T  t  me:  0. 001 200 

Pc  levels  ( ] 0e6 1 


Mesh  plot 
Mesh:  rcc.ttsnel 
DB:  r  ccb06B46 . pdb 
Cyc I e:  604  . 

T  i  me:  0.0  1]  200 


0.1  ms 


1.2  ms 


Mon  Ju I  28  21:45:35  2003 


The  impact  point  7.3”  left  from  the  7-8  T-seal,  25.5”  up  from  the  carrier  panel. 
Flight  direction:  a=5.5°  (bottom  to  top),  P=5.0°  (away  from  center),  30°  clocking. 
The  impact  velocity  was  777  ft/s,  computation  carried  out  to  5  ms. 
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Replicating  RCC  Panel  8  Test 
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Loading  Footprint 
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RCC  panel  6  impact 
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Stresses  RCC  Panels  6  and  8 
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Stress  (ksi) 


Investigation  Conclusions 
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Impact  Velocity  (ft/s) 


The  component  of  velocity 
normal  to  the  impact  surface 
determines  the  local  loading 
pressure. 

Modeling  and  experiments 
showed  that  an  impact  on  the 
underside  of  the  wing  was  not  the 
cause  of  the  accident. 

Forensics,  experiments  and 
modeling  showed  that  the  loss  of 
Columbia  was  due  to  a  foam 
impact  on  RCC  panel  8. 
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The  Physical  Cause 


The  physical  cause  of  the  loss  of  Columbia  and  its  crew  was  a 
breach  in  the  Thermal  Protection  System  on  the  leading  edge  of  the 
left  wing.  The  breach  was  initiated  by  a  piece  of  insulating  foam 
that  separated  from  the  left  bipod  ramp  of  the  External  Tank  and 
struck  the  wing  in  the  vicinity  of  the  lower  half  of  Reinforced 
Carbon-Carbon  panel  8  at  8 1 .9  seconds  after  launch.  During  re¬ 
entry,  this  breach  in  the  Thermal  Protection  System  allowed 
superheated  air  to  penetrate  the  leading-edge  insulation  and 
progressively  melt  the  aluminum  structure  of  the  left  wing, 
resulting  in  a  weakening  of  the  structure  until  increasing 
aerodynamic  forces  caused  loss  of  control,  failure  of  the  wing,  and 
breakup  of  the  Orbiter. 
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SwRI  Three-Pronged  Approach 


Impact  Velocity  (fl/s) 
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CAIB  Recommendations  R3.2-2 

Estimate  Risk 


•  Initiate  a  program  designed  to  increase  the  Orbiter’s 
ability  to  sustain  minor  debris  damage  by  measures  such 
as  improved  impact-resistant  Reinforced  Carbon-Carbon 
and  acreage  tiles.  This  program  should  determine  the 
actual  impact  resistance  of  current  materials  and  the 
effect  of  likely  debris  strikes.  [RTF] 
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CAIB  Recommendation  R3.8-2 
Impact  Damage  Models 


•  Develop,  validate,  and  maintain  physics-based  computer 
models  to  evaluate  Thermal  Protection  System  damage  from 
debris  impacts.  These  tools  should  provide  realistic  and 
timely  estimates  of  any  impact  damage  from  possible  debris 
from  any  source  that  may  ultimately  impact  the  Orbiter. 
Establish  impact  damage  thresholds  that  trigger  responsive 
corrective  action,  such  as  on-orbit  inspection  and  repair, 
when  indicated. 
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Foam  and  Ice  on  the  External  Tank 


LH2  Tank  to  Intertank  Flange 

(BX265) 


LH2  Ice/Frost 
Ramps 

(SLA561/  PPL) 


Intertank  Acreage 

(NCFI24-124) 


Bipod  Closeout 

(BX265/  PDL) 


Umbilical  Ice 


LH2  Tank  Acreage 

(NCFI24-124) 


Bracket/Bellows  Ice 


L02  PAL  Ramp 

(BX250/265/  PDL) 


L02  Tank  to 
Intertank  Flange 

(BX250) 


(2)  L02  /  Intertank 
Ice/Frost  Ramps 

(SLA561/  PDL) 


L02  Tank 
Acreage 

(NCFI24-124) 


LH2  PAL  Ramp 

(BX250/265/  PDL) 
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Physics-Based  Impact  Models: 
Inputs  and  Outputs 


Specific  Impact  Event  Input: 

Impactor  denstity:  pp 

Impactor  Dimensions: 

Length  L,  Width  W,  Height  H 

Impact  Velocity  V  and  Angle  0 

Large  compression  crush-up 
curve  required  (foam) 

Projectile  fracture  (ice) 


Tile  Material  Properties: 

Large  compression  crush-up 
curve  required  (stress  vs.  strain 
in  compression  up  to  90% 
compression) 


Foam,  Ice  and  Ablator 
Impact  Models 


Output:  Computed  Crater  Dimensions 

Depth,  Length,  Width 
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Our  Validation  Triangle 


Experiments 


Large-Scale 

Numerical 

Simulations 


Analytical 

Modeling 


When  experiments, 
large-scale  numerical 
simulations  and  the 


analytical  physics- 
based  model  agree, 
the  physics-based 
model  is  assumed  to 
be  validated. 


^zz  (V) 

PPL 


Both  of  these  models 
are  Physics  Based 


This  is  our  fast-running, 
physics-based  model 
for  flight 
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Powerful  Feature  of  the  Models 


•  For  both  the  fast-running  physics-based  model 
and  CTH 

-  The  model  coding  is  exactly  the  same  for  all 
the  different  cases  that  will  be  shown  -  only 
the  input  material  properties  are  changed. 

•  Thus,  it  is  not  a  different  model  for  ice  or  low 
density  ice,  it  is  not  a  different  model  for  LI-900 
or  FRCI-12;  the  impact  model  is  exactly  the 
same,  only  the  input  material  properties  change. 

•  As  will  be  shown,  excellent  agreement  between 
models  and  test  data  exists  for  all  impact  cases. 
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Thermal  Tiles  (LI-900) 


Silica  foam  (white) 


I'-  _ 


LRSI  White  Tile 


Flown  tile 
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Thermal  Tiles 
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Stress  (10  6Pa) 


Larger  strain  behavior  based  on  Crush-in- 
a-Box  Uniaxial  Strain  Test  Data 


•  Test  performed  with  a  cylinder  from  same  tile  specimen 
confined  with  a  steel  ring. 

•  Data  placed  into  new  tile  model  within  CTH.  Model  includes  a 
yield  at  550  kPa. 
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Material  Properties  for  Impactors 


Ice  Physical  Properties 


•  The  modeling  assumed  the  following  properties  for  ice  (obtained 
from  typical  values  found  in  the  literature): 

-  Modulus  of  Elasticity  (E)  =  8000  MPa 

-  Shear  Modulus  (G)  =  3000  MPa 

-  Density  =  0.914  g/cc  (57  lb/ft3) 

-  Poisson’s  ratio  =  0.33 

-  Flow  stress  =  2.0  MPa 

-  Tensile  strength  =1.0  MPa 

-  Sound  speed  (c0)  =  2954  m/s 
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(DdgOO  SSSJJS 


RT-455  Ablator  Model  Developed  for  CTH 


Tile:  Stress  vs.  Strain  Tile:  Stress  vs.  Strain 


Initial  density  0.66  g/cm3 
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XXSTRESS  (MPa) 


PDL  Model  for  CTH 
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PDL  Stress  vs.  Strain 


Strain(%) 
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Stress  (psi) 


Ice  into  Tile 


Physics-Based  Ice  into  Tile  Model 


Given  the  tile  crush-up  behavior,  the 
model  solves  the  Riemann  problem  to 
find  the  stress  in  the  tile  versus  velocity, 
gzz(v).  A  shear  term  is  also  included. 

Given  this  information,  F  =  ma 
qualitatively  translates  into 

dv  =  czz(v) 

dt  ppL 

(see  the  rest  of  these  charts  for  exactly 
how  the  full  equations  appear). 

The  depth  of  penetration  is  calculated  by 
integrating  this  equation  with  a  numerical 
scheme  to  produce  a  line-of- sight  depth 
of  penetration  PL0S. 
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Tile  Resistance  to  Penetration 


•  The  model  computes  a  penetration  resistance  stress  of  the  form 

7 

^zz  (v)  =  °Hugoniot  (v)  +  -C7crush  \n(a) 

•  Here, 

-  aH  niot(v)  is  the  stress  as  a  function  of  particle  velocity  along 
the  Hugoniot  (the  idealized  one-dimensional  planar  impact). 
This  function  is  computed  based  on  the  large-strain 
compression  curves  and  the  Hugoniot  jump  conditions.  It  is 
stored  as  a  table  for  rapid  look-up  during  the  computation. 

-  ocrush  is  the  crushing  strength  of  the  tile 

-  a(v)  is  the  extent  of  the  deforming  region  within  the  tile,  and  is 
computed  with  a  cavity  expansion  expression.  It  depends  on 
the  material  properties  of  the  tile  and  the  penetration  velocity. 
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Resisting  Stress  (10 3 Psi) 


The  Hugoniot  Contribution  to  the 

Resistance 


0.0  0.4  0.8  1.2  1.6  2.0 

Penetration  Speed  (103ft/s) 


The  model  solves  the  Riemann  problem 
to  find  the  stress  in  the  tile  versus 
velocity,  a Hugonio,(v). 

This  term  is  found  by  simultaneously 
solving  the  Hugoniot  jump  conditions 
using  the  tabular  data  of  tile  properties 

P±  =  l_± 

p  u 

^ Hugoniot  A) 

where  U  is  the  shock  speed  and  u  is  the 
penetration  speed  (U  is  solved  for  also). 

The  plot  shows  the  value  of  the 

°Hugoniot(v)  term  (blue)  as  well  as  the  total 
resisting  stress  including  the  shear 
resistance  term  (red)  in  terms  of 
penetration  velocity  through  the  tile.  > 
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Embedment  Phase 


When  the  projectile  embeds  itself 
into  the  target,  only  part  of  the  face 
of  the  projectile  is  loaded  by  the 
target.  Thus,  the  force  on  the  face  is 
reduced  by  the  appropriate  factor: 

He. 

H 

where 

dv  =  f  °~zz  (v) 

dt  7  ppL 
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Equations  of  Motion 


For  completeness,  we  write  the  full  equations  of  motion,  where  x  is  the  direction 
along  the  surface  of  the  tile  and  y  is  the  normal  direction  into  the  tile: 


<K 

dt 
dv  v 
dt 


-f  —  ^  COS (0)  + 


n 


=-/ 


pPL 

(V)  _ . 


pLHW 


sin(0) 


sin(^)  - 


n 


pLHW 


COS  {6) 


V  =  VVx+Vy 


These  equations  are  integrated  forward  in  time  until  either  the  impactor  speed  v  drops 
below  the  speed  required  to  crush  tile  (i.e.,  it  comes  to  rest)  or  the  impactor  is  forced 
back  to  the  tile  surface  (i.e.,  it  ricochets). 
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Stress  (10 1 2 * * *  6Pa) 


Developing  the  95%  Bounding  Model 


1  2 

1  1 

1  0 
9 
8 
7 

e 

5 
4 
3 

2 

1 

0 

0  20  40  GO  80 


Tile:  Stress  vs+  Strain  (Uniaxial  Strain) 


Test 

With  Yield  M 


We  multiply  each  stress 
term  by  the  scale  factor, 
shifting  the  stress  (hence 
strength  of  the  tile) 
downward) 


Strain  (%) 


Since  the  fast-running  physics-based 
model  is  based  on  material  properties, 
the  only  way  to  adjust  the  model  is  by 
adjusting  those  material  properties. 

We  provide  a  “bounding  model”  by 
adjusting  the  tile  strength  -  that  is, 
making  the  tile  weaker. 

In  particular,  each  point  in  the  data 
curve  is  shifted  down  by  the 
multiplicative  scale  factor,  and  then 
the  Hugoniot  and  shear  terms  are 
recomputed. 
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Confidence  (%) 


Bounding  Model:  95%  Bounding 


Scale  Factor  vs.  Confidence 


Based  on  the  experimental  data,  95%  bounding 
is  achieved  when  only  9  data  points  lie  above  the 
curve. 

With  the  earlier  version  of  the  model,  this  level 
of  bounding  was  achieved  with  a  scale  factor  of 
0.615.  We  have  chosen  to  stick  with  this  scale 
factor.  For  the  updated  model,  at  a  scale  factor 
of  0.615,  7  points  lie  above  the  curve. 

On  page  13  of  Material  Properties  Data,  Volume 
3,  Thermal  Protection  System  Materials  Data, 
the  in-plane  compressive  strength  of  the  LI-900 
tile  material  is  listed: 

-  average  is  70  psi, 

-  90%  minimum  is  53  psi,  corresponding  to  a 
scale  factor  of  0.76 

-  99%  minimum  is  41  psi,  corresponding  to  a 
scale  factor  of  0.59. 

Thus,  the  scale  factors  from  the  material  testing 
are  in  reasonable  agreement  with  that  found 
from  the  ice  impact  tests. 
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Tile  Crater  Depth  (In) 


30° 


F:  L=0.5”  W=  0.5”  H  =0.  5” 


Impact  Speed  (103ft/s) 


Impact  Speed  (10 3 ft /s) 


Depth 


Length 


Z  Position  (cm) 


CTH  calculations  of  Ice  impact  into  LI- 

900  Tile 


•  Work  has  focused  on  3-D 

[3  LI900B 

computations  of  ice  impact 
into  LI-900  tile  using  new 
subroutines  in  CTH. 

•  Approximately  350 
computations  completed  to 
date. 

•  Six  different  projectile 
geometries  were  examined 

•  Velocities  ranged  from  10  to 
300  m/s,  and  impact  angles 
varied  from  10  to  30  degrees. 


r-*j 


ZZ  Stress 


4 


Stress  ZZ 
(d/cm2) 

2.000x10® 


1.241x10® 
9.873x10" 
7.342x10® 
4.810x105 
2.278x10® 
-2.532x10® 
-2.735x10® 
-5.316x10® 
-7.B4Bx10® 
— 1.038x10® 
-1.291x10  s 
—1.544x10® 
— 1.797x10® 


3DR 


-  6 


X  (cm) 


Y  =  5x10“  2 


3-D  CTH  Ice  Impact  Computations 


•  3-D  Computations 

•  2  upright  individual  LI900 
tiles  included  in  the 
computation 

•  4  mm  of  SIP  material 
included 

•  Bottom  support  is  rigid 

•  Cell  size  4  mm  (0.016  in.) 

•  Y=0  is  plane  of  symmetry 


i 


I  lice 
1LI9QOA 

mi  sip 

■  LI900A 


El  SIP 
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CTH  Ablator  Numerical  Simulations 


Ablator  3”  x  0.5”  x  1” 

Tile  6”  x  6”  x  2” 

3-D  calculations 

Y=()  plane  of  symmetry 

Ablator  and  Tile  EOS 
implemented  by  SwRI. 

No  SIP 

1.5  mm  cell  size  (6  cells  across 
ablator  thickness) 

Ran  a  total  of  26  cases 

Two  examples  follow 


ab20deg13  20  degs,  v=  10000. 

GDODAS  G  7/04/0414:34:06  CTHGEN 


0  Time  0.  s 
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30°  Impact  at  350  ft/s  (P2. 11.7-1) 


0.5  ms  1.0  ms  1.5  ms 
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CTH  crush-up  profiles 


Alpha 


T  0  h 

u 

N*  -2 

-4 

-6 

-8 

-1  0 
-2 

3DR 

ab30deg11  30  degs,  v=  10000. 
GBSAEV  7/02/04  20:13:00  CTH 


Alpha 


Alpha  Parameter 


801  T?me=1.0004x10 


-10 
-2 

3DR 

ab30deg11  30  degs,  v=  10000. 
GBSAEV  7/02/04  21:17:05  CTH 


1260  Time=1. 50023x10  3 


DOP  given  by  the  lowest  point  with  alpha  =1.5 
Alpha  is  the  current  density  divided  by  initial  density 


/ 
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Foam  into  Tile 
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BX-265  Foam  Impactors 


A  wide  range  of  projectiles  were  shot: 

-  Return  to  flight: 

•  1.6”xl”xl”;  1  gram  (0.0022  lb) 
.  4”x4”x2”;  20  gram  (0.044  lb) 

•  12”x6”x2”;  90  gram  (0.2  lb) 

-  Columbia  Investigation 

•  22”xl2”x6”;  750  gram  (1.67  lb) 


Loading  Surfaces 


As  the  foam  projectile  impacts  the 
tile,  only  part  of  the  face  of  the 
projectile  is  loaded  by  the  target. 
Thus,  the  force  on  the  face  is  reduced 
by  the  appropriate  factor: 

H 

The  force  on  the  length  of  the 
projectile  acting  upwards  is 


LyW  min(aH  (vy0 )  cos(<9)  +  dcrush_t  sin(<9),  dcrush_t ) 
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Equations  of  Motion 


The  equations  of  motion,  where  x  is  the  direction  along  the  surface  of  the  tile  and  y  is 
the  normal  direction  into  the  tile: 


dt 

civ, 

dt 


-f  —  ^  cos(<9) 


=  -/ 


P,L 


sin(0)  - 


F 


pLHW 


cos  (0) 


V  =  VVX+Vy 


•  These  equations  are  integrated  forward  in  time  until  either  the  impactor  speed  v  drops 
below  the  speed  required  to  crush  tile  (i.e.,  it  comes  to  rest)  or  the  impactor  is  forced 
back  to  the  tile  surface  (i.e.,  it  ricochets). 


(There  is  not  a  sine  term  in  the  first  equation  because  the  compliance  of  the  foam 
leads  to  a  nearly  flat  interface  region  between  the  tile  and  the  foam.) 
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Tile  Resistance  to  Penetration 


•  The  resisting  stress  is  of  the  form 

Czz  (v)  =  °  Hugoniot  (v) 

•  cth  iot(v)  is  the  stress  as  a  function  of  particle  velocity  along  the 
Hugoniot  (the  idealized  one-dimensional  planar  impact).  This 
function  is  computed  based  on  the  large-strain  compression  curves 
and  the  Hugoniot  jump  conditions.  It  is  stored  as  a  table  for  rapid 
look-up  during  the  computation. 
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Particle  Velocity  (m/s) 


Riemann  Solver 


0  100  200  300  400  500 


Impact  Velocity  (m/s) 


•  A  centerpiece  of  the  penetration 
model  is  a  Riemann  solver  that 
calculates  the  penetration  speed  for 
various  impact  speeds,  as  well  as  the 
wave  speeds  in  the  material. 

•  The  solver  solves  the  mass  and 
momentum  forms  of  the  Hugoniot 
jump  conditions  to  compute  the 
stress  and  particle  velocity  at  the 
interface  between  the  two  colliding 
materials 

•  The  solver  uses  exactly  the  same 
equation-of-state  tables  used  by  CTH 
for  modeling  the  tile  and  foam. 

•  The  Riemann  solver  solves  2-,  3-  and 
4-wave  problems,  all  of  which  can 
arise  in  foam  vs.  tile  impacts, 
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Damage-No  Damage  Transition 


Explicit  expression  for  the  transition  velocity  for  a  90°  impact  is 


V crush  ~  Uet  +  Uef  +  < 


1 


1 


\ 


V  Pef  Pf  crush-t  )  J 


crush— t  ^ crush- f  ) 


1/2 


If  there  were  no  edge  catching,  then  the  transition  would  be  given 
by 


V. 


v  = 


crush 


sin(#) 
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Edge  Catch 


The  tile  is  pushed  down  by  the  foam,  which  has  inertial  resistance  as  well  as 
a  spring-like  behavior  of  the  SIP  (represented  by  a  spring  constant  k ): 


y  tile 


Lywa, -M, 


The  current  critical  (transition)  velocity  is  then  computed  according  to 


v. 


y  = 


crush 


sl^me 1  T SIP  sinO  /  4)  +  (1  -  VA v,,/e  /  Ts,p ) sin  6 


If  the  current  velocity  is  above  this  value,  the  edge  catches  and  penetration 
into  the  tile  begins. 
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Tile  Fracture  Algorithm 


Bottom  of  Tile 


There  is  a  tile  fracture  algorithm,  in 
an  effort  to  determine  when  a  tile 
breaks,  as  is  seen  in  the  shear-out  tile 
failure  mode,  for  example. 

The  model  is  still  under  examination. 

A  cross  section  is  taken  through  the 
tile,  with  area  A.  If  the  area  loaded 
by  the  impactor  is  given  by  Ah,  then 
the  tile  breaks  if  the  force  exerted  by 
the  projectile  is  greater  than  the  load 
the  tile  ligament  can  support  (we 
assume  that  the  tensile  strength  of  the 
tile  is  equal  to  c>crush_,): 

AL<y{v)  >  (A-  Al)  <7cmsh_t 
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Tile  Crater  Depth  (In) 


Foam  Model  Examples 


0  . 


<1) 

O 


o 


Impact  Speed  (10 3 ft /s) 


Impact  Speed  (10 3 ft /s) 
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Test  1.1. 9-2,  30  degs,  220  m/s  (718  fps) 
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Tile- Out  or  Pop- Off 
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Tile  Crush-Up  at  Bottom  of  Tile 


•  In  3d  CTH  computations,  we  do  not  quantitatively  see 
as  much  crush-up  at  the  bottom  layer  as  we  see  in  one¬ 
dimensional  computations  (to  be  shown  later)  (i.e.,  a 
is  less):  presumably  this  is  due  to  the  larger  cell  size  in 
3D. 

•  However,  we  do  see  crush-up  all  along  the  bottom  of 
the  tile. 

•  For  example  (next  three  charts)... 
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BX-265  vs.  Tile  CTH  Computations 


•  3-D  Computations 

•  9  upright  individual  LI-900 
tiles  included  in  the 
computation 

•  SIP  simulated  as  a  4  mm 
(0.016  in.)  gap 

•  Bottom  support  is  rigid 

•  Cell  size  4  mm  (0.016  in.) 

•  Y=0  is  plane  of  symmetry 


LO 


in 

OJ 

I 


o 

LO 

I 


Z  (cm) 


Simulation  of  PI.  1.2-4:  V=957  ft/s,  0  =  20° 
Crush-up  plots  in  plane  of  symmetry 
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Simulation  of  PI.  1.2-4:  V=957  ft/s,  0  =  20° 
Crash-up  plots  in  planes  perpendicular  to  impact 

(High  resolution) 


3 

2 

1 

0 

-1 


rvi 

-4 

-5 

-6 

-7 

3Df? 

3d20,ceg5'G 

^WKCWQ 

3 

2 

1 

0 

1  -2 
M 

-4 

5 

-7 

m 

3d20ceg50 

FWKDUQ 


0  2  4  0  8  1  024  6  S 


X  (cm)  Xs-UxlO1 

re  Fined  two  tiles  2D  degs,  v=  30000. 

6/24/04  04:00:31  CTH  2462  T'me=2,Ofl0&KH3" J 

Pressure 


30P  X  (cm)  X--' 

3d2Ocieg50  'efinEid  two  tiles  20  deg&.  v-  30000,  SOP 

FWKDUO  6/24/04  01:09:41  CTH  19BB  Pma-Z.OOQ  3d™—' ™ 

rv 

Pressure 

i  3  - 1 - 1 - 1 - 1 - 1 - j - f - 1 - 


Pressure 


0  2  4  6  8 

X  (cm)  X=- 

yifM  l  HI*.*  in  rJrn^r-  ,  r-  IfWinfl 


Pressure 

3  | - : - 1 - 1 - 1 - 1 - 1 - 1 - r 


0  2  4  6  6 


0  2  4  8  6  1 

3Df?  X  (cm)  3W  X  (cm) 

*  X--6.  3<j2[>C9g5o  refmed  two  Hies  20  degs,  v=  30000.  3d20d»g50  rafined  two  Hies.  20  degs.  vs  30000. 

refined  two  tiles  20  degs.  v=  30000.  ^WKOUQ  6/24/04  01:09:41  CTH  1985  Time=2.O0  FWKD0G  6/24/0*  04:00:31  CTN  2432  Tinted 

6/24/04  04:00:31  CTH  2452  Times  2.0898*10^ 


Pressure 

M$4ib  Suaiiwn. 
iWH^ 

le.’SjH1 

:ac.Bl? 

UMtaltf* 
'itHirt1 

Wh#* 

iMfr.ia* 
US  BMW6 
US?*W6 

0  2  4  6  8  10 

3Df?  X  (cm)  X'~B. 

3d20'ceg50  reFined  two  Hies  20  degs.  v-  30000. 


■?33^5 

-JnCutt" 

sJaWO* 

\*awt' 

UBiW6 

US7>«* 

U9S*#1 

1134 

u}T£.ra* 

UtnirtpS 

10 


3DR  X  (cm)  X=0. 

3d20t!eg50  Fined  two  tiles  20  degs,  vs  30000. 

FWKDUQ  6/24/04  01:09:41  CTH  19B3  Tlme-2.O0Dt4)!lO''‘ 


104 


Simulation  of  PI.  1.2-4:  V=957  ft/s,  0  =  20° 
Crash-up  plots  in  planes  perpendicular  to  impact 

(High  resolution) 


Yellow  at 
bottom  is 
crashed  tile 
material 


Pressure 

Alpha  Parameter  — 

11.994x10° 
1.933x10° 
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|  1.011x10° 

0  2  4  6  8  10 

3DR  X  (cm)  X=0. 

3d20deg50  refined  two  tiles  20  degs,  v=  30000. 

FWKDUQ  6/24/04  01:09:41  CTH  1988  Time=2.00014x10“5 
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RCC  Impact  Testing 
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CAIB  Recommendation  R3.3-4 
RCC  Properties 


•  In  order  to  understand  the  true  material 
characteristics  of  Reinforced  Carbon-Carbon 
components,  develop  a  comprehensive 
database  of  flown  Reinforced  Carbon-Carbon 
material  characteristics  by  destructive  testing 
and  evaluation. 
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A  ram  is  System  Allows  Direct 
Measurement  of  Dynamic  Displacements 


Displacement  Z  at  Max  Displacement 


[in] 

1  0.136 
0.120 

0.100 

0.080 

0.060 

0.040 

0.020 
H  0.000 
[J-0.023 


- 


Strains  can  then  be  computed, 
and  both  displacements  and 
strains  can  be  compared  to 
DYNA  Team  computations. 
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Return  to  Flight  by  Discovery 


STS-114 


•  1 14th  flight  of  the  Space  Shuttle 
Program 

•  3 1 st  flight  of  Discovery 

•  Discovery  first  flew  in  August,  1984 
(it  is  now  the  oldest  operational 
shuttle,  as  well  as  the  most  flown) 

•  Mission:  Tuesday,  July  26  -  Tuesday, 
August  9,  2005 

•  Crewed  by 

-  Eileen  Collins 

-  Jim  Kelly 

-  Charlie  Camarda 

-  Wendy  Lawrence 

-  Steve  Robinson 

-  Andy  Thomas 

-  Soichi  Noguchi 


Prelaunch  Ice  Formation 


111 


•  The  following  plots  from  the  ice  model  are  for  four  masses: 
-  Set  Length  (1.5”  or  2”) 

•  L=1.5”,  H=W=0.767”,  mass  =  13.2  gram  =  0.029  lb 

•  L=2.0”,  H=W=1.02”,  mass  =  3 1 .2  gram  =  0.069  lb 


-  Length/Width  ratio  =  2.5 

•  L=1.77’\  H=W=0.707”,  mass  =  13.2  gram  =  0.029  lb 

•  L=2.35”,  H=W=0.94”,  mass  =  31.1  gram  =  0.069  lb 

•  Impact  speed  in  plots  goes  from  0  to  200  ft/s.  Understand  that  the  model  has  only 
been  baselined  for  speeds  of  80  ft/s  to  1500  ft/s.  The  concern  is  that  ice  may  not 
fracture  at  the  lower  velocities,  and  so  there  may  be  larger  depths  of  penetration  than 
the  model  predicts  (the  model  assumes  the  ice  fractures). 

•  Impact  angles  range  from  5°  to  55°  in  5°  increments. 


•  Blue  curves  are  the  nominal  model,  red  curves  are  the  95%  bounding  curves. 

•  (If  you  want,  you  can  pull  the  picture  out  of  the  presentation  and  enlarge  it.) 
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Tili  Oofer  Depth  fin) 


L=1.5”  W=  0.767”  H  =0.767”  5° 
Mass  =13.2  gram  =  0.029  lb 
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Tile  Croter  Depth  (in) 


L=2.0”  W=  1.02”  H  =1.02”  5°  -  55° 
Mass  =  31.2  gram  =  0.069  lb 
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Tom  Thermal  Blanket 
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Impact  of  Thermal  Blanket  on  Speed 

Brake 


We  performed  computations  to 
examine  the  impact  of  the  thermal 
blanket  were  it  to  tear  free  during  re¬ 
entry  for  STS- 1 14  and  strike  the 
speed  brake. 

Computations  were  performed  with 
the  hydrocode  CTH. 

The  impact  computations  were 
performed  for  a  velocity  of  1600  ft/s. 

“Conservative”  issues  related  to  the 
modeling  were  explored  with  LS- 
DYNA. 
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Ribbon  Fabric  Projectile:  12.7”  x  3.8”  x  0.18”  (.0531b) 
impacting  thermal  blanket  and  speed  brake  at  1600  ft/s 
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□  Mat  4 
U Mat  5 

□  Mat  6 


ruOl  Nextel  12.7  in*  vs,  blanket/honeycomg  ,  v=1SGQ  ff/s 
HCUAAW  G  8/03/05  20:00:22  CTHGEN  0  Time  0.  s 


Ribbon  Fabric  Projectile:  12.7”  x  3.8”  x  0.18”  (.0531b) 
impacting  thermal  blanket  and  speed  brake  at  1600  ft/s 
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□  Mat  6 


ruOl  Mextel  12.7  in.  vs,  blanket/honeycomg  ,  v=1600  ft/s 
HCUAEA  8/03/05  20:48:46  CTH  261  Time  1.253x10  4  s 


Ribbon  Fabric  Projectile:  12.7”  x  3.8”  x  0.18”  (.0531b) 
impacting  thermal  blanket  and  speed  brake  at  1600  ft/s 
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□  Mat  1 
■  Mol  2 
BMat  3 

□  Mat  4 
UMat  5 

□  Mat  6 


ruGI  Nextel  12 J  in.  vs,  blanket/honeycomg  ,  v=1600  ft/s 
HCUAEA  3/03/05  21:25:07  CTH  431  Time  2.5024x10  4  s 


Ribbon  Fabric  Projectile:  12.7”  x  3.8”  x  0.18”  (.0531b) 
impacting  thermal  blanket  and  speed  brake  at  1600  ft/s 
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□  Mat  1 
■  Mai  2 
BMcrt  3 

□  Mat  4 
Cl  Mat  5 

□  Mat  6 


ruGI  Nextel  12,7  in.  vs.  blanket/honeycomg  ,  v=1600  ff/s 
HCUAEA  8/03/05  21:59:46  CTH  692  Time  3.7543x10  4  s 


Ribbon  Fabric  Projectile:  12.7”  x  3.8”  x  0.18”  (.0531b) 
impacting  thermal  blanket  and  speed  brake  at  1600  ft/s 
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□  Mat  1 
■  Mol  2 
BMut  3 

□  Mat  4 
H] Mat  5 

□  Mat  6 


(Last) 


ruOl  Mextel  12,7  in.  vs,  blanket/honeycomg  ,  v=1600  ft/s 
HCUAEA  3/03/05  22:34:28  CTH  892  Time  5.0023x10  4  $ 
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Initial  Geometry 
Impact  speed:  1600  ft/s 


0.053  lb/ 24  g  0.0177  lb/ 8  g  0.0106  lb / 4.8  g  0.0053  lb /  2.4  g 
12.7”x3.8”x0.18”  4.23”x3.8”x0.18”  2.54”x3.8”x0.18”  1.27”x3.8”x0.18” 
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Post  Impact  Geometry 
Impact  speed:  1600  ft/s 


0.053  lb/ 24  g  0.0177  lb/ 8  g  0.0106 lb / 4.8  g 
12.7”x3.8”x0.18”  4.23”x3.8”x0.18”  2.54”x3.8”x0.18” 


0.0053  lb  /  2.4  g 
1.27”x3.8”x0.18” 
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Examining  an  “Conservatism”  in  Our 
Honeycomb  Modeling 


•  Modeling  the  aluminum  honeycomb  is  difficult. 

•  In  the  Eulerian  code  CTH,  we  are  unable  to  resolve  the  structure  of 
the  front  and  back  sheets  of  the  honeycomb,  since  our  computational 
cell  size  is  2.5  mm  and  the  thickness  of  the  aluminum  sheet  is  0.1 1” 
(0.3  mm). 

•  To  see  if  the  honeycomb  fails  much  more  easily  than  it  should,  a 
series  of  computations  were  performed  in  DYNA  to  explore  with  a 
honeycomb  we  built  out  of  shell  elements. 
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Conclusions 


Thermal  blanket  impact  into  the  speed 
brake  was  modeled  using  current  impact 
tools. 

It  was  assumed  the  fabric  behaved 
similarly  to  LRSI  LI-900  tiles  (the 
density  is  similar). 

The  blanket  on  the  outside  of  the  speed 
brake  was  modeled  as  an  HRSI  LI-900 
tile. 

The  honeycomb  skin  of  the  speed  break 
was  modeled  with  newly-created 
homogenized  honeycomb  material 
properties  placed  in  the  foam  model. 

All  impacts  modeled  in  CTH  (at  1600 
ft/s)  resulted  in  penetration  of  the  tile 
and  honeycomb. 

LS-DYNA  honeycomb  model  contained 
a  0.0135  (6.1  g)  impact  (0.053/4).  v 
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Post  Flight  Inspection 


126 


127 


Impact  Damage  to  Tile 


I 


Proj.  Mass/Width  (10  Ib/in) 


Charts  to  Estimate  Debris  Environment 

from  Ice  Damage 


1  .  0 
0 . 9 
0 . 8 
0  *  7 
0  .  S 
0 . 5 
0.4 
0 . 3 
0 . 2 
0  .  1 
0 . 0 

0.0  0.2  0.4  0.6  0.8  1.0  1.2 

Velocity  (lO^ft/s) 


Crater  Length=3.Q00  in,  Angle=2Q  degs 

r 


-  Proj  Length=  0.250 

Proj  Length=  1.500 
Pro|  Length=  2.750 

-  Proj  Length=  4.000 

Pro|  Length=  5.250 

— - Proj  Heigth=  0.250 

Pro|  He1gth=  0.500 
Proj  H«igth=  0.750 

- Pro|  He1gth=  1.000 

Proj  Heigth=  1.750 


Proj.  Mass  (10  lb) 


Charts  to  Estimate  Debris  Environment 

from  Foam  Damage 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


.  0 
.  9 
.  8 

.  7 
.  S 

*  5 
.  4 
.  3 
.  2 

*  1 
.  0 

0.0  0.5  1.0  1.5  2.0  2.5  3.0 

Velocity  (103ft/s) 


Crater  Length^  3,00  Angle=20.00  W=1.0 


CAIB  Recommendations  R3.2-l,2 
Reduce  Damage  and  Risk 


•  Initiate  an  aggressive  program  to  eliminate  all  External 
Tank  Thermal  Protection  System  debris-shedding  at  the 
source  with  particular  emphasis  on  the  region  where  the 
bipod  struts  attach  to  the  External  Tank.  [RTF] 

•  Initiate  a  program  designed  to  increase  the  Orbiter’s 
ability  to  sustain  minor  debris  damage  by  measures  such 
as  improved  impact-resistant  Reinforced  Carbon-Carbon 
and  acreage  tiles.  This  program  should  determine  the 
actual  impact  resistance  of  current  materials  and  the 
effect  of  likely  debris  strikes.  [RTF] 
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STS-121 


•  1 15th  flight  of  the  Space  Shuttle 
Program 

•  32nd  flight  of  Discovery 

•  Discovery  first  flew  in  1984 

•  Projected  Launch:  May  2006 

•  Crewed  by 

-  Steve  Lindsey 

-  Mark  Kelly 

-  Michael  Fossum 

-  Lisa  M.  Nowak 

-  Piers  Sellers 

-  Stephanie  Wilson 

-  Thomas  Reiter 
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Our  Approach  to  Understanding 
and  Our  Validation  Triangle 


Experiments 


Large-Scale 

Numerical 

Simulations 


Analytical 

Modeling 


When  experiments, 
large-scale  numerical 
simulations  and  the 


analytical  physics- 
based  model  agree, 
the  physics-based 
model  is  assumed  to 
be  validated. 


^zz  (y) 


Both  of  these  models 
are  Physics  Based 


This  is  our  fast-running, 
physics-based  model 
for  flight 
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End  of  Charts 
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Introduction  -  Challenge  for  Ceramic  Armor 

Perspective  on  Damage  Diagnostics  &  Cognitive 
Visualization 

Advanced  3D  Voxel  Analysis  &  Visualization 
3D  XCT  Damage  Characterization  &  Visualization 
Summary  Comments 
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•  History:  Application  of 
ceramic  armor  against 
high  L/D  penetrators  is 
in  its’  third  millennium. 


Ancient  Chinese  terra 
cotta  armor  vest 


•  Still  Searching  for 
Best  Ceramic  Armor ! 
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Diagnostics  &  Performance 


Penetration  Analysis: 

-  DOP,  V50,  Field  Ballistic  Tests 


Ballistic 

Impact 


Penetration 

Analysis 


Sj 


i 


Damage 

Diagnostics 


r  m 
m 


Ceramic 

erformance 


DESIGN  BASIS: 

•  Theoretical 

•  Empirical  & 

Numerical 

Computational  Focus 

•  Diagnostic/  Analytical 
&  Mechanistic  Focus 


Damage  Diagnostics  &  Assessment: 

-  Destructive  Sectioning  &  2D  Examination 

-  Traditional  Nondestructive  Examination 


-  High  Resolution  X-ray  Computed 
Tomography,  XCT,  for  3D  Diagnosis 

Ideally  we  want  a  Engineering  Predictive  Modeling  Capability 
addressing  both  penetration  &  damage  considerations. 
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Define  Problem  (Challenge) 

-  Create  Engineering  Approach  (Plan 
Data  (Acquire  &  Process) 

-  Information  (Analyze) 

-  Knowledge  (Understanding) 


“Imagination  is  more  important  than  knowledge. 
Knowledge  is  limited.  Imagination  encircles  the  world. 
Albert  Einstein 

IMAGINATION 


Creativity  (New  Ideas) 

-  innovation  (Putting  Ideas  to  Work) 

-  Applications  (Utilization  of  Technology) 
-  Presentation  &  Reporting 


DATA 


INFORMATION  KNOWLEDGE 
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•Sophisticated  image  analysis  and 
visualization  capability  to  process,  analyze 
and  visualize  voxel/volume  data. 

•Up  to  3  GB  of  memory  utilization  with 
Windows  XP  Professional  OS 

•Multiple  Import/Export  File  Formats 
•Virtual  Metrology  Capabilities 
•Variable  Transparency  &  Virtual  Sectioning 
•Iso-Surface  Extraction 
•Segmentation  &  Grey-Value-Classification 
•Porosity  /  Defect  Analysis 
•Wall  Thickness  Analysis 
•Stereo  Viewing  Tool 

www.volumegraphics.com 
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StudioMax  vl  .2.1 


Encapsulated  TiBo  Experiment  (N.L.  Rupert,  ARL  -1997) 

•Single  Shot  (Full  Penetration  w/o  compressive  ring) 

•Single  Shot  (Partial  Penetration  with  compressive  prestress 
ring) 

•Double  Shot  (Full  Penetration  with  compressive  prestress  ring) 

Summary  Damage  Observations 

Penetration  Decrease  with  Prestress  (17-4  PH  Ring) 

Surface  Topography  -  Ring  Steps,  Radial  Expansion  &  Cracking 
W-alloy  Residual  Fragments 
Complex  Cracking  Modes 
Impact  Induced  Porosity 
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Macro  -Photos 

^■~9. 

Impact  surface 


Single  Shot  Single  Shot  Double  Shot 

(Full  Penetration  (Partial  Penetration  with  (Full  Penetration  with 


compressive  prestress  compressive  prestress  ring) 
ring) 


w/o  compressive 
ring) ,  — 


XCT  Scans  ~ 
mid-thickness 
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(z  =  015)  (z  =  020)  (z  =  025)  (z  =  030 


£»• 


(z  =  035) 


\/ 


(z  =  040) 


/ 


v 


(z  =  045) 


X 

\ 

I*:  "  \ 


¥ 


n 


\ 


:jA- 


/ 

/ 


Front 

Face 


(z  =  050) 


(z  =  055) 


Sequential  XCT  Scans  showing  impact  damage  cracking  features  and  residual 
penetrator  (white)  in  TiB2  SI  wo  Disk  -  near  back  (Z=015)  to  front  face  (Z=55). 
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J -T J I SjziDi  511 

dlnaram  Joel 


OD  Radial 
cracks  on 
Impact 
Surface 


Macro-photograph  -  Normal  View 
Surface  Steps  -  NOT  Visable 
Radial  OD  Cracks  -  ARE  Visable 


XCT  3D  Solid  Object  -  Oblique  View 
Surface  Steps  -  ARE  Visable 
Radial  OD  Cracks  -  ARE  Visable 


Associa^ 
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Pi-1 

Ttl 

m 

Penetrator  &  Ceramic  Rubble 


I  C:/Dqcu  ments  and  Settin  gs/J  oe  Wells/D  eskto  p/XCT  0  vera 


Surface  Steps  -Sectioned 


TiB2  Slwo 


Step  Height 
C  =  ~3.5  mm 
B  =  ~2.5  mm 
A  =  ~1.2  mm 


I 


Step  Heights  Vary  along 
the  ring  circumference 


Note:  Higher  density  (lighter  color)  of  Steps  vs  Bulk  TiB2 
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'UN 


TiB2  Slwo 


l^xiai  ix.y)  Slice  Si 


Axial  Slice  #51 


Dia.  A  =  73.8  mm 
Dia.  B  =  73.4  mm 
Dia.  C  =  72.3  mm 
Dia.  D  =  72.4  mm 


Axial  Slice  #41 


Dia.  A  =  72.7  mm 
Dia.  B  =  72.9  mm 
Dia.  C  =  72.1  mm 
Dia.  D  =  72.0  mm 
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Nonuniform  -  but  localized  radial  expansion  near  impact  surface 


#51 

#41 

L 


Opaque  3D  Solid 
Object  Reconstructions 


TiB2  S2w 


Rings 


Segmented  &  Variable 
Transparency 

Fragments  are  Porous 


Virtual  Metrology 

W  =  -22  mm  H  =  -24.5  mm  W  =  -60  mm  H  =  -25.5  mm 

As  =  4794  mm2  V  =  2076  mm3  As  =  891 1  mm2  V  =  2859  mm3 
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Damage 


Early  Point  Cloud 
I  mages  of  3D  Cracking 

impact  face  — — — ; — — 


TiB2  IS  w/o 
Prestress  Ring 


Concentric5 
Ring  Cracks 


rear  face 


Schematic  of  Concentric 


Hourglass  Ring  Cracking 
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Visualization  of  3D  Cracking 

Damage  in  TiB2 


Orthogonal  Sectioning 
3D  Oblique  View 


Recent  (Preliminary)  3D  I  mages 
of  Segmented  I  mpact  Cracking 


3D  Oblique  View 


Frontal  View 
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count 


Defect  size  histogram 


Porosity  along  Ring  Cracks 


TiB2  IS  w/o 
Prestress  Ring 


mm3 
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Point  Cloud  Visualizations  of  Spiral  Crackin 


in  Ballistic  Impact  Samples 


ivifrain 


Orbital  Meso-Cracking 


APM2  Sample#2 


Projectile  Direction 


Ti-6A1-4V  pc  showing  spiral  cracking 
(Full  Penetration) 


impact  face 


TiC  pc  showing  spiral-  blue 


TiB2  surfaced  pc  showing  spiral 


(No  Penetration) 


(Dual  Impact  -  Full  Penetration) 
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•  The  NDE  Diagnostic  Interrogation  of  Impact 
Damage  in  Armor  Ceramics  is  a  Challenging  Task. 

•  XCT  Diagnostics,  Voxel  Analysis,  and  3D 
Visualization  have  revealed  new  details  &  insights 
into: 

-  Impact  Surface  -  Topography  &  Damage 

-  Internal  Residual  Fragment  Distribution 

-  Internal  Mesoscale  Cracking  Modes 

-  Impact-created  Porosity/Void  Distributions 

-  Volumetric  (3D)  Damage  Perspectives 

•  The  XCT  Diagnostic  approach  to  armor  ceramic 
Damage  Analysis  &  Visualization  is  NOT  yet  widely 
practiced. 

•  Further  Improvements  in  and  Benefits  from  this 
technique  are  possible  and  realistically  anticipated. 
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"The  only 
sure  way  to 
avoid  making 
mistakes  is 
to  have  m 
mm  ideas/8 


for 

Listening  Folks! 
Any  Questions? 


3Tm  Wf? 


As  far  as  the  laws  of  mathematics  refer  to  reality,  they  are  not  certain: 
and  as  far  as  they  are  certain,  they  do  not  refer  to  reality”  -  A.  Einstein 
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Research 
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ight  Characteristics  of 
iectile  With  T rianaular 


oss- 


Dr.  Wenjun  Yi, 


Pr£f.  J  ianping  Qian 


hg_Zhang, 


Ballistic  Research  Laboratory  of  China 


J9Ut 


INTRODUCTION 


EXPERIMENT  RESEARCH 

ANALYSIS  OF  FLIGHT  PEP 

) 

CONCLUSIONS 


NTR0DUCTION 


Recently,  the  projectile  with  non-circular  cross-section  has 


L  A 

. 

L-CU 


/ 


paid  more  and  more  attention,  for  example  guidance  projectile, 
the  changing  of  the  shell  shape  from  common  circular  to  non¬ 


circular  cross-section  makes  it  possible  to  improve  projectiles 
storage,  transport,  discarding  and  aerodynamic  characteristics. 
The  experiments  prove  that  the  non-circular  cross-section  has 


better  rigidity  in  the  same  area  and  lighter  in  the  same  length. 
This  paper  takes  armor-piercing  projectile  for  example, 
contrasting  and  researching  the  aerodynamic  characteristics 
and  trajectory  characteristics  of  triangular  cross-section 
projectile  and  circular  cross-section  projectile  respectively. 


The  structure  of  triangular  cross-section  projectile 


Triangular  cross-section  projectiles 

Armor  Piercing  Fin  Stabilized  Discarding  Sabot 


The  structure  of  circular  cross-section  projectile 


.  g.,  ■  1  mmamm 

liwi 

^Aj 

V( 

EXPERIMENT  RESEARCH 


Experiment  show 


EXPERIMENT  RESEARCH 


Shadowgraphic  image  of  triangular  cross-section  projectile 


EXPERIMENT  RESEARCH 


Shadowgraphic  image  of  circular  cross-section  projectile 


Yaw/deg, 


Yaw  angle  curve  of  triangular  cross-section  projectile 


Yaw  angle  curve  of  circular  cross-section  projectile 


EXPE 
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: RI  M  ENT 


RESEA 

V  / 


EXPERIMENT  RESEARCH 


Table  1.  The  result  of  the  experiment 


item 

No.l 

No  .2 

No  .3 

No  .4 

No  .5 

No  .6 

Shape 

Circular 

Triangular 

Circular 

Triangular 

Circular 

Triangular 

Mach 

1.8514 

1.9302 

1.8996 

1.9101 

1.9034 

1.8954 

*4 

2.472 

2.463 

2.397 

2.517 

2.413 

2.529 

0.0025 

0.0035 

0.0024 

0.0034 

0.0024 

0.0034 

*4 

0.0024 

0.0034 

0.0024 

0.0035 

0.0025 

0.0034 

vo 

648.73 

675.14 

665.20 

670.20 

666.93 

663.76 

1 

62.08 

55.71 

61.79 

54.94 

61.44 

54.17 

6.17 

7.43 

6.46 

6.87 

6.03 

6.62 

ANALYSIS  OF^WGHT  PERFORMANCE 

According  to  the  aerodynamic  coefficients  of  the 
projectile  gained  by  experiments,  contrasting  and 
analyzing  the  triangular  cross-section  projectile 
and  circular  cross-section  projectile  in  fight 
performance  such  as  resistance  characteristic, 
stability  under  low  speed  rotation, 
maneuverability  and  so  on. 


ANALYSIS  OF^LIGHT  PERFORMANCE 

#  \  /  / 

/< —  \  V  / 

T able  2  shows  the  analysis  result  of  triangular  cross- 
section  projectile  and  circular  cross-section 
projectile  in  static  stability. 


Table  2.  Static  stability  contrasts 


Content 

s 

m; 

/<te> 

Tims) 

Mm) 

Triangul 

ar 

-3.612 

5.528 

0.6653 

75.49 

83.23 

55.71 

Circular 

-2.813 

4.622 

0.6086 

69.37 

90.57 

60.13 

(Notes: 


m„ 


m„ 


—  -  —  x  1 00% »  stand  for  static  stability  allowance) 
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#  \  /  / 

\  V  / 

The  analysis  and  results  shows  that  static  stability 
allowance  of  triangular  cross-section  projectile  is  higher 
than  circular  cross-section  projectile,  its  stabilizing 
moment  and  dumping  moment  is  bigger  than  circular 
cross-section  projectile,  especial  the  stability  moment. 

To  triangular  cross-section  projectile,  the  oscillatory 
frequency  is  higher  and  the  oscillatory  wavelength  is 
short.  These  characteristics  represent  that  the  stability  of 
triangular  cross-section  projectile  is  better. 


ANALYSIS  OF^LIGHT  PERFORMANCE 

/  X  #  / 

\  V  / 

Table  3  shows  the  analysis  result  of  triangular  cross- 
section  projectile  and  circular  cross-section  projectile 
in  stability  characteristic  under  low  speed  rotation. 


Table  2 

i.  Stability  charact 

eristic  contra 

sts 

Content 

l^{kg.m2) 

yy}^ 

w  lT 

m  X 

Triangular  0.0986  xl0~5  0.00343  -0.003: 


Circular 


0.0979  X  10 


0.00245 


-0.0032 
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The  analysis  and  results  shows:  the  Magnus  moment 
of  triangular  cross-section  projectile  is  higher  than 
circular  cross-section  projectile.  To  the  projectile  in 
flight,  the  Magnus  moment  is  disturbance  moment. 
The  bigger  of  Magnus  moment,  the  bigger  interaction 
of  pitching  and  yaw  by  rolling,  maybe  it  will  lead  to 
great  distribution.  So  the  effect  of  reduce  offcenter  and 
initial  disturbance  of  triangular  cross-section  projectile 
by  low  speed  rotation  is  not  as  well  as  circular  cross- 
section  projectile. 
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Table  4  shows  the  triangular  cross-section  projectile  and  circular 
cross-section  projectile  in  maneuverability.  The  analysis  and  results 
shows:  the  drag  coefficient  of  triangular  cross-section  projectile  is 
lower  than  circular  cross-section  projectile’s,  the  triangular  cross- 
section  projectile  has  higher  lift-drag  ratio,  and  could  supply  bigger 

normal  overload,  consequently  it  has  better  maneuverability. 


Table  4.  Maneuverability  contrasts 

5 
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Content 

CN 

(L  / 

Triangular 

0.366 

0.  3093 

2.197 

Circular 


0.2686 


1.901 
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■The  projectile  with  triangular  cross-sectio 
smaller  flight  resistance,  less  kinetic  energy 
higher  ratio  of  kinetic  energy  cross  section 
when  impacting. 

/  /  \ 

■The  projectile  with  triangular  cross-section  has 

better  static  stability  and  it’s  static  stability  allowance 
is  bigger  than  the  circular  cross-section  projectile’s. 
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Under  low  speed  rotation,  the  projectile  with 
triangular  cross-sections  has  great  Magnus  moment 
and  more  possible  to  increase  the  projectile 
dispersion. 

❖The  projectile  with  triangular  cross-section  has 
higher  lift-drag  ratio,  could  supply  bigger  normal 
overload  and  better  maneuverability.  Therefore  it’s 
more  suitably  used  for  projectile  missile  shape. 
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